Finite element computation of magneto-hemodynamic flow and heat transfer in a bifurcated artery with saccular aneurysm using the Carreau-Yasuda biorheological model

• Published in: Microvascular research Vol. 138; p. 104221
• Main Authors: Dubey, Ankita, B., Vasu, Bég, O. Anwar, Gorla, R.S.R.
• Format: Journal Article
• Language: English
• Published: United States Elsevier Inc 01-11-2021
• Subjects: Aorta - physiopathology; Aortic Aneurysm - physiopathology; Arterial bifurcation; Biomagnetic therapy; Blood Flow Velocity; Computer Simulation; Electric Conductivity; Energy Transfer; Finite Element Analysis; Finite element method; Hemorheology; Hot Temperature; Humans; Magnetic Fields; Magneto-hemodynamics; Models, Cardiovascular; Non-Newtonian flow; Pulsatile Flow; Regional Blood Flow; Saccular aneurysm; Skin friction; Stress, Mechanical; Time Factors; Finite element method; Biomagnetic therapy; Skin friction; Arterial bifurcation; Non-Newtonian flow; Saccular aneurysm; Magneto-hemodynamics
• ISSN: 0026-2862; 1095-9319
• DOI: 10.1016/j.mvr.2021.104221

“Existing computational fluid dynamics studies of blood flows have demonstrated that the lower wall stress and higher oscillatory shear index might be the cause of acceleration in atherogenesis of vascular walls in hemodynamics. To prevent the chances of aneurysm wall rupture in the saccular aneurysm at distal aortic bifurcation, clinical biomagnetic studies have shown that extra-corporeal magnetic fields can be deployed to regulate the blood flow. Motivated by these developments, in the current study a finite element computational fluid dynamics simulation has been conducted of unsteady two-dimensional non-Newtonian magneto-hemodynamic heat transfer in electrically conducting blood flow in a bifurcated artery featuring a saccular aneurysm. The fluid flow is assumed to be pulsatile, non-Newtonian and incompressible. The Carreau-Yasuda model is adopted for blood to mimic non-Newtonian characteristics. The transformed equations with appropriate boundary conditions are solved numerically by employing the finite element method with the variational approach in the FreeFEM++ code. Hydrodynamic and thermal characteristics are elucidated in detail for the effects of key non-dimensional parameters i.e. Reynolds number (Re = 14, 21, 100, 200), Prandtl number (Pr = 14, 21) and magnetic body force parameter (Hartmann number) (M = 0.6, 1.2, 1.5) at the aneurysm and throughout the arterial domain. The influence of vessel geometry on blood flow characteristics i.e. velocity, pressure and temperature fields are also visualized through instantaneous contour patterns. It is found that an increase in the magnetic parameter reduces the pressure but increases the skin-friction coefficient in the domain. The temperature decreases at the parent artery (inlet) and both the distant and prior artery with the increment in the Prandtl number. A higher Reynolds number also causes a reduction in velocity as well as in pressure. The blood flow shows different characteristic contours with time variation at the aneurysm as well as in the arterial segment. The novelty of the current research is therefore to present a combined approach amalgamating the Carreau-Yasuda model, heat transfer and magnetohydrodynamics with complex geometric features in realistic arterial hemodynamics with extensive visualization and interpretation, in order to generalize and extend previous studies. In previous studies these features have been considered separately and not simultaneously as in the current study. The present simulations reveal some novel features of biomagnetic hemodynamics in bifurcated arterial transport featuring a saccular aneurysm which are envisaged to be of relevance in furnishing improved characterization of the rheological biomagnetic hemodynamics of realistic aneurysmic bifurcations in clinical assessment, diagnosis and magnetic-assisted treatment of cardiovascular disease.”
•Unsteady 2D non-Newtonian magneto-hemodynamic heat transfer in electrically conducting blood flow in bifurcated artery with a saccular aneurysm.•The 6-node Taylor-Hood triangular elements have been deployed in the optimized mesh design for the robust finite element simulations.•The temperature decreases at the parent artery (inlet) and both the distant and prior artery with the increment in the Prandtl number.•A higher Reynolds number also causes a reduction in velocity as well as in pressure.•It has delivered a greater insight of blood flow behavior in saccular aneurysm in the distal aortic artery.