Numerical simulation of non-Newtonian blood flow dynamics in human thoracic 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 and emboli transport patterns during venoarterial extracorporeal membrane oxygenation: A computational fluid dynamics study

PROBLEM

Despite advances in Venoarterial Extracorporeal Membrane Oxygenation (VA-ECMO), a significant mortality rate persists due to complications. The non-physiological blood flow dynamics of VA-ECMO may lead to neurological complications and organ ischemia. Continuous retrograde high-flow oxygenated blood enters through a return cannula placed in the femoral artery which opposes the pulsatile deoxygenated blood ejected by the left ventricle (LV), which impacts upper body oxygenation and subsequent hyperoxemia. The complications underscore the critical need to comprehend the impact of VA-ECMO support level and return cannula size, as mortality remains a significant concern.

AIM

The aim of this study is to predict and provide insights into the complications associated with VA-ECMO using computational fluid dynamics (CFD) simulations. These complications will be assessed by characterising blood flow and emboli transport patterns through a comprehensive analysis of the influence of VA-ECMO support levels and arterial return cannula sizes.

METHODS

Patient-specific 3D aortic and major branch models, derived from a male patient's CT scan during VA-ECMO undergoing respiratory dysfunction, were analyzed using CFD. The investigation employed species transport and discrete particle tracking models to study ECMO blood (oxygenated) mixing with LV blood (deoxygenated) and to trace emboli transport patterns from potential sources (circuit, LV, and aorta wall). Two cannula sizes (15 Fr and 19 Fr) were tested alongside varying ECMO pump flow rates (50%, 70%, and 90% of the total cardiac output).

RESULTS

Cannula size did not significantly affect oxygen transport. At 90% VA-ECMO support, all arteries distal to the aortic arch achieved 100% oxygen saturation. As support level decreased, oxygen transport to the upper body also decreased to a minimum saturation of 73%. Emboli transport varied substantially between emboli origin and VAECMO support level, with the highest risk of cerebral emboli coming from the LV with a 15 Fr cannula at 90% support.

CONCLUSION

Arterial return cannula sizing minimally impacted blood oxygen distribution; however, it did influence the distribution of emboli released from the circuit and aortic wall. Notably, it was the support level alone that significantly affected the mixing zone of VA-ECMO and cardiac blood, subsequently influencing the risk of embolization of the cardiogenic source and oxygenation levels across various arterial branches.

Design of experiment approach to identify the dominant geometrical feature of left coronary artery influencing atherosclerosis

Background and Objective. Coronary artery geometry heavily influences local hemodynamics, potentially leading to atherosclerosis. Consequently, the unique geometrical configuration of an individual by birth can be associated with future risk of atherosclerosis. Although current researches focus on exploring the relationship between local hemodynamics and coronary artery geometry, this study aims to identify the order of influence of the geometrical features through systematic experiments, which can reveal the dominant geometrical feature for future risk assessment.Methods. According to Taguchi's method of design of experiment (DoE), the left main stem (LMS) length (lLMS), curvature (kLMS), diameter (dLMS) and the bifurcation angle between left anterior descending (LAD) and left circumflex (LCx) artery (αLAD-LCx) of two reconstructed patient-specific left coronary arteries (LCA) were varied in three levels to create L9 orthogonal array. Computational fluid dynamic (CFD) simulations with physiological boundary conditions were performed on the resulting eighteen LCA models. Average helicity intensity (h2) and relative atheroprone area (RAA) of near-wall hemodynamic descriptors were analyzed.Results. The proximal LAD (LADproximal) was identified to be the most atheroprone region of the left coronary artery due to higherh2,large RAA of time averaged wall shear stress (TAWSS < 0.4 Pa), oscillatory shear index (OSI ∼ 0.5) and relative residence time (RRT > 4.17 Pa-1). In both patient-specific cases, based onh2and TAWSS,dlmsis the dominant geometric parameter while based on OSI and RRT,αLAD-LCxis the dominant one influencing hemodynamic condition in proximal LAD (p< 0.05). Based on RRT, the rank of the geometrical factors is:αLAD-LCx>dLMS>lLMS>kLMS, indicating thatαLAD-LCxis the most dominant geometrical factor affecting hemodynamics at proximal LAD which may influence atherosclerosis.Conclusion. The proposed identification of the rank of geometrical features of LCA and the dominant feature may assist clinicians in predicting the possibility of atherosclerosis, of an individual, long before it will occur. This study can further be translated to be used to rank the influence of several arterial geometrical features at different arterial locations to explore detailed relationships between the arterial geometrical features and local hemodynamics.

