|
|
||||||||
J Thorac Cardiovasc Surg 2008;136:1528-1535
© 2008 The American Association for Thoracic Surgery
Acquired Cardiovascular Disease |
Department of Human and Engineered Environmental Studies, Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa-shi, Chiba, Japan
Received for publication February 25, 2008; revisions received April 22, 2008; accepted for publication May 19, 2008. * Address for reprints: Seiryo Sugiura, MD, PhD, Department of Human and Engineered Environmental Studies, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwa-no-ha, Kashiwa-shi, Chiba, Japan. (Email: Sugiura{at}k.u-tokyo.ac.jp).
| Abstract |
|---|
|
|
|---|
Methods: The proximal portions of the ascending aorta and aortic valves were modeled with blood flowing inside. Blood flow and the motion of aortic valve leaflets were studied while applying a physiologic pressure waveform using fluid–structure interaction finite element analysis. Waveforms were varied to simulate the change in cardiac contractility.
Results: In the aorta without the sinus, the time during which the valve was open was longer and the rapid valve closing velocity was faster under all conditions studied. In the pseudosinus model, we could clearly observe vortex formation from the early phase of ejection, which seemed to facilitate the gradual but smooth closure of the valve. Valve leaflets without the sinus were subject to greater stress and underwent bending deformation in the longitudinal direction.
Conclusions: Sinuses of Valsalva facilitate the smooth closure of the aortic valve, thereby avoiding the building up of abnormal stress in the leaflet. Such an effect may assure the durability of valve leaflets in aortic grafts with a pseudosinus.
| Introduction |
|---|
|
|
|---|
|
|
Sinuses of Valsalva, with their characteristic morphologic features, have attracted the interest of researchers, and earlier modeling studies have suggested that the sinuses function not only to prevent the contact of valve leaflet with the aortic wall, but also to facilitate valve closure by the formation of vortices inside them.1,2
The emergence of valve-sparing aortic root replacement procedures for the treatment of patients with aortic root disease has provided us with the unique opportunity to test such hypotheses. Leyh and colleagues3
compared patients who had undergone tube replacement of the aortic root (reimplantation [David I] procedure) and those who had undergone separate replacement of the sinuses of Valsalva (remodeling [Yacoub] procedure) to find that the near-normal opening and closing characteristics of valves were achieved by the preservation of the shape and independent mobility of the sinuses. On the other hand, de Oliveira and associates,4
based on a decade of following up surgically treated patients, reported a reduced risk of postoperative aortic insufficiency in patients with reimplantation. They concluded that the implantation procedure secures the entire aortic valve inside the Dacron graft to prevent dilatation of the aortic root and concomitant regurgitation. However, they also recognized the elimination of sinuses as a shortcoming of their technique, and they, as well as other surgeons, have modified their reimplantation techniques to create graft pseudosinuses.5-7
Although it is generally assumed that smooth opening and closing relieved the abnormal stress and strain on leaflets, neither clinical observations3,8-10
nor experimental studies using a mock circulation11,12
could provide us with detailed information on the stress and/or strain distribution in the leaflet to clarify this assumption. An alternative approach to overcome such technical difficulty and obtain stress and strain distribution is a simulation study using finite element analysis. Grande-Allen and associates13
created finite element models of the aortic root and valve with a cylindrical graft, a tailored graft sutured just above the valve, and a pseudosinus graft, to find that the cylindrical graft created the greatest stress. Beck, Thubrikar, and Robicsek14
also compared the finite element models of the tubular aortic root and the root with sinuses under pressure to find that a stress concentration along the leaflet attachment develops only in the tubular root. Although these studies gave strong support to surgeons' beliefs, the simulation was done only by applying pressure to the model, thus totally ignoring the influence of the blood flow, which most of the surgeons and researchers are interested in. The capabilities of the fluid–structure interaction finite element method for the analysis of the aortic valve were first shown by Nicosia and colleagues15
using a commercial software package; subsequently, de Hart and coworkers16
analyzed the normal aortic root using the fictitious domain method. However, owing to the numerical instability of the method, an unphysiologically low Reynolds number flow was introduced.
We17
have developed a numerical approach for fluid–structure interaction analysis based on the arbitrary Lagrangian–Eulerian (ALE) finite element method and applied it to the multiphysics simulation of the heart.18-21
In this study, we applied this method to the analysis of the dynamics of aortic valve opening and closure in two different models of the aortic root, with or without the sinuses of Valsalva. It will be clearly shown that the sinuses, by promoting vortex formation, facilitate the smooth closure of the aortic valve and circumvent the development of abnormal stress in the leaflets.
| Method |
|---|
|
|
|---|
|
Computation
We used a strongly coupled fluid–structure interaction finite element analysis program, which we have developed17
and applied to various problems such as pulsation of the heart.18-21
In this program, the Navier–Stokes equation is described in ALE coordinates, which artificially deform and move according to the instantaneous deformation of the fluid–structure interface. Furthermore, the fluid meshes and structure meshes are generated to coincide with each other on the interface at the beginning of the analysis. Therefore, the geometric compatibility and balance of traction forces are automatically satisfied by merging both meshes on the interface. The disadvantage of the method is the distortion of ALE meshes. Because the fluid–structure interface moves largely in the case of a heart valve problem, excessive ALE mesh distortion occurs, resulting in the degradation of element performance. To avoid such a problem, we applied an automatic mesh reconnecting algorithm, which our group has developed. All of the program codes were written in the laboratory.
| Results |
|---|
|
|
|---|
Motion of the Leaflets
To quantify the motion of the leaflets and compare the simulation results with clinical observations, we traced the edge of the leaflet and plotted the distance from the axis of the aorta as a function of time, as we do in M-mode echocardiography (Figure 2, A
), and calculated the slow closing distance (SCD = [D1 – D2]/D1) and closing time (Figure 2, B). SCD values were greater for pseudosinus models (P) than for cylindrical models (C) under all of the conditions studied, but the difference was pronounced in response to low ejecting pressure (P vs C [%]: 8.9 vs, 5.1 [P110], 7.7 vs 4.6 [P100], and 6.3 vs 0.6 [P90]) (Figure 3, A
). On the contrary, although the differences were small, closing time tended to be longer with the cylindrical model (P vs C [s]: 0.22 vs 0.24 [P110], 0.24 vs 0.25 [P100], and 0.26 vs 0.28 [P90]) (Figure 3, B). We also calculated the rapid valve closing velocity (RVCV) as the ratio between D2 and rapid closing time (RCT in Figure 2, B). RVCV was also faster for cylindrical models (P vs C [cm/s]: 50.6 vs 52.1 [P110], 49.8 vs 58.5 [P100], and 30.8 vs 42.7 [P90]) (Figure 3, C). These results can be taken to indicate that, in the presence of pseudosinuses, the valvular leaflets initiate their motion from the earlier phase of ejection and gradually return to their closing positions. On the other hand, with the cylindrical model, valves are wide open until the late phase of ejection, during which the flow direction is reversed to increase the regurgitant fraction (P vs C [%]: 0.7 vs 4.4 [P110], 1.6 vs 4.2 [P100], and 0.4 vs 5.4 [P90]) (Figure 3, D).
|
|
) parallel (||) and perpendicular (tt) to the fiber orientation. As shown in Figure 4, A
, the leaflets in the cylindrical model experience greater peak strains, especially in the direction perpendicular to the fiber orientation (C vs P: 0.022 vs 0.012 [membrane
||], 0.042 vs 0.024 [membrane
tt], 0.069 vs 0.069 [bending
||], 0.087 vs 0.064 [bending
tt], 0.065 vs 0.060 [total
||], 0.089 vs 0.068 [total
tt]). Interestingly, the bending strains in the fiber direction did not differ appreciably between the two models. Comparison of peak stress (
) values followed similar pattern to the strain values (C vs P [kPa]: 15.9 vs 8.6 [membrane
||], 9.8 vs 5.6 [membrane
tt], 48.0 vs 48.1 [bending
||], 20.3 vs 14.8 [bending
tt], 45.3 vs 41.8 [total
||], 20.8 vs 15.7 [total
tt]) (Figure 4, B).
|
| Discussion |
|---|
|
|
|---|
Simulations of Aortic Root
Compared with studies using mock circulation,1,11,12
simulation studies using the finite element method have advantages in that (1) the shape and the material properties of the model and experimental conditions can be altered and controlled over a wide range, (2) detailed data on the distribution of flow velocity and pressure are available, and (3) stress/strain distribution in the aortic wall and/or valvular leaflet can be calculated. In particular, stress data are important but hard to obtain in clinical settings.
However, so far, we can find only a few simulation studies in which fluid–structure interactions were analyzed, probably owing to the computational difficulties.15,16,24
Among these studies, Nicosia and colleagues15
constructed an anatomically accurate 3-dimensional finite element model in which both the aortic root and valves were represented by Hughes–Liu shell elements. They analyzed the blood flow and the motion of the valve leaflet during ejection by using LS-Dyna—an explicit finite element commercial code. Their pioneering work showed the potential capabilities of the fluid–structure finite element analysis for heart valve problems. However, some unphysiologic conditions, such as the 98.5% reduction in the bending stiffness of the valve leaflets, were introduced to reproduce their pliability, which in turn enforced the scaling down of the magnitude of the peak diastolic pressure to avoid valve element distortion. After Nicosia and colleagues,15
de Hart and associates,16
using the fictitious-domain method, successfully showed that vortex formation in the sinuses of Valsalva is essential for the smooth operation of the aortic valve. In contrast to the present study, however, their computation was carried out with an unphysiologically low Reynolds number (
900) owing to the numerical instabilities inherent in the algorithm.
Furthermore, a comparison of two clinically relevant models, namely, cylindrical and pseudosinus, based on fluid–structure interaction analysis, has been made for the first time. Although we can find the finite element analysis studies comparing the principal tensile stress of leaflets between the cylindrical graft and pseudosinus graft,13,14
in those study, the calculation was carried out for only the diastolic period by applying the pressure, thus completely ignoring the influence of blood flow.
Comparison With Experimental and Clinical Studies
Using the time-resolved 3-dimensional magnetic resonance velocity mapping, Markl and associates10
compared vortex formation among patients who had undergone cylindrical graft (David I procedure) and those who had undergone neosinus graft (David V and David V-Smod [Stanford modification]), with normal volunteers as controls. Even though vorticity was increased in patients who underwent the David V procedure, because the difference did not reach statistical significance, they concluded that, although vortex formation was enhanced by the David V procedure, normal vorticity was preserved even without the sinus creation in the graft. Direct comparison with the current simulation results is difficult, because in this study a rigid tube with perfect cylindrical form was used, but a small degree of vortex formation was also identified in our simulation (Movie E2).
In a study comparing the motion of the leaflets between the patients who had undergone tube graft and those who had undergone remodeling (Yacoub) procedure, Leyh and coworkers3
reported that, in patients who had undergone remodeling, SCD was greater but closing time was shorter, consistent with the present results. A similar tendency was reported for the comparison between patients with a tube graft and those with a newly developed graft with a pseudosinus.8
However, the RVCV differed between the two studies. Leyh and coworkers3
found that the RVCV was slower for tube grafts, but De Paulis and associates8
reported results to the contrary.
In vitro studies may help to resolve this discrepancy. Fries and coworkers11
placed a porcine aortic root operated on with either the David (cylinder) or Yacoub (remodeling) procedure in a mock circulation and recorded the motion of the valve. The difference in SCD was similar to the clinical observation mentioned above, as well as the present results, but they found that the differences in RVCV were dependent on cardiac output and that the RVCV was greater for the David procedure only at low cardiac output, similar to the present result (Figure 3, C). These contradictory results may have arisen because the patients in the study by Leyh and colleagues3
had relatively high cardiac function (ejection fraction > 60%) compared with patients in other studies (about 50%).
As we mentioned earlier, the stresses and/or strains in the leaflet during ejection are hard to measure, but there has been a study in which strains were compared among the various surgical procedures.12
From the images of porcine aortic roots operated on differently and placed in the mock circulation, the authors of this study calculated the cusp-bending deformation index in mid-systole as the fold depth normalized to the size of the leaflet. Although the index of strain used in this study was just an estimate, their results demonstrated a greater strain associated with the reimplantation procedure (tube graft), compatible with our simulation result.
Implications
Although the creation of sinuses in remodeling procedures introduces the smooth closure of the aortic valve,3
a long-term follow-up study demonstrated a better prognosis for patients undergoing the reimplantation (David I) procedure in terms of freedom from aortic regurgitation,4
mainly because of the more reliable annular stabilization with this technique. However, the importance of the sinuses of Valsalva is widely recognized by cardiac surgeons, and various modifications to the reimplantation procedure have been proposed to achieve both annular stability and the creation of sinuses.6-8,25,26
Although these techniques are expected to enhance the long-term durability of the leaflet, greater knowledge of the functional anatomy of the aortic root is required to achieve optimal results.27
The current simulation would serve as a useful tool for designing tailor-made aortic grafts and, in fact, has provided us with new insight into the mechanics of the aortic valve leaflet during ejection.
Shown in Figure 5, A , are the stress distributions in the leaflets during the late phase of ejection for pseudosinus (upper panel) and cylindrical (lower panel) grafts (Movie E3). It can clearly be seen that the leaflet in the cylindrical graft is bent in its middle portion by being pushed by the retrograde flow from both sides (Movie E4). On the other hand, the leaflet in the pseudosinus graft has already returned halfway to the closed position, as evidenced by the large SCD, and accommodates the retrograde flow only on one side to avoid the abnormal bending stresses perpendicular to the fiber orientation (Figure 5, B). So far, the stresses in the leaflet have been discussed mainly during diastole when leaflets are in the closed position, but the current simulation study demonstrated another important time point when abnormal stresses could build up, as well as the importance of fluid–structure interaction in consideration of designing the aortic root grafts.
|
| Conclusion |
|---|
|
|
|---|
| Supplementary data |
|---|
|
|
|---|
| Table E1 |
|---|
|
|
|---|
|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. P. Urbanski, X. Zhan, H. Hijazi, M. Zacher, and A. Diegeler Valve-sparing aortic root repair without down-sizing of the annulus J. Thorac. Cardiovasc. Surg., February 1, 2012; 143(2): 294 - 302.e1. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. M. J. Bernhardt, H. Treede, M. Rybczynski, S. Sheikzadeh, J. F. Kersten, T. Meinertz, Y. von Kodolitsch, and H. Reichenspurner Comparison of aortic root replacement in patients with Marfan syndrome Eur J Cardiothorac Surg, November 1, 2011; 40(5): 1052 - 1057. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Kargar and M. H. Aazami Have we a surgical 'easy-hand' for replacement of the ascending aorta? Eur J Cardiothorac Surg, October 1, 2011; 40(4): 1033 - 1034. [Full Text] [PDF] |
||||
![]() |
A. S. Stewart, H. Takayama, and C. R. Smith Modified Bentall Operation With a Novel Biologic Valved Conduit Ann. Thorac. Surg., March 1, 2010; 89(3): 938 - 941. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Pillai, C. Ratnatunga, J. L. Soon, H. Kattach, A. Khalil, and X. Y. Jin 3f Prosthesis Aortic Cusp Replacement: Implantation Technique and Early Results Asian Cardiovasc Thorac Ann, February 1, 2010; 18(1): 13 - 16. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. E. David, S. Armstrong, M. Maganti, J. Colman, and T. J. Bradley Long-term results of aortic valve-sparing operations in patients with Marfan syndrome. J. Thorac. Cardiovasc. Surg., October 1, 2009; 138(4): 859 - 864. [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 |