Next Article in Journal
The Role of Hydrogel in Cardiac Repair and Regeneration for Myocardial Infarction: Recent Advances and Future Perspectives
Next Article in Special Issue
Influence of Inlet Boundary Conditions on the Prediction of Flow Field and Hemolysis in Blood Pumps Using Large-Eddy Simulation
Previous Article in Journal
Towards Ready-to-Use Iron-Crosslinked Alginate Beads as Mesenchymal Stem Cell Carriers
Previous Article in Special Issue
Computational Modelling of Cerebral Blood Flow Rate at Different Stages of Moyamoya Disease in Adults and Children
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Postoperative Remodeling Characteristics after Modular Inner Branched Stent-Graft Treatment of Aortic Arch Pathologies Using Computational Fluid Dynamics

1
College of Mechanical and Vehicle Engineering, Taiyuan University of Technology, Taiyuan 030024, China
2
Institute of Applied Mechanics, Taiyuan University of Technology, Taiyuan 030024, China
3
First Medical Center of Chinese PLA General Hospital, Department of Vascular Surgery, Beijing 100853, China
4
College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan 030024, China
*
Authors to whom correspondence should be addressed.
Bioengineering 2023, 10(2), 164; https://doi.org/10.3390/bioengineering10020164
Submission received: 29 October 2022 / Revised: 14 January 2023 / Accepted: 24 January 2023 / Published: 27 January 2023

Abstract

:
The modular inner branched stent-graft (MIBSG), a novel interventional therapy, has demonstrated good effects in the endovascular treatment of aortic arch pathologies, especially those involving the supra-aortic branches. Nevertheless, the long-term efficacy of the MIBSG and in-depth quantitative evaluation of postoperative outcomes remain to be examined. Moreover, the regularity of postoperative vascular remodeling induced by MIBSG implantation has yet to be explored. To address these questions, we constructed four models (normal, preoperative, 1 week postoperative, and 6 months postoperative) based on a single patient case to perform computational fluid dynamics simulations. The morphological and hemodynamic characteristics, including the velocity profile, flow rate distribution, and hemodynamic parameter distribution (wall shear stress and its derivative parameters), were investigated. After MIBSG implantation, the morphology of the supra-aortic branches changed significantly, and the branch point moved forward to the proximal ascending aorta. Moreover, the curvature radius of the aortic arch axis continued to change. These changes in morphology altered the characteristics of the flow field and wall shear stress distribution. As a result, the local forces exerted on the vessel wall by the blood led to vessel remodeling. This study provides insight into the vascular remodeling process after MIBSG implantation, which occurs as a result of the interplay between vascular morphological characteristics and blood flow characteristics.

Graphical Abstract

1. Introduction

Aortic arch pathologies (AAPs), such as aortic arch dissection and aortic arch aneurysm, are sometimes associated with severe outcomes. Endovascular aortic repair (EVAR) is widely used as a minimally invasive treatment for AAPs in clinical practice [1,2,3,4]. With advances in instrument technology and the accumulation of clinical experience, endovascular technology for AAPs has made great progress. The branched stent, fenestration stent, and other technologies have been developed [2,5,6,7,8]. However, challenges remain due to the curved structure of the aortic arch and the need to preserve the important brachiocephalic branches of the arch. Branched graft techniques require custom-made devices and complex procedures, which makes them less promising in cases requiring urgent EVAR [9,10,11,12]. For example, the gap between the branch stent and the main stent increases the incidence of postoperative leakage complications [13,14]. Moreover, fenestrated stent technology requires the adjustment of conventional aortic grafts based on the anatomical morphology of the aortic arch. During release, the preopening window of the graft is aligned with the opening of the branch vessels, which increases the procedural difficulty [15,16].
The modular inner branched stent-graft (MIBSG) is a novel stent with an innovative design and a simple release method [17,18,19] that can be used to reconstruct the aortic arch and its upper branches. The MIBSG system comprises three modules, as shown in Figure 1. The first module is a straight stent-graft with double inner branches anchoring to the ascending aorta. The second module comprises covered stents that can be used to reconstruct the superior arch branch. The third module is a flexible thoracic stent-graft that is intended to isolate the aneurysm. The stent-graft is designed using modularization techniques, which makes it more adaptable (the stent is not patient-specific; it is an off-the-shelf application) with an uncomplicated implantation procedure. The technical feasibility of the MIBSG has been verified in animal experiments and clinical trials. Dozens of patients have been treated with this approach, and the early postoperative results are satisfactory. However, to date, the clinical evaluation of MIBSG efficacy has been based on imaging evaluation, which is insufficient for long-term assessment. Nevertheless, computational fluid dynamics (CFD) provide a promising alternative.
While physiological indicators, such as blood perfusion are the concern of physicians, hemodynamic researchers care about the flow field characteristics related to specific pathologies, such as aortic arch aneurysm and type B aortic dissections, as well as the normal flow field characteristics.
CFD is extensively used to assess hemodynamic parameters, such as flow velocity, blood pressure, and wall shear stress (WSS), to analyze the long-term efficacy of EVAR. With CFD, blood vessels before and immediately after surgery are reconstructed using clinical image data. Then, the reconstructed models are used for CFD or fluid—structure interaction computations. Using this approach, the efficacy of surgical treatment or other assessments can be carried out based on the analysis of the numerical results.
In addition, sometimes CFD is the only choice if we need to elucidate the influent of a specific factor on the development of the course of disease. For example, as well known, the functions of stent-graft are to reinforce the wall of an aneurysm to prevent it from rupture, or to seal the tear in the aorta in the case of a dissection. In some cases, the stent-graft does not fit perfectly with the aortic wall, and a so-called bird-beak configuration forms. The bird-beak effect may lead to reduction in the lumen patency, and its long-term influence on the pressure drop and WSS-based parameters is revealed by the usage of CFD-based methods [20,21].
Moreover, CFD research on other important factors, such as the aortic diameter [22] and the beginning position for an octopus endograft [23], are helpful to the clinical predictions.
Similarly, more complicated studies on the surgical intervention techniques are performed to analyze their advantages and disadvantages to a specific patient before surgery [24,25].
While the blood flow field has important physiological meanings, it is attractive to extract an “ideal” structure from real tissues to study the effects of the predominant structure characteristics only [26,27]. This approach is significant for the understanding of typical flow patterns in tissues, such as secondary flow and helical flow.
It is also important to note that vascular remodeling after surgery is a long-term process. The geometric characteristics of blood vessels change continuously over time, which also causes changes in hemodynamic characteristics [28]. Therefore, the characteristics of intraluminal blood flow immediately after interventional therapy are not the sole predictor of long-term outcomes.
In this study, we evaluated the hemodynamic characteristics of a single case of MIBSG implantation. A preoperative model and two postoperative models (immediately and 6 months after MIBSG implantation) were reconstructed from computed tomography angiography (CTA) images. An optimized normal geometry model was used as a control. Using CFD, flow quantities were extracted, and their distributions were examined by referring to geometric changes to explore the interactions between vascular morphology and hemodynamics.