Effect of rheological models on pulsatile hemodynamics in a multiply afflicted descending human aortic network

In the cardiovascular diseased (CVD) conditions, it is essential to choose a suitable rheological model for capturing the correct physics behind the hemodynamic in the multiply afflicted diseased arterial network. This study investigates the effect of blood rheology on hemodynamics in a blood vessel with abdominal aortic aneurysm (AAA) and right internal iliac stenosis (RIIAS). A model with AAA and RIIAS is reconstructed from a human subject's computed tomography (CT) data. Localized mesh generation and pulsatile inflow condition are considered. Non-Newtonian models such as the Power-law, Carreau, Cross, and Herschel Berkley models are used in simulations. The outcome from a validated computational model is compared with the Newtonian model to identify the suitable model for dealing with pathological complications under consideration. The capabilities and significance of various rheological models are also examined via Wall Pressure (WP), Wall Shear Stress (WSS), velocity, Global non-Newtonian importance factor (IG), Vorticity Streamlines, and Swirling Strength. It is noted that during the entire cardiac cycle, the IG factor of the cross model is found to be relatively more significant. Power Law depicts larger IG factor during peak systole and early diastole. Also, the cross model depicts larger WSS, WPS, swirling strength distribution and vorticity during the peak systolic and diastolic phases It is noted that IG ∼0.02 is an appropriate non-Newtonian blood activity cut-off value in the descending abdominal artery having AAA and RIIAS. The critical important WSS values are in the range of 0-9 Pa which is stated in WSS contour plot.

Statistical shape representation of the thoracic aorta: accounting for major branches of the aortic arch

Statistical shape modeling (SSM) is an emerging tool for risk assessment of thoracic aortic aneurysm. However, the head branches of the aortic arch are often excluded in SSM. We introduced an SSM strategy based on principal component analysis that accounts for aortic branches and applied it to a set of patient scans. Computational fluid dynamics were performed on the reconstructed geometries to identify the extent to which branch model accuracy affects the calculated wall shear stress (WSS) and pressure. Surface-averaged and location-specific values of pressure did not change significantly, but local WSS error was high near branches when inaccurately modeled.

Impact of wall displacements on the large-scale flow coherence in ascending aorta

In the context of aortic hemodynamics, uncertainties affecting blood flow simulations hamper their translational potential as supportive technology in clinics. Computational fluid dynamics (CFD) simulations under rigid-walls assumption are largely adopted, even though the aorta contributes markedly to the systemic compliance and is characterized by a complex motion. To account for personalized wall displacements in aortic hemodynamics simulations, the moving-boundary method (MBM) has been recently proposed as a computationally convenient strategy, although its implementation requires dynamic imaging acquisitions not always available in clinics. In this study we aim to clarify the real need for introducing aortic wall displacements in CFD simulations to accurately capture the large-scale flow structures in the healthy human ascending aorta (AAo). To do that, the impact of wall displacements is analyzed using subject-specific models where two CFD simulations are performed imposing (1) rigid walls, and (2) personalized wall displacements adopting a MBM, integrating dynamic CT imaging and a mesh morphing technique based on radial basis functions. The impact of wall displacements on AAo hemodynamics is analyzed in terms of large-scale flow patterns of physiological significance, namely axial blood flow coherence (quantified applying the Complex Networks theory), secondary flows, helical flow and wall shear stress (WSS). From the comparison with rigid-wall simulations, it emerges that wall displacements have a minor impact on the AAo large-scale axial flow, but they can affect secondary flows and WSS directional changes. Overall, helical flow topology is moderately affected by aortic wall displacements, whereas helicity intensity remains almost unchanged. We conclude that CFD simulations with rigid-wall assumption can be a valid approach to study large-scale aortic flows of physiological significance.

Assessing the effect of manufacturing defects and non-Newtonian blood model on flow behaviors of additively manufactured Gyroid TPMS structures

