3-D numerical simulation of blood flow through models of the human aorta

Effect of blood viscosity on the hemodynamics of arteriovenous fistulae based on numerical investigation

Arteriovenous fistula (AVF) is the most commonly used vascular access for hemodialysis in patients with end-stage renal disease. Vascular diseases such as atherosclerosis and thrombosis, triggered by altered hemodynamic conditions, are the main causes of access failure. Changes in blood viscosity accelerate access dysfunction by affecting local velocities and wall shear stress (WSS) distribution in the circulation. Numerical simulation was employed to analyze and compare the hemodynamic behavior of AVF under different blood viscosities (0.001-0.012 Pa∙s). An idealized three-dimensional model with end-to-side anastomosis was established. Transient simulations were conducted using pulsatile inlet velocity and outflow as boundary conditions. The simulation results reveal the blood flow state of AVF under different viscosity physiological conditions and derive the rule of change. When blood viscosity increases, the local velocity in the disturbed region slows down and the stagnation time becomes longer, resulting in increased deposition of substances. As blood viscosity increases, the level of shear stress on the entire wall of the fistula increases accordingly. WSS values at high viscosities above 0.007 Pa∙s showed significantly larger low-shear regions near the anastomosis and increased chances of inducing atheromatous plaques. This research has revealed the correlation between blood dynamic viscosity and the hemodynamic behavior of AVF. Elevated whole blood viscosity increases the incidence of access obstruction and vascular disease leading to fistula failure. The study provides a basis for optimizing the distribution of hemodynamic parameters in the fistula for hemodialysis patients.

Blood flow effects in a patient with a thoracic aortic endovascular prosthesis

This work analyzes hemodynamic phenomena within the aorta of two elderly patients and their impact on blood flow behavior, particularly affected by an endovascular prosthesis in one of them (Patient II). Computational Fluid Dynamics (CFD) was utilized for this study, involving measurements of velocity, pressure, and wall shear stress (WSS) at various time points during the third cardiac cycle, at specific positions within two cross sections of the thoracic aorta. The first cross-section (Cross-Section 1, CS1) is located before the initial fluid bifurcation, just before the right subclavian artery. The second cross-section (Cross-Section 2, CS2) is situated immediately after the left subclavian artery. The results reveal that, under regular aortic geometries, velocity and pressure magnitudes follow the principles of fluid dynamics, displaying variations. However, in Patient II, an endoprosthesis near the CS2 and the proximal border of the endoprosthesis significantly disrupts fluid behavior owing to the pulsatile flow. The cross-sectional areas of Patient I are smaller than those of Patient II, leading to higher flow magnitudes. Although in CS1 of Patient I, there is considerable variability in velocity magnitudes, they exhibit a more uniform and predictable transition. In contrast, CS2 of Patient II, where magnitude variation is also high, displays irregular fluid behavior due to the endoprosthesis presence. This cross-section coincides with the border of the fluid bifurcation. Additionally, the irregular geometry caused by endovascular aneurysm repair contributes to flow disruption as the endoprosthesis adjusts to the endothelium, reshaping itself to conform with the vessel wall. In this context, significant alterations in velocity values, pressure differentials fluctuating by up to 10%, and low wall shear stress indicate the pronounced influence of the endovascular prosthesis on blood flow behavior. These flow disturbances, when compounded by the heart rate, can potentially lead to changes in vascular anatomy and displacement, resulting in a disruption of the prosthesis-endothelium continuity and thereby causing clinical complications in the patient.

Toward a Mesoscopic Modeling Approach of Magnetohydrodynamic Blood Flow in Pathological Vessels: A Comprehensive Review

The investigation of magnetohydrodynamic (MHD) blood flow within configurations that are pertinent to the human anatomy holds significant importance in the realm of scientific inquiry because of its practical implications within the medical field. This article presents an exhaustive appraisal of the diverse applications of magnetohydrodynamics and their computational modeling in biological contexts. These applications are classified into two categories: simple flow and pulsatile flow. An alternative approach of traditional CFD methods called Lattice Boltzmann Method (LBM), a mesoscopic method based on kinetic theory, is introduced to solve complex problems, such as hemodynamics. The results show that the flow velocity reduces considerably by increasing the magnetic field intensity, and the flow separation area is minimized by the increase of magnetic field strength. The LBM with BGK collision model has shown good results in terms of precision. Finally, this literature review has revealed a number of potential avenues for further research. Suggestions for future works are proposed accordingly.

Near-wall hemodynamic changes in subclavian artery perfusion induced by retrograde inner branched thoracic endograft implantation

Objective

Left subclavian artery (LSA)-branched endografts with retrograde inner branch configuration (thoracic branch endoprosthesis [TBE]) offer a complete endovascular solution when LSA preservation is required during zone 2 thoracic endovascular aortic repair. However, the hemodynamic consequences of the TBE have not been well-investigated. We compared near-wall hemodynamic parameters before and after the TBE implantation using computational fluid dynamic simulations.

Methods

Eleven patients who had undergone TBE implantation were included. Three-dimensional aortic arch geometries were constructed from the pre- and post-TBE implantation computed tomography images. The resulting 22 three-dimensional aortic arch geometries were then discretized into finite element meshes for computational fluid dynamic simulations. Inflow boundary conditions were prescribed using normal physiological pulsatile circulation. Outlet boundary conditions consisted of Windkessel models with previously published values. Blood flow, modeled as Newtonian fluid, simulations were performed with rigid wall assumptions using SimVascular's incompressible Navier-Stokes solver. We compared well-established hemodynamic descriptors: pressure, flow rate, time-averaged wall shear stress (TAWSS), the oscillatory shear index (OSI), and percent area with an OSI of >0.2. Data were presented on the stented portion of the LSA.

Results

TBE implantation was associated with a small decrease in peak LSA pressure (153 mm Hg; interquartile range [IQR], 151-154 mm Hg vs 159 mm Hg; IQR, 158-160 mm Hg; P = .005). No difference was observed in peak LSA flow rates before and after implantation: 40.4 cm³/ (IQR, 39.5-41.6 cm³/s) vs 41.3 cm³/s (IQR, 37.2-44.8 cm³/s; P = .59). There was a significant postimplantation increase in TAWSS (15.2 dynes/cm² [IQR, 12.2-17.7 dynes/cm²] vs 6.2 dynes/cm² [IQR, 5.7-10.3 dynes/cm²]; P = .003), leading to decreases in both the OSI (0.088 [IQR, 0.063 to -0.099] vs 0.1 [IQR, 0.096-0.16]; P = .03) and percentage of area with an OSI of >0.2 (10.4 [IQR, 5.8-15.8] vs 15.7 [IQR, 10.7-31.9]; P = .13). Neither LSA side branch angulation (median, 81°, IQR, 77°-109°) nor moderate compression (16%-58%) seemed to have an impact on the pressure, flow rate, TAWSS, or percentage of area with an OSI of >0.2 in the stented LSA.

Conclusions

The implantation of TBE produces modest hemodynamic disturbances that are unlikely to result in clinically relevant changes.