2. Materials and Methods

2.1. Image-Based Model Reconstruction

The patient, a 74-year-old male, suffered from hypertension and presented with a cystic aortic arch aneurysm. The aneurysm was excluded, and two branches of the innominate artery (IA) and the left common carotid artery (LCCA) were reconstructed by MIBSG implantation. The left subclavian artery (LSA) was occluded after surgery, and LCCA-LSA bypass could take effect when needed (Figure 2).
Based on CTA images, the three-dimensional geometry was reconstructed using Mimics commercial software (version 19.0; Materialise, Plymouth, MI, USA). The geometry of the MIBSG was not reconstructed separately, but the composite structure of MIBSG-aortic arch was reconstructed as a whole. This is due to the fact that the simulation used in this research is pure CFD. On the other hand, the reconstruction of stent-graft from postoperative CTA datasets is technically difficult. The preoperative and postoperative geometry (1 week and 6 months after MIBSG implantation) were reconstructed with a supra-arch branch wall thickness of 1.0 mm and an aortic wall thickness of 2.0 mm. To illustrate the utility of MIBSG implantation in restoring the normal flow pattern, a normal geometry model was constructed with the aneurysm removed (Figure 2). All of the models started above the sinus of the artery and ended at the bifurcation point of the renal artery.

2.2. Computational Framework

The blood flow through the aortic arch was modeled as laminar flow and simulated using the laminar flow module of commercial finite element analysis software (COMSOL 5.5). Generally, turbulent flow patterns are always observed in the blood flow inside the aneurysm and the proximal or distal portion of the aortic arch due to the complex geometry structure and relatively high velocity, i.e., large Reynolds number. Therefore, the use of turbulent models, such as Reynolds-averaged Navier-Stokes (RANS) or large eddy simulation (LES) and even direct numerical simulation (DNS) is natural in hemodynamic researches [29,30]. However, physiologically, the turbulent flow in aorta is harmful, and the high degree of pulsatility of blood flow depresses the development of turbulence [31]. Moreover, research by Xu et al. [23] showed that the laminar flow model with adequately fine mesh resolutions near the walls is capable of capturing the main flow patterns. Therefore, in this research, the laminar flow model is employed with six boundary layers to capture the near-wall blood flow behaviors. Mesh-independent solutions were validated using sufficiently fine elements, for example, 795,698 tetrahedral elements were used for the normal model. The time step used by the implicit solver is adaptive with the maximum value of 0.05 s and the relative error is set to 0.005 between two inner iterations.
The blood is naturally a non-Newtonian fluid; however, in aortas, its viscosity changes slightly and thus a Newtonian constitutive relationship is used for the blood. The incompressible Navier-Stokes equation was used to describe its flow behavior. The density and viscosity were set at 1060 kg/m3 and 0.0035 Pa·s, respectively.

2.3. Boundary Conditions

In the simulation, the outer wall of the aortic arch was fixed, which induced displacement of the order μm of the inner wall, which had a trivial effect on the characteristics of the flow field. Therefore, the inner walls were maintained as rigid in all four simulations, and the two-way fluid—structure interaction between the blood and vessels was dismissed. Blood flow in the aorta is pulsatile, thus the boundary conditions of the inlet and outlet must be set accordingly. Considering the convergence of CFD problems and physiological properties, velocity condition was used as the inlet and pressure condition was used as the outlet in all four models. The inlet and outlet were extruded outward to five-times their diameter to develop the flow pattern. The velocity condition profile was based on the Doppler ultrasound measurement, and the pressure condition was derived from the clinical examinations. The expressions for velocity condition and pressure condition were as follows:
v t = 0.3825 cos 6.6667 π t + 0.6775 ,                                 0 < t 0.3   s 0.1405 cos 3.3333 π t 0.3 + 0.4355 ,         0.3 < t 0.6   s 0.1405 cos 2.5 π t 1 + 0.4355 ,                           0.6 < t 1   s
p t = 25 cos 4 π t + 115 ,     0 < t 0.35   s p 0.35 90 cos 0.7692 π t 0.35 + π 1.5 + p 0.35 ,     0.35 < t 1   s
where v and p are velocity and pressure, respectively. Their profiles are shown in Figure 3.

2.4. Related Indicators

Hemodynamic parameters, such as time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT), were derived from the WSS. The distributions and trends of these parameters provide clues to the evolution of the restored flow patterns, which then cause the physiological response of cells to the shear stress and strain exerted by blood. Therefore, hemodynamic parameters carry important physiological meanings, and they serve as a bridge between CFD and clinical practice. In this study, hemodynamic parameters were calculated in the post-process to assess the treatment efficacy.
TAWSS is defined as a scalar that indicates the time-averaged magnitude of the surface shear force, as follows:
TAWSS = τ abs = 1 T 0 T τ w d t
Taylor, Hughes, and Zarins [32] defined another scalar as the magnitude of the time-averaged mean shear force, as follows:
τ mean = 1 T 0 T τ w d t
Then, the OSI was calculated using the above two scalars, as follows:
OSI = 1 2 1 τ mean τ abs
Finally, the RRT was evaluated using the TAWSS and OSI, as follows:
RRT = 1 TAWSS 1 2   OSI