In the field of medical engineering, Triply Periodic Minimal Surfaces (TPMS) structures have been studied widely owing to their physical attributes similar to those of human bones. Computational Fluid Dynamics (CFD) is often used to reveal the interaction between structural architectures and flow fields. Nevertheless, a comprehensive study on the effect of manufacturing defects and non-Newtonian behavior on the fluid responses in TPMS scaffolds is still lacking. Therefore, the present study fabricated Gyroid TPMS with four relative densities from 0.1 to 0.4. Non-destructive techniques were used to examine surface roughness and geometric deviation. We found that the manufacturing defects had a minor effect on fluid responses. The pressure drop comparison between defect-containing and defect-free models could be differed up to 7%. The same comparison for the average shear stress showed a difference up to 23%, in which greater deviation between both models was observed at higher relative density. On the contrary, the viscosity model played a significant role in flow prediction. By comparing the Newtonian model with Carreau-Yasuda non-Newtonian model, the resulting pressure drop and average wall shear stress from non-Newtonian viscosity could be higher than those of the Newtonian model by more than a factor of two. In addition, we matched the fluid-induced shear stress from both viscosity models with desirable ranges of shear stresses for tissue growth obtained from the literature. Up to 70% from the Newtonian model fell within the desirable range while the matching stress reduced to lower than 8% for the non-Newtonian results. Furthermore, by correlating geometric features with physical outputs, the geometric deviation was seen associated with surface curvature while the local shear stress revealed a strong correlation with inclination angle. Overall, the present work emphasized the importance of the viscosity model for CFD analysis of the scaffolds, especially when resulting fluid-induced wall shear stress is of interest. In addition, the geometric correlation has introduced the alternative consideration of structural architectures from local perspectives, which could assist the further comparison and optimization among different porous scaffolds in the future.

The effect of arterial cannula tip position on differential hypoxemia during venoarterial extracorporeal membrane oxygenation

Interaction between native ventricular output and venoarterial extracorporeal membrane oxygenation (VA ECMO) flow may hinder oxygenated blood flow to the aortic arch branches, resulting in differential hypoxemia. Typically, the arterial cannula tip is placed in the iliac artery or abdominal aorta. However, the hemodynamics of a more proximal arterial cannula tip have not been studied before. This study investigated the effect of arterial cannula tip position on VA ECMO blood flow to the upper extremities using computational fluid dynamics simulations. Four arterial cannula tip positions (P1. common iliac, P2. abdominal aorta, P3. descending aorta and P4. aortic arch) were compared with different degrees of cardiac dysfunction and VA ECMO support (50%, 80% and 90% support). P4 was able to supply oxygenated blood to the arch vessels at all support levels, while P1 to P3 only supplied the arch vessels during the highest level (90%) of VA ECMO support. Even during the highest level of support, P1 to P3 could only provide oxygenated VA-ECMO flow at 0.11 L/min to the brachiocephalic artery, compared with 0.5 L/min at P4. This study suggests that cerebral perfusion of VA ECMO flow can be increased by advancing the arterial cannula tip towards the aortic arch.

Patient-Specific Image-Based Computational Fluid Dynamics Analysis of Abdominal Aorta and Branches

The complicated abdominal aorta and its branches are a portion of the circulatory system prone to developing atherosclerotic plaque and aneurysms. These disorders are closely connected to the changing blood flow environment that the area’s complicated architecture produces (between celiac artery and iliac artery bifurcation); this phenomenon is widespread at arterial bifurcations. Based on computed tomography angiography (CTA) scans, this current work offers a numerical analysis of a patient-specific reconstruction of the abdominal aorta and its branches to identify and emphasize the most likely areas to develop atherosclerosis. The simulations were run following the heart cycle and under physiological settings. The wall shear stress (WSS), velocity field, and streamlines were examined. According to the findings, complex flow is primarily present at the location of arterial bifurcations, where abnormal flow patterns create recirculation zones with low and fluctuating WSS (<0.5 Pa), which are known to affect endothelial homeostasis and cause adverse vessel remodeling. The study provides a patient-specific hemodynamic analysis model, which couples in vivo CT imaging with in silico simulation under physiological circumstances. The study offers quantitative data on the range fluctuations of important hemodynamic parameters, such as WSS and recirculation region expansion, which are directly linked to the onset and progression of atherosclerosis. The findings could also help drug targeting at this vascular level by understanding blood flow patterns in the abdominal aorta and its branches.

Numerical investigation of different viscosity models on pulsatile blood flow of thoracic aortic aneurysm (TAA) in a patient-specific model