Analysis of blood flow of unsteady Carreau-Yasuda nanofluid with viscous dissipation and chemical reaction under variable magnetic field

Blood flow analysis through arterial walls depicts unsteady non-Newtonian fluid flow behavior. Arterial walls are impacted by various chemical reactions and magnetohydrodynamic effects during treatment of malign and tumors, cancers, drug targeting and endoscopy. In this regard, current manuscript focuses on modeling and analysis of unsteady non-Newtonian Carreau-Yasuda fluid with chemical reaction, Brownian motion and thermophoresis under variable magnetic field. The main objective is to simulate the effect of different fluid parameters, especially variable magnetic field, chemical reaction and viscous dissipation on the blood flow to help medical practitioners in predicting the changes in blood to make diagnosis and treatment more efficient. Suitable similarity transformations are used for the conversion of partial differential equations into a coupled system of ordinary differential equations. Homotopy analysis method is used to solve the system and convergent results are drawn. Effect of different dimensionless parameters on the velocity, temperature and concentration profiles of blood flow are analyzed in shear thinning and thickening cases graphically. Analysis reveals that chemical reaction increases blood concentration which enhance the drug transportation. It is also observed that magnetic field elevates the blood flow in shear thinning and thickening scenarios. Furthermore, Brownian motion and thermophoresis increases temperature profile.

Comparison of thrombosis risk in an abdominal aortic dissection aneurysm with a double false lumen using computational fluid dynamic simulation method

BACKGROUND: Aneurysms are associated with a mortality rate of 81% or more in cases of rupture. OBJECTIVE: To analyse the haemodynamic indices and compare the thrombosis risk in a double false lumen abdominal aortic dissection aneurysm using computational fluid dynamics (CFD). METHODS: Computer tomography angiography (CTA) imaging data were collected from a patient with a double false lumen abdominal aortic dissection aneurysm, and three different lesion morphology aneurysm models were established, double false lumen abdominal aortic dissection aneurysm, single false lumen abdominal aortic dissection aneurysm and saccular abdominal aortic aneurysm, in order to analyse the flow velocity, time-averaged shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT) of blood flow, and endothelial cell activation potential (ECAP). RESULTS: All three aneurysms were in a low-flow state within the body, and the low-flow velocity flow in the proximal vessel wall extended to the right common iliac artery; the vortex intensity was more intense in the abdominal aortic dissection aneurysm than in the saccular abdominal aortic aneurysm. The risk area for thrombosis was concentrated in the expansion part of the aneurysm and the false lumen. The RRT and ECAP maxima of the double false lumen abdominal aortic dissection aneurysm were much greater than those of the single false lumen dissection aneurysm and saccular aortic aneurysm. CONCLUSION: Low-velocity blood flow, high OSI, low TAWSS, high RRT, and high ECAP regions correlate with the risk of thrombosis. The double false lumen type of abdominal aortic dissection aneurysm had some specificity in this case. The risk of thrombosis in the patient was extremely high, and the largest risk zone was within the smaller false lumen, which could be because the smaller false lumen was connected to the true lumen by only one breach. The results of the study provide some guidance in the early screening and development of treatment plans.

Non-invasive diagnostics of blockage growth in the descending aorta-computational approach

Atherosclerosis causes blockages to the main arteries such as the aorta preventing blood flow from delivering oxygen to the organs. Non-invasive diagnosis of these blockages is difficult, particularly in primary healthcare. In this paper, the effect of arterial blockage development and growth is investigated at the descending aorta on some possible non-invasive assessment parameters including the blood pressure waveform, wall shear stress (WSS), time-average WSS (TAWSS) and the oscillation shear index (OSI). Blockage severity growth is introduced in a simulation model as 25%, 35%, 50% and 65% stenosis at the descending aorta based on specific healthy control aorta data clinically obtained. A 3D aorta model with invasive pulsatile waveforms (blood flow and pressure) is used in the CFD simulation. Blockage severity is assessed by using blood pressure measurements at the left subclavian artery. An arterial blockage growth more than 35% of the lumen diameter significantly affects the pressure. A strong correlation is also observed between the ascending aorta pressure values, pressure at the left subclavian artery and the relative residence time (RRT). An increase of RRT downstream from the stenosis indicates a 35% stenosis at the descending aorta which results in high systolic and diastolic pressure readings. The findings of this study could be further extended by transferring the waveform reading from the left subclavian artery to the brachial artery.

Graphical abstract

Hemodynamic numerical simulation of aortic arch modular inner branched stent-graft in eight early patients from the first-in-human case series

Background

The modular inner branched stent-graft (MIBSG) (WeFlow-Arch™) is an emerging device for challenging aortic arch pathologies. Hemodynamic numerical simulation is conducive to predicting long-term outcomes as well as optimizing the stent-graft design.

Objective

This study aims to analyze the hemodynamic characteristics of the MIBSG devices based on numerical simulation analyses.

Methods

From June 2019 to June 2021, MIBSGs were utilized in eight cases. Numerical simulation analyses of branch perfusion and indicators including the time-averaged wall shear stress, oscillatory shear index, and relative residence time were performed.

Results

Lesions involved Zone 1 (n = 2), Zone 2 (n = 4), and Zone 3 (n = 2). Branched stent-grafts were deployed in the innominate artery and left common carotid artery (n = 5) or in the innominate artery and left subclavian artery (n = 3). The hemodynamic change in common was increased perfusion in the descending aorta and left common carotid artery. Half of the patients had increased cerebral perfusion of 8.7% at most, and the other half of the patients showed a reduction of 5.3% or less. Case 3 was considered to have acquired the greatest improvement in hemodynamic features.

Conclusion

The MIBSG showed improved hemodynamic features in most cases. The design of the MIBSG could be partly modified to acquire better hemodynamic performance.

The Dean Effect: An Aortic Arch Flow Artifact Mimicking Dissection

The unique hemodynamics of the aortic arch create conditions for potential formation of a flow-related artifact that mimics disease on CT angiographic images. The hemodynamic basis for this artifact can be explained by fluid mechanics incorporating a mathematical principle known as the Dean number. Therefore, in this review, the artifact is referred to as the Dean effect. It is important for radiologists and other clinicians to recognize this artifact when encountered. It is also helpful for the interpreting radiologist to have a basic understanding of the relevant hemodynamic principles. This review provides an example of the artifact, reviews the basic underlying hemodynamics, and presents methods of how to prevent this artifact and distinguish it from pathologic mimics in clinical practice. Keywords: CT Angiography, Vascular, Thorax, Aorta, Artifacts, Blood, Dissection, Hemodynamics/Flow Dynamics © RSNA, 2022.

Numerical investigation of patient-specific thoracic aortic aneurysms and comparison with normal subject via computational fluid dynamics (CFD)