3. Results

3.1. Morphological Description

The morphology of the branches above the aortic arch changed after MIBSG implantation. As shown in Figure 2, the IA and LCCA branches were reconstructed, and the LSA was sealed after surgery (Figure 2). The ostia of supra-arch branches were located in the first module of the MIBSG, moved from the preoperative arch to the nearby proximal end.
To explore the geometric characteristics of the aortic arch at different stages, Mimics software was used to capture the vessel centerlines of the models (the branches were uncontained). The morphological characteristics of the aorta were analyzed at the three observational stages and in the normal model. The highest points of the height direction of the models were set to overlap each other for comparison (Figure 4). The centerline direction of the aortic arch changed greatly owing to the existence of the stent-graft. The radius of the aortic arch curvature significantly decreased 1 week after MIBSG implantation, but it increased after 6 months.

3.2. Blood Perfusion

The perfusion rate has always been a consideration for cardiovascular surgery clinicians. In this study, the inlet and outlet flow rates were calculated via surface integration. Figure 5 shows the average perfusion rate of each outlet within one pulsation cycle. Figure 5 shows that, compared with the perfusion rate of the normal model, the existence of an arch aneurysm reduced the total perfusion rate of the superior arch branch. There was also a corresponding decrease in the total volume of the two branches (right common carotid artery [RCCA] + LCCA), which supplies blood to the brain. Six months later, the total perfusion rate of the superior arch branch significantly decreased, especially for the right subclavian artery (RSA) branch. The perfusion rate of the RCCA significantly decreased, and the RSA was compensatory.

3.3. Hemodynamic Description

Flow pattern
The profiles of peak systolic streamlines visibly highlight the flow patterns of the models at different stages after MIBSG implantation (Figure 6). Compared with Case A (the normal model), in Case B (preoperative), the aneurysm had a larger impact on the blood flow pattern, there was a clear turbulent flow in the aneurysm, and the disturbed flow in the proximal lesser curvature of the descending aorta was strengthened (Figure 6). Additionally, a clear backflow is observed inside the aneurysm. After treatment, the branched stent overlapped with the aortic stent in the first module and was compressed. The influence of this geometric feature on blood flow was as follows. First, the location and lengthening of the branches strengthened the disturbed blood flow within the branches. Second, clear blood flow disturbance could be seen at the proximal end of the third module. Third, there were more vortices of varying intensity in the lesser curvature of the descending aorta at both postoperative stages for the increased curvature, and thus strengthened the secondary flow here (1 week [Case C] and 6 months [Case D]; Figure 6).
Pressure
No significant pressure distribution differences were observed between the model with the presence of an arch aneurysm and the normal model; however, a uniform decrease in pressure existed in both models. One week after MIBSG implantation, the decrease in pressure of Case C increased, and an area with an adverse pressure gradient in the lesser curvature of the proximal descending aorta appeared. However, the overall decrease in pressure of the model 6 months after MIBSG implantation was smaller than observed 1 week after MIBSG implantation, and the adverse pressure gradient on the side of the lesser curvature of the proximal descending aorta was also relieved (Figure 7).
WSS-based hemodynamic indicators
WSS-based hemodynamic indicators, such as time-averaged TAWSS, OSI, and RRT, have important physiological meanings. According to Malek [33], regions with a TAWSS of <0.4 Pa, an OSI of >0.25, and an RRT of >5 Pa−1 show a stronger tendency to develop arterial atherosclerosis. In an attempt to clearly distinguish the tendency of different postoperative complications, high-risk regions with abnormal values for these factors were classified.
Figure 8 shows the distribution of areas with a TAWSS of <0.4 Pa in the four models. Compared with the model with a normal structure, the area with a TAWSS of <0.4 Pa increased in the presence of an aneurysm, and the affected area was located mainly in the aneurysm and lesser curvature of the descending aorta. One week after MIBSG implantation, this area decreased in size and was even smaller than in the normal model. However, 6 months after MIBSG implantation, this area increased in size and was primarily distributed in the third module stent and descending aorta. One week and 6 months after MIBSG implantation, the area with a TAWSS of <0.4 Pa appeared on every branch, which was not apparent in the preoperative and normal models.
The higher the OSI, the higher the intensity of the shear stress oscillation due to fluctuations of flow. The regional distribution of an OSI of >0.25 and the area statistics in the four models are shown in Figure 9. The area with an OSI of >0.25 was smaller in the normal model, and it was predominantly distributed in regions with disordered blood flow, such as areas proximal to the great curvature of the ascending aorta and vascular branch bifurcation (Case A in Figure 9). The presence of an aneurysm decreased the overall area with an OSI of >0.25 but increased the area with an OSI of >0.25 on the branching vessels. The area with an OSI of >0.25 increased significantly after stent implantation. The distribution characteristics of these regions were analogous in the 1-week and 6-month postoperative models. The region with a high OSI increased clearly in the proximal part of the third module, the proximal part of the lesser curvature of the descending aorta, and the two branch vessels of the two postoperative models (Cases C and D of Figure 9).
A higher RRT (particle retention time) indicates that the physical material in the blood vessel stays there longer and may eventually form deposits. As shown in Figure 10, the distribution characteristics of regions with an RRT of >5 Pa−1 are similar to those with a TAWSS of <0.4 Pa. This suggests a long residence time of substances in the blood in these regions.

4. Discussion