Aortic aneurysm is one of the most common aortic diseases that can lead to unfortunate consequences. Numerical simulations have an important role in the prediction of the aftereffects of vascular diseases including aneurysm. In this research, numerical simulation of pulsatile blood flow is performed for a 3-dimensional patient-specific model of a thoracic aortic aneurysm (TAA). Since the choice of blood viscosity model may have a significant impact on the simulation results, the effects of four non-Newtonian models of blood viscosity namely Carreau, Casson, Herschel-Bulkley, power low, and the Newtonian model on the wall shear stress (WSS) distribution, shear rate, and oscillatory shear index (OSI) have been analyzed. Simulation results showed that all the non-Newtonian and Newtonian models generally, predict similar patterns for blood flow and shear rate. At high flow rates in the cardiac cycle, the WSS value for all the models are similar to each other except for the power-law model due to the shear thinning behavior. All models predict high values of OSI on the inner wall of the ascending aorta and broad areas of the inner wall of the aneurysm sac. However, the Newtonian model predicts the OSI less than the non-Newtonian models in some areas of the aneurysm sac. Results indicated that the Newtonian model generally can predict the hemodynamic parameters of the blood flow similar to the non-Newtonian but for more precise analysis and to predict the regions prone to rupture and atherosclerosis, choosing a proper non-Newtonian model is recommended.

Pulsatile flow of thixotropic blood in artery under external body acceleration

This work presents a numerical technique for simulating non-Newtonian blood flow in human's arteries driven by an oscillating pressure gradient. The blood is considered as a thixotropic fluid and its structural properties are considered to obey Moore's thixotropic model as a constitutive equation. The equations of motion are simplified considering the flow laminar, axisymmetric and the fluid incompressible. A numerical solution is presented using finite difference method in order to compute the velocity field and wall shear stress distribution. The numerical results obtained have been validated with the analytical solution available in the literature. Furthermore, the effect of the structural properties, the average of the pressure gradient and the external acceleration on the velocity and wall shear stress distribution is investigated. These results reveal the influence of the different parameters studied on the pipe flow response of the thixotropic fluid.

Simulation of Mechanical Heart Valve Dysfunction and the Non-Newtonian Blood Model Approach

The mechanical heart valve (MHV) is commonly used for the treatment of cardiovascular diseases. Nonphysiological hemodynamic in the MHV may cause hemolysis, platelet activation, and an increased risk of thromboembolism. Thromboembolism may cause severe complications and valve dysfunction. This paper thoroughly reviewed the simulation of physical quantities (velocity distribution, vortex formation, and shear stress) in healthy and dysfunctional MHV and reviewed the non-Newtonian blood flow characteristics in MHV. In the MHV numerical study, the dysfunction will affect the simulation results, increase the pressure gradient and shear stress, and change the blood flow patterns, increasing the risks of hemolysis and platelet activation. The blood flow passes downstream and has obvious recirculation and stagnation region with the increased dysfunction severity. Due to the complex structure of the MHV, the non-Newtonian shear-thinning viscosity blood characteristics become apparent in MHV simulations. The comparative study between Newtonian and non-Newtonian always shows the difference. The shear-thinning blood viscosity model is the basics to build the blood, also the blood exhibiting viscoelastic properties. More details are needed to establish a complete and more realistic simulation.

A Detailed Study to Discover the Trade between Left Atrial Blood Flow, Expression of Calcium-Activated Potassium Channels and Valvular Atrial Fibrillation

Background: The present study aimed to explore the correlation between calcium-activated potassium channels, left atrial flow field mechanics, valvular atrial fibrillation (VAF), and thrombosis. The process of transforming mechanical signals into biological signals has been revealed, which offers new insights into the study of VAF. Methods: Computational fluid dynamics simulations use numeric analysis and algorithms to compute flow parameters, including turbulent shear stress (TSS) and wall pressure in the left atrium (LA). Real-time PCR and western blotting were used to detect the mRNA and protein expression of IKCa2.3/3.1, ATK1, and P300 in the left atrial tissue of 90 patients. Results: In the valvular disease group, the TSS and wall ressure in the LA increased, the wall pressure increased in turn in all disease groups, mainly near the mitral valve and the posterior portion of the LA, the increase in TSS was the most significant in each group near the mitral valve, and the middle and lower part of the back of the LA and the mRNA expression and protein expression levels of IKCa2.3/3.1, AKT1, and P300 increased (p < 0.05) (n = 15). The present study was preliminarily conducted to elucidate whether there might be a certain correlation between IKCa2.3 and LA hemodynamic changes. Conclusions: The TSS and wall pressure changes in the LA are correlated with the upregulation of mRNA and protein expression of IKCa2.3/3.1, AKT1, and P300.

Steady-state solutions of one-dimensional equations of non-Newtonian hemodynamics