Vascular hemodynamics play an important role in cardiovascular diseases. This work aimed to investigate the effects of an increase in ascending aortic diameter (AAD) on hemodynamics throughout a cardiac cycle for real patients. In this study, two scans of thoracic aortic aneurysm (TAA) subject with different AADs (42.94 mm and 48.01 mm) and a scan of a normal subject (19.81 mm) were analyzed to assess the effects of hemodynamics on the progression of TAA with the same flow rate. Real-patient aortic geometries were scanned by computed tomography angiography (CTA), and steady and pulsatile flow conditions were used to simulate real patient aortic geometries. Aortic arches were obtained from routine clinical scans. Computational fluid dynamics (CFD) simulations were performed with in vivo boundary conditions, and 3D Navier-Stokes equations were solved by a UDF (user-defined function) code defining a real cardiac cycle of one patient using Fourier series (FS). Wall shear stress (WSS) and pressure distributions were presented from normal subject to TAA cases. The results show that during the peak systolic phase pressure load increased by 18.56% from normal subject to TAA case 1 and by 23.8% from normal subject to TAA case 2 in the aneurysm region. It is concluded that although overall WSS increased in aneurysm cases but was low in dilatation areas. As a result, abnormal changes in WSS and higher pressure load may lead to rupture and risk of further dilatation. CFD simulations were highly effective to guide clinical predictions and assess the progress of aneurysm regions in case of early surgical intervention. Graphical abstract.

Analysis and manegement of laminar blood flow inside a cerebral blood vessel using a finite volume software program for biomedical engineering

BACKGROUND AND OBJECTIVE

Hemodynamic blood flow analysis in the cerebrovascular is has become one of the important research topics in the bio-mechanic in recent decades. The primary duty of the cerebral blood vessel is supplying Glucose and oxygen for the brain.

METHODS

In this investigation, the non-Newtonian blood flow in the cerebral blood vessels studied. For modeling the geometry of this problem, we used Magnetic Resonance Image (MRI) approach to take Digital Imaging and Communications in Medicine (DICOM) images and using an open-source software package to construct the geometry, which is a complicated one. The power-law indexes, heat flux, and Reynolds number range in the investigation are 0.6 ≤ n ≤ 0.8, 5 ≤ q ≤ 15Wm^-2 and 160≤Re≤310. Effects of Reynolds number, power-law indexes and heat fluxes are investigated.

RESULTS

We found that the pressure drop increase with increasing the Reynolds number and power-law index. The maximum Nusselt number in the cerebral blood vessels accrued in the running position of the body in n = 0.8. Also, the highest average wall shear stress occurs in maximum power-law indexes and Reynolds number.

CONCLUSION

By increasing the power-law index and Reynolds number, the wall shear stress increases.

The application of computational modeling for risk prediction in type B aortic dissection

OBJECTIVE

New tools are urgently needed to help with surgical decision-making in type B aortic dissection (TBAD) that is uncomplicated at the time of initial presentation. This narrative review aims to answer the clinical question, Can computational modeling be used to predict risk in acute and chronic Stanford TBAD?

METHODS

The review (PROSPERO 2018 CRD42018104472) focused on risk prediction in TBAD. A comprehensive search of the Ovid MEDLINE database, using terms related to computational modeling and aortic dissection, was conducted to find studies of any form published between 1998 and 2018. Cohort studies, case series, and case reports of adults (older than 18 years) with computed tomography or magnetic resonance imaging diagnosis of TBAD were included. Computational modeling was applied in all selected studies.

RESULTS

There were 37 studies about computational modeling of TBAD identified from the search, and the findings were synthesized into a narrative review. Computational modeling can produce numerically calculated values of stresses, pressures, and flow velocities that are difficult to measure in vivo. Hemodynamic parameters-high or low wall shear stress, high pressure gradient between lumens during the cardiac cycle, and high false lumen flow rate-have been linked to the pathogenesis of branch malperfusion and aneurysm formation by numerous studies. Considering the major outcomes of end-organ failure, aortic rupture, and stabilization and remodeling, hypotheses have been generated about inter-relationships of measurable parameters in computational models with observable anatomic and pathologic changes, resulting in specific clinical outcomes.

CONCLUSIONS

There is consistency in study findings about computational modeling in TBAD, although a limited number of patients have been analyzed using various techniques. The mechanistic patterns of association found in this narrative review should be investigated in larger cohort prospective studies to further refine our understanding. It highlights the importance of patient-specific computational hemodynamic parameters in clinical decision-making algorithms. The current challenge is to develop and to test a risk assessment method that can be used by clinicians for TBAD.

Numerical Investigation of Blood Flow Characteristics through Cannulated Aorta

Cannulation of the aorta is done in order to provide oxygenation and circulatory function through the use of the heart lung machine during cardio-pulmonary bypass (CPB). The nature of the blood flow through the aorta and its ramifications during CPB is mostly linear as compared to the physiological flow, which is pulsatile in nature. This leads to the development of multiple morbidities caused by the development of emboli and atheromas. Perioperative postoperative care is necessitated by these conditions. As such the understanding of the blood flow characteristics is necessitated in order to effectively prevent the formation of emboli and to prevent the "Sandblasting" effect. The authors in this work seek to investigate the nature of blood flow through the aorta under such circumstances. The results obtained show the nature of blood flow in the cannulated aorta as well as the optimum angle of placement of the cannula with respect to the aortic wall.

Hemodynamics in diabetic human aorta using computational fluid dynamics

Three-dimensional (3D) computational aortic models have been established to reproduce aortic diseases such as aortic aneurysm and dissection; however, no such models have been developed to study diabetes mellitus (DM). To characterize biomechanical properties of the human aorta with DM, reconstructed aortic CT images were converted into DICOM format, and imported into the 3D segmentation using Mimics software. This resulted in a 3D reconstruction of the complete aorta, including three branches. We applied a pulsatile blood pressure waveform for the ascending aorta to provide a biomimetic environment using COMSOL Multiphysics software. Hemodynamics were compared between the control and DM models. We observed that mean blood flow velocity, aortic pressure, and von Mises stress values were lower in the DM model than in the control model. Furthermore, the range of aortic movement was lower in the DM model than in the control model, suggesting that the DM aortic wall is more susceptible to rupture. When comparing biomechanical properties in discrete regions of the aorta, all values were higher in the ascending aorta for both control and DM models, corresponding to the location of most aortic lesions. We have developed a compute based that integrates advanced image processing strategies and computational techniques based on finite element method to perform hemodynamics analysis based on CT images. Our study of image-based CFD analysis hopes to provide a better understanding of the relationship between aortic hemodynamic and developing pathophysiology of aortic diseases.

Comparison between a generalized Newtonian model and a network-type multiscale model for hemodynamic behavior in the aortic arch: Validation with 4D MRI data for a case study

Blood is a complex fluid in which the presence of the various constituents leads to significant changes in its rheological properties. Thus, an appropriate non-Newtonian model is advisable; and we choose a Modified version of the rheological model of Phan-Thien and Tanner (MPTT). The different parameters of this model, derived from the rheology of polymers, allow characterization of the non-Newtonian nature of blood, taking into account the behavior of red blood cells in plasma. Using the MPTT model that we implemented in the open access software OpenFOAM, numerical simulations have been performed on blood flow in the thoracic aorta for a healthy patient. We started from a patient-specific model which was constructed from medical images. Exiting flow boundary conditions have been developped, based on a 3-element Windkessel model to approximate physiological conditions. The parameters of the Windkessel model were calibrated with in vivo measurements of flow rate and pressure. The influence of the selected viscosity of red blood cells on the flow and wall shear stress (WSS) was investigated. Results obtained from this model were compared to those of the Newtonian model, and to those of a generalized Newtonian model, as well as to in vivo dynamic data from 4D MRI during a cardiac cycle. Upon evaluating the results, the MPTT model shows better agreement with the MRI data during the systolic and diastolic phases than the Newtonian or generalized Newtonian model, which confirms our interest in using a complex viscoelastic model.