As a new stent technology, MIBSG with novel design was proposed to address two main problems: Poor tolerance of brain tissue to ischemia and the reconstruction difficulty caused by a large arch curvature. Despite good clinical imaging results after MIBSG, there is still no detailed analysis of local hemodynamics to predict the outcomes of endovascular surgery. To fill this gap in knowledge, we examined the postoperative blood flow field using image-based CFD. The blood flow characteristics pre- and post-EVAR have been previously analyzed to evaluate the pathological development and postoperative efficacy [21,34,35,36,37,38,39]. Here, we combined morphological characteristics and computational hemodynamics to provide a detailed analysis of the impact of MIBSG implantation in the follow-up period of a single patient.
For numerical simulation, the geometry profile of the solution domain provides the boundary of a fluid dynamic problem, and thus imposes a remarkable influence on the flow characteristics, e.g., the flow pattern [40,41,42]. This drives the need of an accompanying analysis of the morphology evolution when the analysis of flow patterns is performed.
The aortic arch is complex for its curvature in three dimensions, and the aneurysm enlarges the local diameter of the aortic arch. According to fluid dynamic theory, the bend of the arch generates a secondary flow [43,44]. The secondary flow occurs due to the velocity difference between the lesser and greater curvature of the arch and thus the difference of centrifugal forces when the blood flows through the aortic arch. This indicates that the strength of the secondary flow is related directly with the curvature of the aortic arch. This is why the vortex region was more evident on the lesser curvature side of the descending aorta in the Case C model (Figure 6), compared to the preoperative model. This is one of the factors that gives rise to the risky areas, which consists of great physiological meanings.
In addition to the bending, the aortic arch has significant torsion, which leads to the helpful helical flow.
Another notable morphology feature is the ostia of supra-arch branches, which also consists of the bifurcation sites of the blood flow. As blood enters the bifurcation sites, the flux distributes among the branches. The ostia of supra-arch branches moved from the preoperative arch to the nearby proximal end after surgery (Figure 2), i.e., the bifurcation sites changed greatly. Another change in the ostia is that the LSA is blocked. These factors cause the blood perfusion to change significantly, as shown in Section 3.2. When the blood enters the branch, flow separation may occur, resulting in low wall shear stress regions [31].
With respect to the aneurysm, for its growth deviation from the main flow direction, an adverse pressure gradient arises near the neck of the aneurysm, which gives birth to a vortex [22,44,45], and results in the slow speed of blood flow. Additionally, the longer resident time of blood particles in this region (Case B) provides another formative factor for the risky areas.
As discussed above, after the treatment of AAPs, the morphologies of the aortic arch change significantly, which alters the existing flows. This provides a clue for the assessment of different treatment approaches before operation [13,20,22,23,26,28].
Another aspect to be noticed in the treatment of AAPs with stent-graft is that a difference in the rigidity of the stents and adjacent vessels creates a mismatch between their displacement [46,47,48,49,50]. Moreover, stent-graft placement inside the aortic arch modifies the aortic curvature [21,51], and the mechanical interaction between the vascular wall and the stent-graft may cause changes in anatomical features [52]. The change in vascular morphology and the consequent hemodynamics, which are greatly affected by the anatomy, greatly influence the long-term outcomes of stent placement [21,53]. In this study, we found that the central axis of the aortic arch changed after stent implantation. The radius of the aortic arch curvature was significantly smaller 1 week after MIBSG implantation than in the normal model and the preoperative model. This geometrical characteristic increases the blood flow resistance, which is reflected by the larger pressure decrease in the whole model of Case C (Figure 7). The change in blood flow characteristics; the mismatch in mechanical properties, such as Young’s modulus and Poisson’s ratio; and displacement between the stent and blood vessel can cause continuous blood vessel remodeling. As a result, the radius of the curvature of the central axis of the blood vessel increased 6 months after surgery (Case D) compared with 1 week after surgery (Case C), which indicates that the stress within the vessels and MIBSG is released to some extent. Additionally, clear changes in the blood perfusion rate of each supra-aortic branch were observed as shown in Figure 5.

5. Limitations

This study has some limitations that should be noted. First, only a single case of MIBSG implantation was analyzed. Therefore, more cases should be assessed to explore the regularity of vascular remodeling after MIBSG implantation. Second, only the blood flow field is investigated using the CFD technique. The deformation of MIBSG and the aortic arch, the interaction between MIBSG and aortic arch, i.e., the mechanical mismatch between them, the interaction between MIBSG and blood flow, and between the aortic arch and blood flow are left out of consideration. This neglection is to some extent reasonable if our focus is on the main flow characteristics. On the other hand, if we need to assess the long-term durability of MIBSG, its interactions with the aortic arch and the blood flow are necessary to be taken into account. Third, the same boundary conditions were used in all four models, which is ideal and cannot capture the real clinical scenario. This is attributed to the difficulty of measuring the velocity profile and pressure distribution of arbitrary vessel cross-sections in vivo. Without doubt, velocity and pressure conditions at the cross-section of the inlets and outlets are crucial for the numerical results of CFD. Therefore, more non-invasive measurement methodologies and instruments with high resolutions are needed, such as 4D-flow MRI. It is possible to retrieve the velocity of any point within the scan, whereas the resolution near the vessel walls is still not satisfying. Compared with other measurements, the flow field visualized by 4D-flow MRI provides more reliable boundary conditions and can be used to calibrate the results of CFD. We have tried to apply this methodology in this research in the early stage, but image artifacts restricted its reliability. This is to be addressed in further research.

6. Conclusions

In this study, CFD was used to assess the long-term efficacy of MIBSG implantation. Compared with the preoperative model, the flow pattern of the postoperative models changed greatly, with more vortices in the lesser curvature of the descending aorta. Compared with the trend in hemodynamic parameters in the model immediately after surgery, the hemodynamic parameters of subsequent models showed an apparent fallback, which indicates that remodeling of the aortic arch morphology was present after surgery.

Author Contributions

F.L. contributed to the design of the study and reconstructed the CADs. Y.Z. contributed to the reconstruction and analysis of CFD results. H.S. performed the CFD simulation. H.Z. validated the numerical results. L.C. supervised the research. W.G. co-supervised the research. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by National Natural Science Foundation of China [12002231], the Beijing Municipal Natural Science Foundation [7212095], “1331 project” Key Innovation Teams of Shanxi Province, and Central Guidance on Local Science and Technology Development Found of Shanxi Province (YDZJSX2021A021).