This paper is devoted to obtaining and analysis of steady-state solutions of one-dimensional equations for the simulation of blood flow when the non-Newtonian nature of blood is taken into account. The models, based on the rheological relations, widely used for the blood, are considered. The expressions for the nonlinear frictional term are presented. For the Power Law, Simplified Cross, and Quemada models, the exact integrals of the nonlinear ordinary differential equation, obtained from the averaged momentum equation, are obtained. It is demonstrated that several solutions exist for every rheological model, but the physically relevant solutions can be selected by the appropriate value of Mach number. The effects of the velocity profile and the value of hematocrit on the steady-state solutions are analyzed. It is demonstrated that the flattening of the velocity profile, which is typical for the blood, leads to the diminishing of the length of the interval, where the solution exists. The same effect is observed when the hematocrit value is increased.

Tunable Blood Shunt for Neonates With Complex Congenital Heart Defects

Despite advancements in procedures and patient care, mortality rates for neonatal recipients of the Norwood procedure, a palliation for single ventricle congenital malformations, remain high due to the use of a fixed-diameter blood shunt. In this study, a new geometrically tunable blood shunt was investigated to address limitations of the current treatment paradigm (e.g., Modified Blalock-Taussig Shunt) by allowing for controlled modulation of blood flow through the shunt to accommodate physiological changes due to the patient’s growth. First, mathematical and computational cardiovascular models were established to investigate the hemodynamic requirements of growing neonatal patients with shunts and to inform design criteria for shunt diameter changes. Then, two stages of prototyping were performed to design, build and test responsive hydrogel systems that facilitate tuning of the shunt diameter by adjusting the hydrogel’s degree of crosslinking. We examined two mechanisms to drive crosslinking: infusion of chemical crosslinking agents and near-UV photoinitiation. The growth model showed that 15–18% increases in shunt diameter were required to accommodate growing patients’ increasing blood flow; similarly, the computational models demonstrated that blood flow magnitudes were in agreement with previous reports. These target levels of diameter increases were achieved experimentally with model hydrogel systems. We also verified that the photocrosslinkable hydrogel, composed of methacrylated dextran, was contact-nonhemolytic. These results demonstrate proof-of-concept feasibility and reflect the first steps in the development of this novel blood shunt. A tunable shunt design offers a new methodology to rebalance blood flow in this vulnerable patient population during growth and development.

Development of Inspired Therapeutics Pediatric VAD: Computational Analysis and Characterization of VAD V3

PURPOSE

Pediatric heart failure patients remain in critical need of a dedicated mechanical circulatory support (MCS) solution as development efforts for specific pediatric devices continue to fall behind those for the adult population. The Inspired Pediatric VAD is being developed as a pediatric specific MCS solution to provide up to 30-days of circulatory or respiratory support in a compact modular package that could allow for patient ambulation during treatment.

METHODS

Hydrodynamic performance (flows, pressures), impeller/rotor mechanical properties (torques, forces), and flow shear stress and residence time distributions of the latest design version, Inspired Pediatric VAD V3, were numerically predicted and investigated using computational fluid dynamics (CFD) software (SolidWorks Flow Simulator).

RESULTS

Hydrodynamic performance was numerically predicted, indicating no change in flow and pressure head compared to the previous device design (V2), while displaying increased impeller/rotor torques and translation forces enabled by improved geometry. Shear stress and flow residence time volumetric distributions are presented over a range of pump rotational speeds and flow rates. At the lowest pump operating point (3000 RPM, 0.50 L/min, 75 mmHg), 79% of the pump volume was in the shear stress range of 0-10 Pa with < 1% of the volume in the critical range of 150-1000 Pa for blood damage. At higher speed and flow (5000 RPM, 3.50 L/min, 176 mmHg), 65% of the volume resided in the 0-10 Pa range compared to 2.3% at 150-1000 Pa.

CONCLUSIONS

The initial computational characterization of the Inspired Pediatric VAD V3 is encouraging and future work will include device prototype testing in a mock circulatory loop and acute large animal model.

In vitro flow study in a compliant abdominal aorta phantom with a non-Newtonian blood-mimicking fluid