Modeling and analysis of biomagnetic blood Carreau fluid flow through a stenosis artery with magnetic heat transfer: A transient study

We present a numerical investigation of tapered arteries that addresses the transient simulation of non-Newtonian bio-magnetic fluid dynamics (BFD) of blood through a stenosis artery in the presence of a transverse magnetic field. The current model is consistent with ferro-hydrodynamic (FHD) and magneto-hydrodynamic (MHD) principles. In the present work, blood in small arteries is analyzed using the Carreau-Yasuda model. The arterial wall is assumed to be fixed with cosine geometry for the stenosis. A parametric study was conducted to reveal the effects of the stenosis intensity and the Hartman number on a wide range of flow parameters, such as the flow velocity, temperature, and wall shear stress. Current findings are in a good agreement with recent findings in previous research studies. The results show that wall temperature control can keep the blood in its ideal blood temperature range (below 40°C) and that a severe pressure drop occurs for blockages of more than 60 percent. Additionally, with an increase in the Ha number, a velocity drop in the blood vessel is experienced.

Computer-Simulated Fluid Dynamics of Arterial Perfusion in Extracorporeal Circulation: From Reality to Virtual Simulation

Background

Atheroembolism due to aortic manipulation remains an unsolved problem in surgery for thoracic aortic aneurysm. The goal of the present study is to create a computer simulation (CS) model with which to analyze blood flow in the diseased aorta.

Method

A three-dimensional glass model of the aortic arch was constructed from CT images of a normal, healthy person and a patient with transverse aortic arch aneurysm. Separately, a CS model of the curved end-hole cannula was created, and flow from the aortic cannula was recreated using a numerical simulation.

Results

Comparison of the data obtained by the glass model analyses revealed that the flow velocity and the vector of the flow around the exit of the cannula were similar to that in the CS model. A high-velocity area was observed around the cannula exit in both the glass model and the CS model. The maximum flow velocity was as large as 1.0 m/s at 20 mm from the cannula exit and remained as large as 0.5 to 0.6 m/s within 50 mm of the exit. In the aortic arch aneurysm models, the rapid jet flow from the cannula moved straight toward the lesser curvature of the transverse aortic arch. The locations and intensities of the calculated vortices were slightly different from those obtained for the glass model.

Conclusions

The proposed CS method for the analysis of blood flow from the aortic cannulae during extracorporeal circulation can reproduce the flow velocity and flow pattern in the proximal and transverse aortic arches.

Red blood cell motion and deformation in a curved microvessel

The flow of cells through curved vessels is often encountered in various biomedical and bioengineering applications, such as red blood cells (RBCs) passing through the curved arteries in circulation, and cells sorting through a shear-induced migration in a curved channels. Most of past numerical studies focused on the cell deformation in small straight microvessels, or on the flow pattern in large curved vessels without considering the cell deformation. However, there have been few attempts to study the cell deformation and the associated flow pattern in a curved microvessel. In this work, a particle-based method, smoothed dissipative particle dynamics (SDPD), is used to simulate the motion and deformation of a RBC in a curved microvessel of diameter comparable to the RBC diameter. The emphasis is on the effects of the curvature, the type and the size of the curved microvessel on the RBC deformation and the flow pattern. The simulation results show that a small curved shape of the microvessel has negligible effect on the RBC behavior and the flow pattern which are similar to those in a straight microvessel. When the microvessel is high in curvature, the secondary flow comes into being with a pair of Dean vortices, and the velocity profile of the primary flow is skewed toward the inner wall of the microvessel. The RBC also loses the axisymmetric deformation, and it is stretched first and then shrinks when passing through the curved part of the microvessel with the large curvature. It is also found that a pair of Dean vortices arise only under the condition of De>1 (De is the Dean number, a ratio of centrifugal to viscous competition). The Dean vortices are more easily observed in the larger or more curved microvessels. Finally, it is observed that the velocity profile of primary flow is skewed toward the inner wall of curved microvessel, i.e., the fluid close to the inner wall flows faster than that close to the outer wall. This is contrary to the common sense in large curved vessels. This velocity skewness was found to depend on the curvature of the microvessel, as well as the viscous and inertial forces.

Elimination of Transcoarctation Pressure Gradients Has No Impact on Left Ventricular Function or Aortic Shear Stress After Intervention in Patients With Mild Coarctation

OBJECTIVES

This study sought to investigate the impact of transcatheter intervention on left ventricular function and aortic hemodynamics in patients with mild coarctation of the aorta (COA).

BACKGROUND

The optimal method and timing of transcatheter intervention for COA remains unclear, especially when the severity of COA is mild (peak-to-peak transcoarctation pressure gradient <20 mm Hg). Debate rages regarding the risk/benefit ratio of intervention versus long-term effects of persistent minimal gradient in this heterogeneous population with differing blood pressures, ventricular function, and peripheral perfusion.

METHODS

We developed a unique computational fluid dynamics and lumped parameter modeling framework based on patient-specific hemodynamic input parameters and validated it against patient-specific clinical outcomes (before and after intervention). We used clinically measured hemodynamic metrics and imaging of the aorta and the left ventricle in 34 patients with mild COA to make these correlations.

RESULTS

Despite dramatic reduction in the transcoarctation pressure gradient (catheter and Doppler echocardiography pressure gradients reduced by 75% and 47.3%, respectively), there was only modest effect on aortic flow and no significant impact on aortic shear stress (the maximum time-averaged wall shear stress in descending aorta was reduced 5.1%). In no patient did transcatheter intervention improve left ventricular function (e.g., stroke work and normalized stroke work were reduced by only 4.48% and 3.9%, respectively).

CONCLUSIONS

Transcatheter intervention that successfully relieves mild COA pressure gradients does not translate to decreased myocardial strain. The effects of the intervention were determined to the greatest degree by ventricular-vascular coupling hemodynamics and provide a novel valuable mechanism to evaluate patients with COA that may influence clinical practice.

The risk of stanford type-A aortic dissection with different tear size and location: a numerical study

Background

This study is to investigate the influence of hemodynamics on Stanford type-A aortic dissection with different tear size and location, to provide some support for the relationships between the risks (rupture, reverse tearing and further tearing) and tear size and location for clinical treatment.

Methods

Four numerical models of Stanford type-A aortic dissection were established, with different size and location of the tears. The ratio of the area between the entry and re-entry tears(RA) is various within the model; while, the size and the location of the re-entry in the distal descending aorta are fixed. In model A11 and A21, the entry tears are located near the ascending aorta. The RA in these models are 1 and 2, respectively; In the model B11 and B21, the entry tears are located near the proximal descending aorta and the RA in these models are again assigned to 1 and 2, respectively. Then hemodynamics in these models was solved with numerically and the flow patterns and loading distributions were investigated.