Institutional Review Board Statement

The study was approved by the Ethics Committee of the Chinese People’s Liberation Army General Hospital (S2018-230-172 01) and adhered to the principles of the Declaration of Helsinki.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The original data were available from the corresponding author upon an appropriate request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CFDComputational fluid dynamics
AAPsAortic arch pathologies
CTAComputed tomography angiography
IAInnominate artery
LCCALeft common carotid artery
RCCARight common carotid artery
LSALeft subclavian artery
RSARight subclavian artery
MIBSGModular inner branched stent-graft
OSIOscillatory shear index
RRTRelative residence time
TAWSSTime-averaged wall shear stress
EVAREndovascular aortic repair
WSSWall shear stress

References

  1. Kamman, A.V.; de Beaufort, H.W.L.; van Bogerijen, G.H.W.; Nauta, F.J.H.; Heijmen, R.H.; Moll, F.L.; van Herwaarden, J.A.; Trimarchi, S. Contemporary Management Strategies for Chronic Type B Aortic Dissections: A Systematic Review. PLoS ONE 2016, 11, e0154930. [Google Scholar] [CrossRef] [PubMed]
  2. Lioupis, C.; Abraham, C.Z. Results and Challenges for the Endovascular Repair of Aortic Arch Aneurysms. Perspect. Vasc. Surg. Endovasc. Ther. 2011, 23, 202–213. [Google Scholar] [CrossRef] [PubMed]
  3. Schenning, R.C.; Al-Hakim, R. Aortic Dissection: Branched, Fenestrated, and Parallel Aortic Stent Grafts in the Ascending Aorta and Arch. Tech. Vasc. Interv. Radiol. 2021, 24, 100754. [Google Scholar] [CrossRef]
  4. Manetta, F.; Newman, J.; Mattia, A. Indications for Thoracic EndoVascular Aortic Repair (TEVAR): A Brief Review. Int. J. Angiol. 2018, 27, 177–184. [Google Scholar] [CrossRef]
  5. Maurel, B.; Mastracci, T.M.; Spear, R.; Hertault, A.; Azzaoui, R.; Sobocinski, J.; Haulon, S. Branched and Fenestrated Options to Treat Aortic Arch Aneurysms. J. Cardiovasc. Surg. 2016, 57, 686–697. [Google Scholar]
  6. Blanco Amil, C.L.; Mestres Alomar, G.; Guarnaccia, G.; Luoni, G.; Yugueros Castellnou, X.; Chiara Vigliotti, R.; Ramses, R.; Riambau, V. The Initial Experience on Branched and Fenestrated Endografts in the Aortic Arch. A Systematic Review. Ann. Vasc. Surg. 2021, 75, 29–44. [Google Scholar] [CrossRef] [PubMed]
  7. Yuri, K.; Kimura, N.; Hori, D.; Yamaguchi, A.; Adachi, H. A Challenging Treatment for Aortic Arch Aneurysm with Fenestrated Stent Graft. Ann. Thorac. Surg. 2017, 104, 1915–1922. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Hughes, G.C.; Vekstein, A. Current State of Hybrid Solutions for Aortic Arch Aneurysms. Ann. Cardiothorac. Surg. 2021, 10, 731–743. [Google Scholar] [CrossRef]
  9. Milne, C.P.E.; Haulon, S.; Oderich, G.S. Technical Aspects and Results of Branched Endografts for Repair of Aortic Arch Aneurysms. In Endovascular Aortic Repair: Current Techniques with Fenestrated, Branched and Parallel Stent-Grafts; Oderich, G.S., Ed.; Springer International Publishing: Cham, Switzerland, 2017; pp. 507–529. ISBN 978-3-319-15192-2. [Google Scholar]
  10. Inoue, K.; Hosokawa, H.; Iwase, T.; Sato, M.; Yoshida, Y.; Ueno, K.; Akiyoshi, T.; Terumitsu, T.; Shunichi, T.; Takahiko, S. Aortic Arch Reconstruction by Transluminally Placed Endovascular Branched Stent Graft. Circulation 1999, 100, II–316–II–321. [Google Scholar] [CrossRef]
  11. van Bakel, T.M.; de Beaufort, H.W.; Trimarchi, S.; Marrocco-Trischitta, M.M.; Bismuth, J.; Moll, F.L.; Patel, H.J.; van Herwaarden, J.A. Status of Branched Endovascular Aortic Arch Repair. Ann. Cardiothorac. Surg. 2018, 7, 406–413. [Google Scholar] [CrossRef] [Green Version]
  12. Hauck, S.R.; Kupferthaler, A.; Kern, M.; Rousseau, H.; Ferrer, C.; Iwakoshi, S.; Sakaguchi, S.; Stelzmuller, M.-E.; Ehrlich, M.; Loewe, C.; et al. Branched versus Fenestrated Thoracic Endovascular Aortic Repair in the Aortic Arch: A Multicenter Comparison. J. Thorac. Cardiovasc. Surg. 2022, 164, 1379–1389.e1. [Google Scholar] [CrossRef]
  13. Boufi, M.; Aouini, F.; Guivier-Curien, C.; Dona, B.; Loundou, A.D.; Deplano, V.; Alimi, Y.S. Examination of Factors in Type I Endoleak Development after Thoracic Endovascular Repair. J. Vasc. Surg. 2015, 61, 317–323. [Google Scholar] [CrossRef] [Green Version]
  14. Spear, R.; Haulon, S.; Ohki, T.; Tsilimparis, N.; Kanaoka, Y.; Milne, C.P.E.; Debus, S.; Takizawa, R.; Kölbel, T. Subsequent Results for Arch Aneurysm Repair with Inner Branched Endografts. J. Vasc. Surg. 2016, 63, 1129. [Google Scholar] [CrossRef]
  15. Kasprzak, P.M.; Mueller, T.; Loibnegger, A.; Kobuch, R.; Zorger, N.; Pfister, K. Fenestrated and Branched Stent Prostheses in the Aortic Arch: First Experiences. Gefasschirurgie 2009, 14, 198–205. [Google Scholar] [CrossRef]
  16. Simons, J.P.; Crawford, A.S.; Flanagan, C.P.; Aiello, F.A.; Arous, E.J.; Judelson, D.R.; Messina, L.M.; Robichaud, D.; Valliere, S.A.; Schanzer, A. Evolution of Fenestrated/Branched Endovascular Aortic Aneurysm Repair Complexity and Outcomes at an Organized Center for the Treatment of Complex Aortic Disease. J. Vasc. Surg. 2021, 73, 1148. [Google Scholar] [CrossRef]
  17. Wei, G.; Xin, J.; Yang, D.; Liu, X.; Tai, Y.; Zhang, H.; Liang, F.; Zhang, G. A New Modular Stent Graft to Reconstruct Aortic Arch. Eur. J. Vasc. Endovasc. Surg. 2009, 37, 560–565. [Google Scholar] [CrossRef] [Green Version]
  18. Guo, W.; Zhang, H.; Liu, X.; Ren, W.; Wang, Y. Endovascular Repair of Aortic Arch Aneurysm with a New Modular Double Inner Branch Stent Graft. Ann. Vasc. Surg. 2021, 77, 347.e1–347.e5. [Google Scholar] [CrossRef]
  19. Zhu, Y.; Li, F.; Zhang, H.; Song, H.; Ma, X.; Cao, L.; Zhang, W.; Guo, W. Hemodynamic Numerical Simulation of Aortic Arch Modular Inner Branched Stent-Graft in Eight Early Patients from the First-in-Human Case Series. Front. Cardiovasc. Med. 2022, 9, 981546. [Google Scholar] [CrossRef]
  20. van Bogerijen, G.H.W.; Auricchio, F.; Conti, M.; Lefieux, A.; Reali, A.; Veneziani, A.; Tolenaar, J.L.; Moll, F.L.; Rampoldi, V.; Trimarchi, S. Aortic Hemodynamics After Thoracic Endovascular Aortic Repair, With Particular Attention to the Bird-Beak Configuration. J. Endovasc. Ther. 2014, 21, 791–802. [Google Scholar] [CrossRef]
  21. Gallo, D.; Lefieux, A.; Morganti, S.; Veneziani, A.; Reali, A.; Auricchio, F.; Conti, M.; Morbiducci, U. A Patient-Specific Follow up Study of the Impact of Thoracic Endovascular Repair (TEVAR) on Aortic Anatomy and on Post-Operative Hemodynamics. Comput. Fluids 2016, 141, 54–61. [Google Scholar] [CrossRef]
  22. Etli, M.; Canbolat, G.; Karahan, O.; Koru, M. Numerical Investigation of Patient-Specific Thoracic Aortic Aneurysms and Comparison with Normal Subject via Computational Fluid Dynamics (CFD). Med. Biol. Eng. Comput. 2021, 59, 71–84. [Google Scholar] [CrossRef] [PubMed]
  23. Xu, H.; Mei, Y.; Han, X.; Wei, J.; Watton, P.N.; Jia, W.; Li, A.; Chen, D.; Xiong, J. Optimization Schemes for Endovascular Repair with Parallel Technique Based on Hemodynamic Analyses. Int. J. Numer. Methods Biomed. Eng. 2019, 35, e3197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Nardi, A.; Even-Chen, B.; Avrahami, I. Experimental and Numerical Study of the Flow Dynamics in Treatment Approaches for Aortic Arch Aneurysms. In Aortic Aneurysm; IntechOpen: London, UK, 2017. [Google Scholar] [CrossRef] [Green Version]
  25. Qiao, Y.; Mao, L.; Ding, Y.; Fan, J.; Zhu, T.; Luo, K. Hemodynamic Consequences of TEVAR with in Situ Double Fenestrations of Left Carotid Artery and Left Subclavian Artery. Med. Eng. Phys. 2020, 76, 32–39. [Google Scholar] [CrossRef] [PubMed]
  26. Gur, H.B.; Halak, M.; Brand, M. Blood Flow in the Abdominal Aorta Post’ Chimney’ Endovascular Aneurysm Repair. In Proceedings of the 9th EUROSIM Congress on Modelling and Simulation, EUROSIM 2016, The 57th SIMS Conference on Simulation and Modelling SIMS2016, Oulu, Finland, 13–15 September 2016; pp. 667–672. [Google Scholar] [CrossRef]
  27. Liu, M.; Sun, A.; Deng, X. Hemodynamic Performance within Crossed Stent Grafts: Computational and Experimental Study on the Effect of Cross Position and Angle. BioMed. Eng. OnLine 2018, 17, 85. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Duraiswamy, N.; Schoephoerster, R.T.; Moreno, M.R.; Moore, J.E. Stented Artery Flow Patterns and Their Effects on the Artery Wall. Annu. Rev. Fluid Mech. 2007, 39, 357–382. [Google Scholar] [CrossRef]
  29. Numata, S.; Itatani, K.; Kanda, K.; Doi, K.; Yamazaki, S.; Morimoto, K.; Manabe, K.; Ikemoto, K.; Yaku, H. Blood Flow Analysis of the Aortic Arch Using Computational Fluid Dynamics. Eur. J. Cardio-Thorac. Surg. 2016, 49, 1578–1585. [Google Scholar] [CrossRef] [Green Version]
  30. Mirfendereski, S.; Park, J.S. Direct Numerical Simulation of a Pulsatile Flow in a Stenotic Channel Using Immersed Boundary Method. Eng. Rep. 2022, 4, e12444. [Google Scholar] [CrossRef]
  31. Chandran, K.B.; Rittgers, S.E.; Yoganathan, A.P. Biofluid Mechanics: The Human Circulation, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2012; ISBN 978-1-4398-4516-5. [Google Scholar]
  32. Taylor, C.A.; Hughes, T.J.R.; Zarins, C.K. Effect of Exercise on Hemodynamic Conditions in the Abdominal Aorta. J. Vasc. Surg. 1999, 29, 1077–1089. [Google Scholar] [CrossRef] [Green Version]
  33. Malek, A.M. Hemodynamic Shear Stress and Its Role in Atherosclerosis. JAMA 1999, 282, 2035. [Google Scholar] [CrossRef]
  34. Qiao, Y.; Fan, J.; Ding, Y.; Zhu, T.; Luo, K. A Primary Computational Fluid Dynamics Study of Pre- and Post-TEVAR With Intentional Left Subclavian Artery Coverage in a Type B Aortic Dissection. J. Biomech. Eng.-Trans. Asme 2019, 141, 111002. [Google Scholar] [CrossRef]
  35. Tseng, F.S.; Soong, T.K.; Syn, N.; Ong, C.W.; Liang, L.H.; Choong, A. Computational Fluid Dynamics in Complex Aortic Surgery: Applications, Prospects and Challenges. J. Surg. Simul. 2017, 4, 1–4. [Google Scholar] [CrossRef]
  36. Tran, K.; Kaladji, A.; Yang, W.; Marsden, A.L.; Lee, J. Patient-Specific Computational Flow Simulation Reveals Differences in Aortic Hemodynamics Between Two versus Four Vessel Fenestrated Endovascular Aneurysm Repair. Circulation 2021, 144, A13545. [Google Scholar] [CrossRef]
  37. Auricchio, F.; Conti, M.; Lefieux, A.; Morganti, S.; Reali, A.; Sardanelli, F.; Secchi, F.; Trimarchi, S.; Veneziani, A. Patient-Specific Analysis of Post-Operative Aortic Hemodynamics: A Focus on Thoracic Endovascular Repair (TEVAR). Comput. Mech. 2014, 54, 943–953. [Google Scholar] [CrossRef]
  38. Midulla, M.; Moreno, R.; Negre-Salvayre, A.; Beregi, J.-P.; Haulon, S.; Loffroy, R.; Dake, M.; Rousseau, H. Impact of Thoracic Endografting on the Hemodynamics of the Native Aorta: Pre- and Postoperative Assessments of Wall Shear Stress and Vorticity Using Computational Fluid Dynamics. J. Endovasc. Ther. 2021, 28, 63–69. [Google Scholar] [CrossRef] [PubMed]
  39. Xu, H.; Li, Z.; Dong, H.; Zhang, Y.; Wei, J.; Watton, P.N.; Guo, W.; Chen, D.; Xiong, J. Hemodynamic Parameters That May Predict False-Lumen Growth in Type-B Aortic Dissection after Endovascular Repair: A Preliminary Study on Long-Term Multiple Follow-Ups. Med. Eng. Phys. 2017, 50, 12–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Dottori, J.; Casciaro, M.; Craiem, D.; El-Batti, S.; Mousseaux, E.; Alsac, J.-M.; Larrabide, I. Regional Assessment of Vascular Morphology and Hemodynamics: Methodology and Evaluation for Abdominal Aortic Aneurysms after Endovascular Repair. Comput. Methods Biomech. Biomed. Eng. 2020, 23, 1060–1070. [Google Scholar] [CrossRef]
  41. Stefanov, F.; McGloughlin, T.; Morris, L. A Computational Assessment of the Hemodynamic Effects of Crossed and Non-Crossed Bifurcated Stent-Graft Devices for the Treatment of Abdominal Aortic Aneurysms. Med. Eng. Phys. 2016, 38, 1458–1473. [Google Scholar] [CrossRef]
  42. Kokkalis, E.; Aristokleous, N.; Houston, J.G. Haemodynamics and Flow Modification Stents for Peripheral Arterial Disease: A Review. Ann. Biomed. Eng. 2016, 44, 466–476. [Google Scholar] [CrossRef] [Green Version]
  43. Arab, M.I.; Bouzit, M.; Ameur, H.; Kamla, Y. Numerical Study on Thoracic Aortic Aneurysms: The Aneurysm Aggravation Effects on the Secondary Flow Motion. Mechanics 2020, 26, 407–415. [Google Scholar] [CrossRef]
  44. Matsuzawa, T.; Gao, F.; Qiao, A.; Ohta, O.; Okada, H. Numerical Simulation in Aortic Arch Aneurysm. Etiol. Pathog. Pathophysiol. Aortic Aneurysms Aneurysm Rupture 2011, 12, 208–222. [Google Scholar] [CrossRef] [Green Version]
  45. Qiao, A.-K.; Fu, W.-Y.; Liu, Y.-J. Study on Hemodynamics in Patient-Specific Thoracic Aortic Aneurysm. Theor. Appl. Mech. Lett. 2011, 1, 14001. [Google Scholar] [CrossRef] [Green Version]
  46. Rolland, P.H.; Charifi, A.-B.; Verrier, C.; Bodard, H.; Friggi, A.; Piquet, P.; Moulin, G.; Bartoli, J.-M. Hemodynamics and Wall Mechanics after Stent Placement in Swine Iliac Arteries: Comparative Results from Six Stent Designs1. Radiology 1999, 213, 229–246. [Google Scholar] [CrossRef] [PubMed]
  47. Moore, J.E.; Berry, J.L. Fluid and Solid Mechanical Implications of Vascular Stenting. Ann. Biomed. Eng. 2002, 30, 498–508. [Google Scholar] [CrossRef] [PubMed]
  48. Nardi, A.; Avrahami, I. Approaches for Treatment of Aortic Arch Aneurysm, a Numerical Study. J. Biomech. 2017, 50, 158–165. [Google Scholar] [CrossRef]
  49. Romarowski, R.M.; Conti, M.; Morganti, S.; Grassi, V.; Marrocco-Trischitta, M.M.; Trimarchi, S.; Auricchio, F. Computational Simulation of TEVAR in the Ascending Aorta for Optimal Endograft Selection: A Patient-Specific Case Study. Comput. Biol. Med. 2018, 103, 140–147. [Google Scholar] [CrossRef]
  50. Qiao, Y.; Mao, L.; Zhu, T.; Fan, J.; Luo, K. Biomechanical Implications of the Fenestration Structure after Thoracic Endovascular Aortic Repair. J. Biomech. 2020, 99, 109478. [Google Scholar] [CrossRef]
  51. Mestres, G.; Garcia, M.E.; Yugueros, X.; Urrea, R.; Tripodi, P.; Gomez, F.; Maeso, J.; Riambau, V. Aortic Arch and Thoracic Aorta Curvature Remodeling after Thoracic Endovascular Aortic Repair. Ann. Vasc. Surg. 2017, 38, 233–241. [Google Scholar] [CrossRef]
  52. Hirotsu, K.; Suh, G.-Y.; Lee, J.T.; Dake, M.D.; Fleischmann, D.; Cheng, C.P. Changes in Geometry and Cardiac Deformation of the Thoracic Aorta after Thoracic Endovascular Aortic Repair. Ann. Vasc. Surg. 2018, 46, 83–89. [Google Scholar] [CrossRef]
  53. Fung, G.S.K.; Lam, S.K.; Cheng, S.W.K.; Chow, K.W. On Stent-Graft Models in Thoracic Aortic Endovascular Repair: A Computational Investigation of the Hemodynamic Factors. Comput. Biol. Med. 2008, 38, 484–489. [Google Scholar] [CrossRef]