In vitro aortic flow simulators allow studying hemodynamics with a wider range of flow visualization techniques compared to in vivo medical imaging and without the limitations of invasive examinations. This work aims to develop an experimental bench to emulate the pulsatile circulation in a realistic aortic phantom. To mimic the blood shear thinning behavior, a non-Newtonian aqueous solution is prepared with glycerin and xanthan gum polymer. The flow is compared to a reference flow of Newtonian fluid. Particle image velocimetry is carried out to visualize 2D velocity fields in a phantom section. The experimental loop accurately recreates flowrates and pressure conditions and preserves the shear-thinning properties of the non-Newtonian fluid. Velocity profiles, shear rate, and shear stress distribution maps show that the Newtonian fluid tends to dampen the observed velocities. Preferential asymmetrical flow paths are observed in a diameter narrowing region and amplified in the non-Newtonian case. Wall shear stresses are about twice higher in the non-Newtonian case. This study shows new insights on flow patterns, velocity and shear stress distributions compared to rigid and simplified geometry aorta phantom with Newtonian fluid flows studies. The use of a non-Newtonian blood analog shows clear differences in flows compared to the Newtonian one in this compliant patient-specific geometry. The development of this aortic simulator is a promising tool to better analyze and understand aortic hemodynamics and to aid in clinical decision-making.

Comparison of Non-Newtonian Models of One-Dimensional Hemodynamics

The paper is devoted to the comparison of different one-dimensional models of blood flow. In such models, the non-Newtonian property of blood is considered. It is demonstrated that for the large arteries, the small parameter is observed in the models, and the perturbation method can be used for the analytical solution. In the paper, the simplified nonlinear problem for the semi-infinite vessel with constant properties is solved analytically, and the solutions for different models are compared. The effects of the flattening of the velocity profile and hematocrit value on the deviation from the Newtonian model are investigated.

Simulation of non-Newtonian flow of blood in a modified laparoscopic forceps used in minimally invasive surgery

During surgeries, blood often oozes out of the operated tissue and this has to be sucked out by the S-I device. Blood is more viscous than saline, the cleaning fluid is used in the S-I process. Therefore, for a more comprehensive CFD flow analysis of the improved forceps is simulated in the present work for different driving pressures. The resulting flow rate of blood is compared among the prospective designs and the S-I device currently in use. The new surgical forceps eliminates re-insertion of dissector with suction-irrigator and is reusable, multi-functional, non-toxic, corrosion resistant, toughened, and cost-effective.

Transient numerical simulation of the right coronary artery originating from the left sinus and the effect of its acute take-off angle on hemodynamics

Background

An anomalous origin of the right coronary artery from the left coronary artery sinus is usually characterized by an acute take-off angle. Most affected patients have no clinical symptoms; however, some patients have decreased blood flow into the right coronary artery during exercise, which can lead to symptoms such as myocardial ischemia. Most researchers who have studied an anomalous origin of the right coronary artery from the left coronary artery sinus have done so through clinical cases. In this study, we used numerical simulation to evaluate the hemodynamics of this condition and the effect of an acute take-off angle on hemodynamic parameters. We expect that the results of this study will help in further understanding the clinical symptoms of this anomaly and the hemodynamic impact of an acute take-off angle.

Methods

Three-dimensional models were reconstructed based on the computed tomography images from 16 patients with a normal right coronary artery and 26 patients with an anomalous origin of the right coronary artery from the left coronary artery sinus. A numerical simulation of a two-way fluid-structure interaction was executed with ANSYS Workbench software. The blood was assumed to be an incompressible Newtonian fluid, and the vessel was assumed to be an isotropic, linear elastic material. Hemodynamic parameters and the effect of an acute take-off angle were statistically analyzed.

Results

During the systolic period, the wall pressure in the right coronary artery was significantly reduced in patients with an anomalous origin of the right coronary artery (t =1.32 s, P=0.0001; t =1.34-1.46 s, P<0.0001). The wall shear stress in the abnormal group was higher at the beginning of the systolic period (t =1.24 s, P=0.0473; t =1.26 s, P=0.0193; t =1.28 s, P=0.0441). The acute take-off angle was smaller in patients with clinical symptoms (27.81°±4.406°) than in patients without clinical symptoms (31.86°±2.789°; P=0.017). In the symptomatic group, pressure was negatively correlated with the acute take-off angle (P=0.0185-0.0341, r=-0.459 to -0.4167).

Conclusions

This study shows that an anomalous origin of the right coronary artery from the left coronary artery sinus causes changes in hemodynamic parameters, and that an acute take-off angle in patients with this anomaly is associated with terminal ischemia of the right coronary artery.

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.

Patient-Specific Numerical Analysis of Coronary Flow in Children With Intramural Anomalous Aortic Origin of Coronary Arteries

Unroofing surgery for anomalous aortic origin of a coronary artery (AAOCA) alters coronary anatomy by opening the intramural segment so that the anomalous coronary orifice arises perpendicularly from appropriate aortic sinus. Computational fluid dynamics modeling (CFD) allows for quantification of hemodynamics linked to morbidity such as wall shear stress (WSS), relative to patient-specific features like the angle of origin (AO). We hypothesize that CFD will reveal abnormal WSS indices in unroofed arteries that are related to AO. Six AAOCA patients (3 left, 3 right) status post unroofing (median = 13.5 years, range 9-17) underwent cardiac magnetic resonance imaging. CFD models were created from pre (n = 2) and postunroofing (n = 6) cardiac magnetic resonance imaging data, for the anomalous and contralateral normally-arising arteries. Downstream vasculature was represented by lumped parameter networks. Time-averaged WSS (TAWSS) and oscillatory shear index (OSI) were quantified relative to AO and measured hemodynamics. TAWSS was elevated along the outer wall of the normally-arising left vs right coronary arteries, as well as along unroofed left vs right coronary arteries (n = 6/group). No significant differences were noted when comparing unroofed and same-sided normally-arising coronaries. TAWSS was reduced after unroofing (eg, 276 ± 28 dyne/cm² vs 91 ± 15 dyne/cm²; n = 2/group). Models with more acute preoperative AO indicated lower TAWSS at the proximity of ostium. Differences in OSI were not significant. Different flow patterns exist natively between right and left coronary arteries. Unroofing may normalize TAWSS but with variance related to the AO. This study suggests CFD may help stratify risk in AAOCA.

Assessment of surface roughness and blood rheology on local coronary haemodynamics: a multi-scale computational fluid dynamics study

The surface roughness of the coronary artery is associated with the onset of atherosclerosis. The study applies, for the first time, the micro-scale variation of the artery surface to a 3D coronary model, investigating the impact on haemodynamic parameters which are indicators for atherosclerosis. The surface roughness of porcine coronary arteries have been detailed based on optical microscopy and implemented into a cylindrical section of coronary artery. Several approaches to rheology are compared to determine the benefits/limitations of both single and multiphase models for multi-scale geometry. Haemodynamic parameters averaged over the rough/smooth sections are similar; however, the rough surface experiences a much wider range, with maximum wall shear stress greater than 6 Pa compared to the approximately 3 Pa on the smooth segment. This suggests the smooth-walled assumption may neglect important near-wall haemodynamics. While rheological models lack sufficient definition to truly encompass the micro-scale effects occurring over the rough surface, single-phase models (Newtonian and non-Newtonian) provide numerically stable and comparable results to other coronary simulations. Multiphase models allow for phase interactions between plasma and red blood cells which is more suited to such multi-scale models. These models require additional physical laws to govern advection/aggregation of particulates in the near-wall region.

Non-Newtonian Effects on Patient-Specific Modeling of Fontan Hemodynamics

The Fontan procedure is a common palliative surgery for congenital single ventricle patients. In silico and in vitro patient-specific modeling approaches are widely utilized to investigate potential improvements of Fontan hemodynamics that are related to long-term complications. However, there is a lack of consensus regarding the use of non-Newtonian rheology, warranting a systematic investigation. This study conducted in silico patient-specific modeling for twelve Fontan patients, using a Newtonian and a non-Newtonian model for each patient. Differences were quantified by examining clinically relevant metrics: indexed power loss (iPL), indexed viscous dissipation rate (iVDR), hepatic flow distribution (HFD), and regions of low wall shear stress (A_WSS). Four sets of "non-Newtonian importance factors" were calculated to explore their effectiveness in identifying the non-Newtonian effect. No statistical differences were observed in iPL, iVDR, and HFD between the two models at the population-level, but large inter-patient variations exist. Significant differences were detected regarding A_WSS, and its correlations with non-Newtonian importance factors were discussed. Additionally, simulations using the non-Newtonian model were computationally faster than those using the Newtonian model. These findings distinguish good importance factors for identifying non-Newtonian rheology and encourage the use of a non-Newtonian model to assess Fontan hemodynamics.

Computational Fluid Dynamics Modelling and Simulation of an Inclined Horizontal Axis Hydrokinetic Turbine