Results

The flow velocity of the true lumen in model A21, B21 is lower than that in A11, B11, respectively; the time-averaged wall shear stress (TAWSS) of the false lumen in model A21 and B21 is higher, and for ascending aorta false lumen, A11, A21 are higher than B11, B21, respectively. False lumen intimal wall pressure of A11, A21 are always higher than the true lumen ones.

Conclusion

The variation of the RA can significantly affect the dynamics of blood within the aortic dissection. When the entry tear size is larger than the re-entry tear ones, the false lumen, proximal descending aorta and the wall near re-entry tear are prone to cracking. Entry tear location can significantly alter the hemodynamics of aortic dissection as well. When entry tear location is closer to proximal ascending aorta, false lumen continues to expand and compress the true lumen resulting in the true lumen reduction. For proximal ascending aorta, high pressure in false lumen predicts a higher risk of reverse tear.

Laminar-to-turbulence and relaminarization zones detection by simulation of low Reynolds number turbulent blood flow in large stenosed arteries

Laminar, turbulent, transitional, or combine areas of all three types of viscous flow can occur downstream of a stenosis depending upon the Reynolds number and constriction shape parameter. Neither laminar flow solver nor turbulent models for instance the k-ω (k-omega), k-ε (k-epsilon), RANS or LES are opportune for this type of flow. In the present study attention has been focused vigorously on the effect of the constriction in the flow field with a unique way. It means that the laminar solver was employed from entry up to the beginning of the turbulent shear flow. The turbulent model (k-ω SST Transitional Flows) was utilized from starting of turbulence to relaminarization zone while the laminar model was applied again with onset of the relaminarization district. Stenotic flows, with 50 and 75% cross-sectional area, were simulated at Reynolds numbers range from 500 to 2000 employing FLUENT (v6.3.17). The flow was considered to be steady, axisymmetric, and incompressible. Achieving results were reported as axial velocity, disturbance velocity, wall shear stress and the outcomes were compared with previously experimental and CFD computations. The analogy of axial velocity profiles shows that they are in acceptable compliance with the empirical data. As well as disturbance velocity and wall shear stresses anticipated by this new approach, part by part simulation, are reasonably valid with the acceptable experimental studies.

Numerical simulation of blood flow in femoral perfusion: comparison between side-armed femoral artery perfusion and direct femoral artery perfusion

Computational numerical analysis was performed to elucidate the flow dynamics of femoral artery perfusion. Numerical simulation of blood flow was performed from the right femoral artery in an aortic model. An incompressible Navier-Stokes equation and continuity equation were solved using computed flow dynamics software. Three different perfusion models were analyzed: a 4.0-mm cannula (outer diameter 15 French size), a 5.2-mm cannula (18 French size) and an 8-mm prosthetic graft. The cannula was inserted parallel to the femoral artery, while the graft was anastomosed perpendicular to the femoral artery. Shear stress was highest with the 4-mm cannula (172 Pa) followed by the graft (127 Pa) and the 5.2-mm cannula (99 Pa). The cannula exit velocity was high, even when the 5.2-mm cannula was used. Although side-armed perfusion with an 8-mm graft generated a high shear stress area near the point of anastomosis, flow velocity at the external iliac artery was decreased. The jet speed decreased due to the Coanda effect caused by the recirculation behind sudden expansion of diameter, and the flow velocity maintains a constant speed after the reattachment length of the flow. This study showed that iliac artery shear stress was lower with the 5.2-mm cannula than with the 4-mm cannula when used for femoral perfusion. Side-armed graft perfusion generates a high shear stress area around the anastomotic site, but flow velocity in the iliac artery is slower in the graft model than in the 5.2-mm cannula model.

An experimental assessment of catheter trackability forces with tortuosity parameters along patient-specific coronary phantoms

Coronary artery disease is one of the leading causes of cardiovascular deaths worldwide. Approximately 70% of patients requiring coronary revascularisation receive endovascular stents. The endovascular procedure is the preferred option due to its minimally invasive nature when compared to open heart surgery. Stent delivery is paramount for the success of the endovascular procedure. Catheter delivery forces within tortuous blood vessels can produce vasoconstriction and injury, resulting in reactive intimal proliferation or distal embolisation associated with end-organ ischaemia and infarction. Trackability is evaluated by most medical device companies for further development of their delivery systems. Relevant device design attributes must be tested in settings which simulate aspects of the intended use conditions, such as vessel geometry and compliance. Various tortuosity parameters are used to facilitate endovascular intervention planning. This study assessed the significance and correlation between the trackability forces for a coronary stent system with various geometrical parameters based on patient-specific geometries. A motorised delivery system delivered a commercially available coronary stent system and monitored the trackability forces along three phantom patient-specific thin-walled, compliant coronary vessels supported by a cardiac phantom model. The maximum trackability forces, curvature and torsion values ranged from 0.31 to 0.87 N, 0.06 to 0.22 mm(-1) and -11.1 to 5.8 mm(-1), respectively. The trackability forces were significantly different between all vessels (p < 0.002), while the tortuosity parameters were not significantly different (p > 0.05). A new tortuosity parameter-coined tracking curvature which considers the lumen radius as well as the curvature along the centreline was statistically different (p < 0.002) for all vessels and correlated with the trackability forces. There was a strong correlation between the cumulative trackability force and the cumulative tracking curvature. Tracking curvature could be used as a predictive clinical tool to aid stent delivery to the vicinity of the lesion.

Effects of arterial blood flow on walls of the abdominal aorta: distributions of wall shear stress and oscillatory shear index determined by phase-contrast magnetic resonance imaging

Although abdominal aortic aneurysms (AAAs) occur mostly inferior to the renal artery, the mechanism of the development of AAA in relation to its specific location is not yet clearly understood. The objective of this study was to evaluate the hypothesis that even healthy volunteers may manifest specific flow characteristics of blood flow and alter wall shear or oscillatory shear stress in the areas where AAAs commonly develop. Eight healthy male volunteers were enrolled in this prospective study, aged from 24 to 27. Phase-contrast magnetic resonance imaging (MRI) was performed with electrocardiographic triggering. Flow-sensitive four-dimensional MR imaging of the abdominal aorta, with three-directional velocity encoding, including simple morphological image acquisition, was performed. Information on specific locations on the aortic wall was applied to the flow encodes to calculate wall shear stress (WSS) and oscillatory shear index (OSI). While time-framed WSS showed the highest peak of 1.14 ± 0.25 Pa in the juxtaposition of the renal artery, the WSS plateaued to 0.61 Pa at the anterior wall of the abdominal aorta. The OSI peaked distal to the renal arteries at the posterior wall of the abdominal aorta of 0.249 ± 0.148, and was constantly elevated in the whole abdominal aorta at more than 0.14. All subjects were found to have elevated OSI in regions where AAAs commonly occur. These findings indicate that areas of constant peaked oscillatory shear stress in the infra-renal aorta may be one of the factors that lead to morphological changes over time, even in healthy individuals.