Figure 1. Diagram of the MIBSG.
Figure 1. Diagram of the MIBSG.
Bioengineering 10 00164 g001
Figure 2. Top: Three-dimensional model reconstruction and mesh generation. (Case A) Normal geometry; (Case B) Preoperative geometry; (Case C) 1-week postoperative geometry; (Case D) 6-month postoperative geometry.
Figure 2. Top: Three-dimensional model reconstruction and mesh generation. (Case A) Normal geometry; (Case B) Preoperative geometry; (Case C) 1-week postoperative geometry; (Case D) 6-month postoperative geometry.
Bioengineering 10 00164 g002
Figure 3. Inlet and outlet boundary conditions.
Figure 3. Inlet and outlet boundary conditions.
Bioengineering 10 00164 g003
Figure 4. Vessel centerline (contains no branches). (A) Coronal plane view; (B) 90° anti-clockwise rotation of coronal plane view; (C) Horizontal plane view.
Figure 4. Vessel centerline (contains no branches). (A) Coronal plane view; (B) 90° anti-clockwise rotation of coronal plane view; (C) Horizontal plane view.
Bioengineering 10 00164 g004
Figure 5. Average perfusion rate of each outlet within one pulsation cycle.
Figure 5. Average perfusion rate of each outlet within one pulsation cycle.
Bioengineering 10 00164 g005
Figure 6. Intravascular flow field at peak systole. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model. The dashed circle is proximal to the third module and adjacent to the overlapping area of the branched stent and aortic stent. The red arrow indicates an area of disturbed blood flow along the lesser curvature of the aorta.
Figure 6. Intravascular flow field at peak systole. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model. The dashed circle is proximal to the third module and adjacent to the overlapping area of the branched stent and aortic stent. The red arrow indicates an area of disturbed blood flow along the lesser curvature of the aorta.
Bioengineering 10 00164 g006
Figure 7. The pressure distribution of the four models at peak systole. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model.
Figure 7. The pressure distribution of the four models at peak systole. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model.
Bioengineering 10 00164 g007
Figure 8. Changes in the regional distribution of the TAWSS of <0.4 Pa and the area statistics of the four models. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model.
Figure 8. Changes in the regional distribution of the TAWSS of <0.4 Pa and the area statistics of the four models. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model.
Bioengineering 10 00164 g008
Figure 9. Changes in the regional distribution of the OSI of >0.25 and the area statistics of the four models. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model.
Figure 9. Changes in the regional distribution of the OSI of >0.25 and the area statistics of the four models. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model.
Bioengineering 10 00164 g009
Figure 10. Changes in the regional distribution of the RRT of 5 Pa−1 and the area statistics of the four models. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model.
Figure 10. Changes in the regional distribution of the RRT of 5 Pa−1 and the area statistics of the four models. (Case A) Normal model; (Case B) Preoperative model; (Case C) 1-week postoperative model; (Case D) 6-week postoperative model.
Bioengineering 10 00164 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, F.; Zhu, Y.; Song, H.; Zhang, H.; Chen, L.; Guo, W. Analysis of Postoperative Remodeling Characteristics after Modular Inner Branched Stent-Graft Treatment of Aortic Arch Pathologies Using Computational Fluid Dynamics. Bioengineering 2023, 10, 164. https://doi.org/10.3390/bioengineering10020164

AMA Style

Li F, Zhu Y, Song H, Zhang H, Chen L, Guo W. Analysis of Postoperative Remodeling Characteristics after Modular Inner Branched Stent-Graft Treatment of Aortic Arch Pathologies Using Computational Fluid Dynamics. Bioengineering. 2023; 10(2):164. https://doi.org/10.3390/bioengineering10020164

Chicago/Turabian Style

Li, Fen, Yating Zhu, Hui Song, Hongpeng Zhang, Lingfeng Chen, and Wei Guo. 2023. "Analysis of Postoperative Remodeling Characteristics after Modular Inner Branched Stent-Graft Treatment of Aortic Arch Pathologies Using Computational Fluid Dynamics" Bioengineering 10, no. 2: 164. https://doi.org/10.3390/bioengineering10020164

APA Style

Li, F., Zhu, Y., Song, H., Zhang, H., Chen, L., & Guo, W. (2023). Analysis of Postoperative Remodeling Characteristics after Modular Inner Branched Stent-Graft Treatment of Aortic Arch Pathologies Using Computational Fluid Dynamics. Bioengineering, 10(2), 164. https://doi.org/10.3390/bioengineering10020164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop