Characteristics of Transition to Turbulence in a Healthy Thoracic Aorta Using Large Eddy Simulation

Overview

Journal Sci Rep

Date 2025 Jan 25

PMID 39863653

Authors

Kuiyu Cheng

Shehnaz Akhtar

Kwan Yong Lee

Sang Wook Lee

Sang-Wook Lee

Affiliations

Soon will be listed here.

Abstract

This study employed large eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) model to investigate transitional flow characteristics in an idealized model of a healthy thoracic aorta. The OpenFOAM solver pimpleFoam was used to simulate blood flow as an incompressible Newtonian fluid, with the aortic walls treated as rigid boundaries. Simulations were conducted for 30 cardiac cycles and ensemble averaging was employed to ensure statistically reliable results. Main hemodynamic parameters, such as velocity fields, turbulence intensity turbulent kinetic energy (TKE), oscillatory shear index (OSI) and wall shear stress (WSS), were analyzed throughout the circulatory system. Through 3D computational fluid dynamics (CFD) visualization, we explained the transition from laminar to turbulent flow and its development throughout the cardiac cycle. The results demonstrated that turbulence originates in the aortic arch following the peak systole phase and further develops in the aortic arch and descending aorta during the mid-deceleration and end-systole phases, with the maximum turbulence intensity exceeding 25%. WSS reached up to 30 Pa during the peak systole, with an average WSS of 6.5 Pa across the cardiac cycle. Low and oscillatory WSS were observed during diastole which can potentially contribute to the development of vascular diseases including, aortic dissection and atherosclerosis.

References

Saqr K, Zidane I . On non-Kolmogorov turbulence in blood flow and its possible role in mechanobiological stimulation. Sci Rep. 2022; 12(1):13166. PMC: 9343407. DOI: 10.1038/s41598-022-16079-5. View

Zhu Y, Mirsadraee S, Rosendahl U, Pepper J, Xu X . Fluid-Structure Interaction Simulations of Repaired Type A Aortic Dissection: a Comprehensive Comparison With Rigid Wall Models. Front Physiol. 2022; 13:913457. PMC: 9237394. DOI: 10.3389/fphys.2022.913457. View

Chiu J, Usami S, Chien S . Vascular endothelial responses to altered shear stress: pathologic implications for atherosclerosis. Ann Med. 2008; 41(1):19-28. DOI: 10.1080/07853890802186921. View

Ku D, Giddens D, Zarins C, Glagov S . Pulsatile flow and atherosclerosis in the human carotid bifurcation. Positive correlation between plaque location and low oscillating shear stress. Arteriosclerosis. 1985; 5(3):293-302. DOI: 10.1161/01.atv.5.3.293. View

Lantz J, Gardhagen R, Karlsson M . Quantifying turbulent wall shear stress in a subject specific human aorta using large eddy simulation. Med Eng Phys. 2012; 34(8):1139-48. DOI: 10.1016/j.medengphy.2011.12.002. View

Qiao Y, Luo K, Fan J . Heat transfer mechanism in idealized healthy and diseased aortas using fluid-structure interaction method. Biomech Model Mechanobiol. 2023; 22(6):1953-1964. DOI: 10.1007/s10237-023-01745-y. View

Wen J, Yan T, Su Z, Huang H, Gao Q, Chen X . Risk evaluation of type B aortic dissection based on WSS-based indicators distribution in different types of aortic arch. Comput Methods Programs Biomed. 2022; 221:106872. DOI: 10.1016/j.cmpb.2022.106872. View

Ha H, Ziegler M, Welander M, Bjarnegard N, Carlhall C, Lindenberger M . Age-Related Vascular Changes Affect Turbulence in Aortic Blood Flow. Front Physiol. 2018; 9:36. PMC: 5788974. DOI: 10.3389/fphys.2018.00036. View

Fukuda I, Minakawa M, Fukui K, Taniguchi S, Daitoku K, Suzuki Y . Breakdown of atheromatous plaque due to shear force from arterial perfusion cannula. Ann Thorac Surg. 2007; 84(4):e17-8. DOI: 10.1016/j.athoracsur.2007.06.013. View

10.

Vasava P, Jalali P, Dabagh M, Kolari P . Finite element modelling of pulsatile blood flow in idealized model of human aortic arch: study of hypotension and hypertension. Comput Math Methods Med. 2012; 2012:861837. PMC: 3287084. DOI: 10.1155/2012/861837. View

11.

Dimakopoulos Y, Bogaerds A, Anderson P, Hulsen M, Baaijens F . Direct numerical simulation of a 2D-stented aortic heart valve at physiological flow rates. Comput Methods Biomech Biomed Engin. 2011; 15(11):1157-79. DOI: 10.1080/10255842.2011.581238. View

12.

Alimohammadi M, Agu O, Balabani S, Diaz-Zuccarini V . Development of a patient-specific simulation tool to analyse aortic dissections: assessment of mixed patient-specific flow and pressure boundary conditions. Med Eng Phys. 2013; 36(3):275-84. DOI: 10.1016/j.medengphy.2013.11.003. View

13.

Kilner P, Yang G, Mohiaddin R, Firmin D, LONGMORE D . Helical and retrograde secondary flow patterns in the aortic arch studied by three-directional magnetic resonance velocity mapping. Circulation. 1993; 88(5 Pt 1):2235-47. DOI: 10.1161/01.cir.88.5.2235. View

14.

Yamaguchi T, Kikkawa S, Yoshikawa T, Tanishita K, Sugawara M . Measurement of turbulence intensity in the center of the canine ascending aorta with a hot-film anemometer. J Biomech Eng. 1983; 105(2):177-87. DOI: 10.1115/1.3138403. View

15.

Sundin J, Engvall J, Nylander E, Ebbers T, Bolger A, Carlhall C . Improved Efficiency of Intraventricular Blood Flow Transit Under Cardiac Stress: A 4D Flow Dobutamine CMR Study. Front Cardiovasc Med. 2020; 7:581495. PMC: 7724031. DOI: 10.3389/fcvm.2020.581495. View

16.

Manchester E, Pirola S, Pirola S, Mastroiacovo G, Polvani G, Pontone G . Aortic valve neocuspidization and bioprosthetic valves: Evaluating turbulence haemodynamics. Comput Biol Med. 2024; 171:108123. DOI: 10.1016/j.compbiomed.2024.108123. View

17.

Cheng C, Tempel D, van Haperen R, van der Baan A, Grosveld F, Daemen M . Atherosclerotic lesion size and vulnerability are determined by patterns of fluid shear stress. Circulation. 2006; 113(23):2744-53. DOI: 10.1161/CIRCULATIONAHA.105.590018. View

18.

Chi Q, He Y, Luan Y, Qin K, Mu L . Numerical analysis of wall shear stress in ascending aorta before tearing in type A aortic dissection. Comput Biol Med. 2017; 89:236-247. DOI: 10.1016/j.compbiomed.2017.07.029. View

19.

Ming Tse K, Chang R, Lee H, Lim S, Venkatesh S, Ho P . A computational fluid dynamics study on geometrical influence of the aorta on haemodynamics. Eur J Cardiothorac Surg. 2012; 43(4):829-38. DOI: 10.1093/ejcts/ezs388. View

20.

Takeda R, Sato F, Yokoyama H, Sasaki K, Oshima N, Kuroda A . Investigations into the Potential of Using Open Source CFD to Analyze the Differences in Hemodynamic Parameters for Aortic Dissections (Healthy versus Stanford Type A and B). Ann Vasc Surg. 2021; 79:310-323. DOI: 10.1016/j.avsg.2021.08.007. View