The Effects That Cardiac Motion has on Coronary Hemodynamics and Catheter Trackability Forces for the Treatment of Coronary Artery Disease: An In Vitro Assessment

The coronary arterial tree experiences large displacements due to the contraction and expansion of the cardiac muscle and may influence coronary haemodynamics and stent placement. The accurate measurement of catheter trackability forces within physiological relevant test systems is required for optimum catheter design. The effects of cardiac motion on coronary flowrates, pressure drops, and stent delivery has not been previously experimentally assessed. A cardiac simulator was designed and manufactured which replicates physiological coronary flowrates and cardiac motion within a patient-specific geometry. A motorized delivery system delivered a commercially available coronary stent system and monitored the trackability forces along three phantom patient-specific thin walled compliant coronary vessels supported by a dynamic cardiac phantom model. Pressure drop variation is more sensitive to cardiac motion than outlet flowrates. Maximum pressure drops varied from 7 to 49 mmHg for a stenosis % area reduction of 56 to 90%. There was a strong positive linear correlation of cumulative trackability force with the cumulative curvature. The maximum trackability forces and curvature ranged from 0.24 to 0.87 N and 0.06 to 0.22 mm(-1) respectively for all three vessels. There were maximum and average percentage differences in trackability forces of (23-49%) and (1.9-5.2%) respectively when comparing a static pressure case with the inclusion of pulsatile flow and cardiac motion. Cardiac motion with pulsatile flow significantly altered (p value <0.001) the trackability forces along the delivery pathways with high local percentage variations and pressure drop measurements.

Effects of severity and location of stenosis on the hemodynamics in human aorta and its branches

Pulsatile blood flow is studied in a three-dimensional model of human thoracic aorta at different stages of atherosclerotic lesion growth, taking into account the effect of atherosclerotic plaque location and peripheral symmetry. The model is reconstructed from the computed tomography images. The wall shear stress (WSS), time-averaged WSS, and the oscillatory shear index are applied to determine susceptible sites for the onset of early atherosclerosis. Then, two different degrees of stenosis severity, 50 and 80 %, are introduced to vulnerable areas of the healthy aorta geometry. The overriding issue addressed is that the WSS distribution and magnitude are strongly affected by the atherosclerotic plaque size, its symmetric features, and the location, i.e., the branch it is formed. The present study, for the first time, is capable of providing information on the high shear environment that may exist upon the rupture of plaque surface and any thrombosis due to platelet deposition. The magnitude of WSS and its distribution at the throat of 50 % stenosed aortic arch are in agreement with the previous numerical study (Huang et al. in Exp Fluids 48(3):497-508, 2010). Results show that WSS values exceed 50 Pa at the throat of 80 % stenosed left common carotid and brachiocephalic arteries.

The influence of computational assumptions on analysing abdominal aortic aneurysm haemodynamics

The variation in computational assumptions for analysing abdominal aortic aneurysm haemodynamics can influence the desired output results and computational cost. Such assumptions for abdominal aortic aneurysm modelling include static/transient pressures, steady/transient flows and rigid/compliant walls. Six computational methods and these various assumptions were simulated and compared within a realistic abdominal aortic aneurysm model with and without intraluminal thrombus. A full transient fluid-structure interaction was required to analyse the flow patterns within the compliant abdominal aortic aneurysms models. Rigid wall computational fluid dynamics overestimates the velocity magnitude by as much as 40%-65% and the wall shear stress by 30%-50%. These differences were attributed to the deforming walls which reduced the outlet volumetric flow rate for the transient fluid-structure interaction during the majority of the systolic phase. Static finite element analysis accurately approximates the deformations and von Mises stresses when compared with transient fluid-structure interaction. Simplifying the modelling complexity reduces the computational cost significantly. In conclusion, the deformation and von Mises stress can be approximately found by static finite element analysis, while for compliant models a full transient fluid-structure interaction analysis is required for acquiring the fluid flow phenomenon.

An innovative numerical approach to resolve the pulse wave velocity in a healthy thoracic aorta model

Aortic dissection and atherosclerosis are highly fatal diseases. The development of both diseases is closely associated with highly complex haemodynamics. Thus, in predicting the onset of cardiac disease, it is desirable to obtain a detailed understanding of the flowfield characteristics in the human cardiovascular circulatory system. Accordingly, in this study, a numerical model of a normal human thoracic aorta is constructed using the geometry information obtained from a phase-contrast magnetic resonance imaging (PC-MRI) technique. The interaction between the blood flow and the vessel wall dynamics is then investigated using a coupled fluid-structure interaction (FSI) analysis. The simulations focus specifically on the flowfield characteristics and pulse wave velocity (PWV) of the blood flow. Instead of using a conventional PC-MRI method to measure PWV, we present an innovative application of using the FSI approach to numerically resolve PWV for the assessment of wall compliance in a thoracic aorta model. The estimated PWV for a normal thoracic aorta agrees well with the results obtained via PC-MRI measurement. In addition, simulations which consider the FSI effect yield a lower predicted value of the wall shear stress at certain locations in the cardiac cycle than models which assume a rigid vessel wall. Consequently, the model provides a suitable basis for the future development of more sophisticated methods capable of performing the computer-aided analysis of aortic blood flows.

THE IMPACT OF THE NUMBER OF TEARS IN PATIENT-SPECIFIC STANFORD TYPE B AORTIC DISSECTING ANEURYSM: CFD SIMULATION

It is believed that the progression of Stanford type B aortic dissection is closely associated with vascular geometry and hemodynamic parameters. The hemodynamic differences owing to the presence of greater than two tears have not been explored. The focus of the present study is to investigate the impact of an additional re-entry tear on the flow, pressure and wall shear stress distribution in the dissected aorta. A 3D aorta model with one entry and one re-entry tear was generated from computed tomography (CT) angiographic images of a patient with Stanford Type B aortic dissection. To investigate the hemodynamic effect of more than two tear locations, an additional circular re-entry tear was added 24 mm above the original re-entry tear. Our simulation results showed that the presence of an additional re-entry tear provided an extra return path for blood back to the true lumen during systole, and an extra outflow path into the false lumen during diastole. The presence of this additional path led to a decrease in the false lumen pressure, particularly at the distal region. Meanwhile, the presence of this additional tear causes no significant difference on the time average wall shear stress (TAWSS) distribution except at regions adjacent to re-entry tear 2. Moderate and concentrated TAWSS was observed at the bottom region of this additional tear which may lead to further extension of the tear distally.

Investigation of risks for cerebral embolism associated with the hemodynamics of cardiopulmonary bypass cannula: a numerical model

Cerebral emboli originating in the ascending aorta are a major cause of noncardiac complications following cardiac surgery. The hemodynamics of the aortic cannula has been proven to play a significant role in emboli generation and distribution. The aim of the current study was to perform a thorough numerical investigation in order to examine the effect of the design and orientation of the cannula used during cardiopulmonary bypass on the risk to develop cerebral embolism. Hemodynamic analyses compared numerical models of 27 cases consisting of six different cannula orientations, four aortic anatomies, and three cannula designs. The cannula designs included a straight-tip (ST) cannula, a moderately curved tip cannula (TIP1 ), and a sharp-angle curved cannula (TIP2 ). Outcome measures included hemodynamic parameters such as emanating jet velocity, jet velocity drop, maximal shear stress, aortic wall reaction, emboli pathlines and distribution between upper and lower vessels, and stagnation regions. Based on these parameters, the risks for hemolysis, atheroembolism, and cerebral embolism were evaluated and compared. On one hand, the jet emerging from the ST cannula generated large wall-shear stress at the aortic wall; this may have triggered the erosion and distribution of embolic atheromatous debris from the aortic arch. On the other hand, it diverted more emboli from the clamp region to the descending aorta and thus reduced the risk for cerebral embolism. The TIP1 cannula demonstrated less shear stress on the aortic wall and diverted more emboli from the clamp region toward the upper vessels. The TIP2 cannula exhibited a stronger emanating jet, higher shear stress inside the cannula, and highly disturbed flow, which was more stagnant near the clamp region. Current findings support the significant impact of the cannula design and orientation on emboli generation and distribution. Specifically, the straight tip cannula demonstrated a reduced risk of cerebral embolism, which may be pivotal in the clinical setting.

Fluid dynamics of coarctation of the aorta and effect of bicuspid aortic valve

Up to 80% of patients with coarctation of the aorta (COA) have a bicuspid aortic valve (BAV). Patients with COA and BAV have elevated risks of aortic complications despite successful surgical repair. The development of such complications involves the interplay between the mechanical forces applied on the artery and the biological processes occurring at the cellular level. The focus of this study is on hemodynamic modifications induced in the aorta in the presence of a COA and a BAV. For this purpose, numerical investigations and magnetic resonance imaging measurements were conducted with different configurations: (1) normal: normal aorta and normal aortic valve; (2) isolated COA: aorta with COA (75% reduction by area) and normal aortic valve; (3) complex COA: aorta with the same severity of COA (75% reduction by area) and BAV. The results show that the coexistence of COA and BAV significantly alters blood flow in the aorta with a significant increase in the maximal velocity, secondary flow, pressure loss, time-averaged wall shear stress and oscillatory shear index downstream of the COA. These findings can contribute to a better understanding of why patients with complex COA have adverse outcome even following a successful surgery.

Numerical simulation of Dean number and curvature effects on magneto-biofluid flow through a curved conduit

A numerical study is performed to investigate the magnetohydrodynamic viscous steady biofluid flow through a curved pipe with circular cross section under various conditions. A spectral method is applied as the principal tool for the numerical simulation with Fourier series, Chebyshev polynomials, collocation methods and an iteration method as secondary tools. The combined effects of Dean number, Dn, magnetic parameter, Mg, and tube curvature, δ, are studied. The flow patterns have been shown graphically for large Dean numbers as well as magnetic parameter and a wide range of curvatures, 0.01 ≤ δ≤ 0.2. Two-vortex solutions have been found. Axial velocity has been found to increase with an increase of Dean number, whereas it is suppressed with greater curvature and magnetic parameters. For high magnetic parameter and Dean number and low curvature, almost all the fluid vortex strengths are weak. The study is relevant to magnetohydrodynamic blood flow in the cardiovascular system.

A patient-specific study of type-B aortic dissection: evaluation of true-false lumen blood exchange

Background

Aortic dissection is a severe pathological condition in which blood penetrates between layers of the aortic wall and creates a duplicate channel – the false lumen. This considerable change on the aortic morphology alters hemodynamic features dramatically and, in the case of rupture, induces markedly high rates of morbidity and mortality.

Methods

In this study, we establish a patient-specific computational model and simulate the pulsatile blood flow within the dissected aorta. The k-ω SST turbulence model is employed to represent the flow and finite volume method is applied for numerical solutions. Our emphasis is on flow exchange between true and false lumen during the cardiac cycle and on quantifying the flow across specific passages. Loading distributions including pressure and wall shear stress have also been investigated and results of direct simulations are compared with solutions employing appropriate turbulence models.

Results

Our results indicate that (i) high velocities occur at the periphery of the entries; (ii) for the case studied, approximately 40% of the blood flow passes the false lumen during a heartbeat cycle; (iii) higher pressures are found at the outer wall of the dissection, which may induce further dilation of the pseudo-lumen; (iv) highest wall shear stresses occur around the entries, perhaps indicating the vulnerability of this region to further splitting; and (v) laminar simulations with adequately fine mesh resolutions, especially refined near the walls, can capture similar flow patterns to the (coarser mesh) turbulent results, although the absolute magnitudes computed are in general smaller.

Conclusions

The patient-specific model of aortic dissection provides detailed flow information of blood transport within the true and false lumen and quantifies the loading distributions over the aorta and dissection walls. This contributes to evaluating potential thrombotic behavior in the false lumen and is pivotal in guiding endovascular intervention. Moreover, as a computational study, mesh requirements to successfully evaluate the hemodynamic parameters have been proposed.

The effect of aortic wall and aortic leaflet stiffening on coronary hemodynamic: a fluid-structure interaction study

Pathologies of the aortic valve such as aortic sclerosis are thought to impact coronary blood flow. Recent clinical investigations have observed simultaneous structural and hemodynamic variations in the aortic valve and coronary arteries due to regional pathologies of the aortic valve. The goal of the present study is to elucidate this observed and yet unexplained phenomenon, in which a local pathology in the aortic valve region could potentially lead to the initiation or progression of coronary artery disease. Results revealed a considerable impact on the coronary flow, velocity profile, and consequently shear stress due to an increase in the aortic wall or aortic leaflet stiffness and thickness which concur with clinical observations. The cutoff value of 0.75 for fractional flow reserve was reached when the values of leaflet thickness and aortic wall stiffness were approximately twice and three times their normal value, respectively. Variations observed in coronary velocity profiles as well as wall shear stress suggest a possible link for the initiation of coronary artery disease.

A longitudinal study of Type-B aortic dissection and endovascular repair scenarios: computational analyses

Conservative medical treatment is commonly first recommended for patients with uncomplicated Type-B aortic dissection (AD). However, if dissection-related complications occur, endovascular repair or open surgery is performed. Here we establish computational models of AD based on radiological three-dimensional images of a patient at initial presentation and after 4-years of best medical treatment (BMT). Computational fluid dynamics analyses are performed to quantitatively investigate the hemodynamic features of AD. Entry and re-entries (functioning as entries and outlets) are identified in the initial and follow-up models, and obvious variations of the inter-luminal flow exchange are revealed. Computational studies indicate that the reduction of blood pressure in BMT patients lowers pressure and wall shear stress in the thoracic aorta in general, and flattens the pressure distribution on the outer wall of the dissection, potentially reducing the progressive enlargement of the false lumen. Finally, scenario studies of endovascular aortic repair are conducted. The results indicate that, for patients with multiple tears, stent-grafts occluding all re-entries would be required to effectively reduce inter-luminal blood communication and thus induce thrombosis in the false lumen. This implicates that computational flow analyses may identify entries and relevant re-entries between true and false lumen and potentially assist in stent-graft planning.

The concept of aortic replacement based on computational fluid dynamic analysis: patient-directed aortic replacement

OBJECTIVES

Aortic replacement is based on the aortic diameter in the absence of dissection or connective tissue diseases. Frequently, a number of different aortic-to-prosthetic anastomotic positions are possible depending on patient factors and surgeon preferences. High stress on residual aortic tissue may result in aneurysm formation or aneurysmal dilatation. Utilizing a computational fluid dynamic evaluation, we aimed to define possible optimal operative interventions with regard to the extent of aortic replacement.

METHODS

For proof of principle, a computational fluid dynamic (CFD) analysis, using Fluent 6.2 (Ansys UK Ltd, Sheffield, UK), was performed on a simplified ascending arch and descending aortic geometry. Wall shear stress in three dimensions was assessed for the standard operations: ascending aortic replacement, arch replacement and proximal descending aortic replacement.

RESULTS

Hermiarch replacement is superior to isolated ascending aortic replacement with regard to residual stress analysis on tissues (up to a 10-fold reduction). Aortic arch replacement with island implantation of the supra-aortic vessels may potentially result in high stress on the residual aorta (10-fold increase). Aortic arch replacement with individual supra-aortic vessel implantation may result in areas of high stress (10-fold increase) on native vessels if an inadequate length of supra-aortic tissue is not resected, regardless of it being aneurysmal.

CONCLUSIONS

Computational fluid dynamic evaluation, which will have to be patient-specific, 3D anatomical and physiological, potentially has enormous implications for operative strategy in aortic replacement surgery. CFD analysis may direct the replacement of normal-diameter aortas in the future.

Hemodynamic variations due to spiral blood flow through four patient-specific bifurcated stent graft configurations for the treatment of abdominal aortic aneurysms

Endovascular repair is now a recognised procedure for treating abdominal aortic aneurysms. However, post-operative complications such as stent graft migration and thrombus may still occur. To assess these complications numerically, the correct input boundary conditions, which include the full human aorta with associated branching, should be included. Four patient-specific computed tomography scanned bifurcated stent grafts (SGs) were modelled and attached onto a full human aorta, which included the ascending, aortic arch and descending aortas. Two of the SG geometries had a twisted leg configuration, while the other two had conventional nontwisted leg configurations. Computational fluid dynamics was completed for both geometries and the hemodynamics assessed. The complexity of the flow patterns and secondary flows were influenced by the inclusion of the full human aorta at the SG proximal section. During the decelerating phase significant recirculations occurred along the main body of all SG configurations. The inclusion of the full human aorta did not impact the velocity contours within the distal legs and there was no difference in drag forces with the SG containing the full human aorta and those without. A twisted leg configuration further promoted a spiral flow formation along its distal legs.

Stent graft performance in the treatment of abdominal aortic aneurysms: the influence of compliance and geometry

The long-term success of the endovascular procedure for the treatment of Abdominal Aortic Aneurysms (AAAs ) depends on the secure fixation of the proximal end and the geometry of the stent-graft (SG) device. Variations in SG types can affect proximal fixation and SG hemodynamics. Such hemodynamic variations can have a catastrophic effect on the vascular system and may result from a SG/arterial wall compliance mismatch and the sudden decrease in cross-sectional area at the bifurcation, which may result in decreased distal perfusion, increased pressure wave reflection and increased stress at the interface between the stented and non-stented portion of the vessel. To examine this compliance mismatch, a commercial SG device was tested experimentally under a physiological pressure condition in a silicone AAA model based on computed tomography scans. There was a considerable reduction in compliance of 54% and an increase in the pulse wave velocity of 21%, with a significant amount of the forward pressure wave being reflected. To examine the SG geometrical effects, a commercial bifurcated geometry was compared computationally and experimentally with a geometrical taper in the form of a blended section, which provided a smooth transition from the proximal end to both iliac legs. The sudden contraction of commercial SG at the bifurcation region causes flow separation within the iliac legs, which is known to cause SG occlusion and increased proximal pressure. The blended section along the bifurcation region promotes a greater uniformity of the fluid flow field within the distal legs, especially, during the deceleration phase with reduced boundary layer reversal. In order to reduce the foregoing losses, abrupt changes of cross-section should be avoided. Geometrical tapers could lead to improved clinical outcomes for AAA SGs.

Embolic events caused by aortic thrombi: an underestimated entity?

Stroke and other thromboembolic events are mainly caused by emboli from heart, aorta and other arteries. In this paper we describe a group of 5 middle-aged patients suffering from emboli caused by large thrombi in the aorta. Since the development of giant thrombi under high flow conditions in the aorta is a pathophysiological process which is not well understood, a model of flow distribution by numerically simulating the Navier-Stokes equation for an incompressible fluid was generated. This model simulated how such thrombi may develop in the aorta. We hypothesize that large thrombi issuing from the aortic vessel wall represent a underestimated entity in middleaged persons and are probably overlooked as the cause of stroke or other embolic events in some cases.

A computational fluid dynamics study on geometrical influence of the aorta on haemodynamics

OBJECTIVES

Cardiovascular diseases, such as atherosclerosis and aneurysm, are closely associated with haemodynamic factors that are governed by luminal geometry. The present work aimed to study the effect of geometrical variation of aging aortas on haemodynamics.

METHODS

Six aged subjects with intricate geometrical features, such as bulging or twisted supra-aortic arteries, sharply curved arch and double-curved descending aorta, were chosen from our medical database. These six geometrically variant aortas were reconstructed and the pulsatile nature of the blood flow of these subject-specific aorta models investigated using computational fluid dynamics simulations. Realistic time-dependent boundary conditions are prescribed for various arteries of the aorta models.

RESULTS

This study suggests that haemodynamics in the human aorta is highly dependent on geometrical features. The positioning and contouring of the supra-aortic arteries may be associated with the skewness of velocity profiles. The flow profiles in the aortic arch or bends are generally skewed towards the inner curvature wall and this skewness may give rise to the formation of secondary flow in the inner curvature wall of the distal arch. The degree of vorticity in the distal aortic arch is found to be related to the arch curvature. The helical nature of aortic haemodynamics is predominant in the systole phrase when it begins with a left-handed rotation and then vanishes in the ascending aorta, whereas a right-handed rotation persists in the distal aortic arch. Lower wall shear stress is also found in the ascending regions where secondary flow is present.

CONCLUSIONS

The aorta with an irregular contour and large degree of curvature at its arch favours the development of the intra-aortic secondary flow that subsequently relates to the pathogenesis of atheroma. The present study identifies the general trend of haemodynamic behaviours associated with various local geometrical features. Combining the knowledge of the correlation between haemodynamics and the underlying risks in the development of cardiovascular diseases, our study hopes to provide a better understanding of the relationship between aortic morphology and developing pathobiology of cardiovascular diseases. As such, early medical planning as well as surgical interventions can be designed to retard or prevent the development of cardiovascular diseases.