In this contribution, unsteady three-dimensional numerical simulations of the water flow through a horizontal axis hydrokinetic turbine (HAHT) of the Garman type are performed. This study was conducted in order to estimate the influence of turbine inclination with respect to the incoming flow on turbine performance and forces acting on the rotor, which is studied using a time-accurate Reynolds-averaged Navier-Stokes (RANS) commercial solver. Changes of the flow in time are described by a physical transient model based on two domains, one rotating and the other stationary, combined with a sliding mesh technique. Flow turbulence is described by the well-established Shear Stress Transport (SST) model using its standard and transitional versions. Three inclined operation conditions have been analyzed for the turbine regarding the main stream: 0° (SP configuration, shaft parallel to incoming velocity), 15° (SI15 configuration), and 30° (SI30 configuration). It was found that the hydrodynamic efficiency of the turbine decreases with increasing inclination angles. Besides, it was obtained that in the inclined configurations, the thrust and drag forces acting on rotor were lower than in the SP configuration, although in the former cases, blades experience alternating loads that may induce failure due to fatigue in the long term. Moreover, if the boundary layer transitional effects are included in the computations, a slight increase in the power coefficient is computed for all inclination configurations.

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.

Computational fluid dynamics in abdominal aorta bifurcation: non-Newtonian versus Newtonian blood flow in a real case study

Hemodynamic in abdominal aorta bifurcation was investigated in a real case using computational fluid dynamics. A Newtonian and non-Newtonian (Walburn-Schneck) viscosity models were compared. The geometrical model was obtained by 3D reconstruction from CT-scan and hemodynamic parameters obtained by laser-Doppler. Blood was assumed incompressible fluid, laminar flow in transient regime and rigid vessel wall. Finite volume-based was used to study the velocity, pressure, wall shear stress (WSS) and viscosity throughout cardiac cycle. Results obtained with Walburn-Schneck's model, during systole, present lower viscosity due to shear thinning behavior. Furthermore, there is a significant difference between the results obtained by the two models for a specific patient. During the systole, differences are more pronounced and are preferably located in the tortuous regions of the artery. Throughout the cardiac cycle, the WSS amplitude between the systole and diastole is greater for the Walburn-Schneck's model than for the Newtonian model. However, the average viscosity along the artery is always greater for the non-Newtonian model, except in the systolic peak. The hemodynamic model is crucial to validate results obtained with CFD and to explore clinical potential.

Pulsatile magneto-hydrodynamic blood flows through porous blood vessels using a third grade non-Newtonian fluids model

In this paper, the unsteady pulsatile magneto-hydrodynamic blood flows through porous arteries concerning the influence of externally imposed periodic body acceleration and a periodic pressure gradient are numerically simulated. Blood is taken into account as the third-grade non-Newtonian fluid. Besides the numerical solution, for small Womersley parameter (such as blood flow through arterioles and capillaries), the analytical perturbation method is used to solve the nonlinear governing equations. Consequently, analytical expressions for the velocity profile, wall shear stress, and blood flow rate are obtained. Excellent agreement between the analytical and numerical predictions is evident. Also, the effects of body acceleration, magnetic field, third-grade non-Newtonian parameter, pressure gradient, and porosity on the flow behaviors are examined. Some important conclusions are that, when the Womersley parameter is low, viscous forces tend to dominate the flow, velocity profiles are parabolic in shape, and the center-line velocity oscillates in phase with the driving pressure gradient. In addition, by increasing the pressure gradient, the mean value of the velocity profile increases and the amplitude of the velocity remains constant. Also, when non-Newtonian effect increases, the amplitude of the velocity profile.

The numerical analysis of non-Newtonian blood flow in human patient-specific left ventricle

Recently, various non-invasive tools such as the magnetic resonance image (MRI), ultrasound imaging (USI), computed tomography (CT), and the computational fluid dynamics (CFD) have been widely utilized to enhance our current understanding of the physiological parameters that affect the initiation and the progression of the cardiovascular diseases (CVDs) associated with heart failure (HF). In particular, the hemodynamics of left ventricle (LV) has attracted the attention of the researchers due to its significant role in the heart functionality. In this study, CFD owing its capability of predicting detailed flow field was adopted to model the blood flow in images-based patient-specific LV over cardiac cycle. In most published studies, the blood is modeled as Newtonian that is not entirely accurate as the blood viscosity varies with the shear rate in non-linear manner. In this paper, we studied the effect of Newtonian assumption on the degree of accuracy of intraventricular hemodynamics. In doing so, various non-Newtonian models and Newtonian model are used in the analysis of the intraventricular flow and the viscosity of the blood. Initially, we used the cardiac MRI images to reconstruct the time-resolved geometry of the patient-specific LV. After the unstructured mesh generation, the simulations were conducted in the CFD commercial solver FLUENT to analyze the intraventricular hemodynamic parameters. The findings indicate that the Newtonian assumption cannot adequately simulate the flow dynamic within the LV over the cardiac cycle, which can be attributed to the pulsatile and recirculation nature of the flow and the low blood shear rate.