Hemodynamics of the aorta
The hemodynamics of the aorta is an ongoing field of research in which the goal is to identify what flow patterns and subsequent forces occur within the thoracic aorta. These patterns and forces are used to identify the presence and severity of cardiovascular diseases such as aortic aneurysm and atherosclerosis.[1] Some of the methods used to study the hemodynamics of aortic flow are patient scans, computational fluid dynamics models, and particle tracking velocimetry (PTV). The information gathered through these studies can be used for surgery planning and the development of implants.[2] Greater understanding of this topic reduces mortality rates associated with cardiovascular disease.
General flow patterns
The mean velocity in the aorta varies over the cardiac cycle. During systole the mean velocity rises to a peak, then it falls during diastole. This pattern is repeated with each squeezing pulse of the heart. The highest velocities are found at the exit of the valve during systole. At this stage the majority of the flow can be described with velocity vectors normal to the entrance, but in plane velocities tangent to the flow are present.[3] As the path starts to curve in the ascending aorta, the blood towards the outside of the arch tends to rotate towards the inner wall, causing a helical pattern that is observed in most individuals. As the blood moves into the aortic arch, the area with the highest velocity tends to be on the inner wall. Helical flow within the ascending aorta and aortic arch help to reduce flow stagnation and increase oxygen transport.[4] As the blood moves into the descending aorta, rotations in the flow are less present. Physiological abnormalities due to plague formation or aneurysm lead to helical flows and high velocity flows in locations where they would not normally be present or as prominent. The abnormal high velocity areas generate a higher amount of wall shear stress than normal and contribute to stenosis and further plaque formation.[5] Abnormal helical structures expose tissue to low wall shear stresses that it would not normally experience. Simulations of these flow patterns seek to identify what normal wall shear stress conditions and helical flows are present at specific location within the aorta.
Effect of age and sex
When evaluating the significance of the hemodynamics of a patient their age and gender play a role.[6] Each individual will have specific aortic geometries but trends can be identified when observed as a group. As age increases the aortic diameter tends to increase and the peak velocity of systolic flow tends to decrease until patients reach an age greater than 60 years old.[6] Patients over the age of 60 tend to have an increase of peak systolic velocity.[6] While both sexes experience the same pattern of velocity change with age, men tend to experience a wider range and higher peak velocity with age.
Effect of diabetes
Diabetes mellitus (diabetes) is a significant risk factor for cardiovascular diseases.[7] The presence of diabetes affects the dynamic viscosity of blood and the compliance of the aortic walls.[8] The dynamic viscosity of blood where diabetes is present is higher than that of healthy blood making it slightly less resistive to flow. The Young's modulus of the aortic walls where diabetes is present is higher than that of a healthy patient making it more stiff. When comparing CFD models of normal blood and wall properties with CFD models where the blood and wall properties to replicate that of an individual with diabetes, it is found that the models with diabetes have a lower mean velocity.[8] It is also observed that the outlet velocity of the descending aorta is lower in the diabetes model.[8] The blood pressure in the diabetes model is lower than that of the control model, but the mean pressures of the entire aorta are similar between both models.[8]
Modeling of aortic flow
CFD modeling of the aorta
CFD models allow for researchers to recreate flows happening within the aorta and evaluate factors that cannot be obtained through normal patient scans. These factors include wall shear stress and helicity. These factors are then used to evaluate the progression and severity of cardiovascular diseases.
Patient specific information
In order to replicate patient-specific geometries a CT scan or an MRI is taken.[2] From this scan the inlet, various outlets, and walls can be digitally reconstructed to create the control volume. A common software used to construct the geometry and discretize the mesh is ANSYS. The inlet is identified as the cross section directly above the aortic valve. The outlets are identified as the brachiocephalic artery, left and right common carotid artery, subclavian artery, and the descending aorta.
In order to replicate the flow velocities that occur in individual patients a PC-MRI is taken. The PC-MRI can be taken be 1D, 3D, or 4D. 1D PC-MRIs only capture the velocity in one direction, typically axially with the inlet. A 4D PC-MRI can capture the axial through plane velocity, as well as orthogonal in plane velocities. Although 4D PC-MRIs provide more accurate and useful information on the flow, 1D PC-MRIs are more commonly taken and used in CFD modeling of the aorta.[9] Wall shear stress and helicity of the flow tend to be influenced by which type of velocity information os used in the model.[9]
Boundary conditions
There are a variety of flows that have been modeled and studied as the inlet boundary condition. Some simplified flows include plug flow, parabolic flow, linear shear flow, and skewed cubic flow.[2] 1D and 3D flows generated from patient scans can be used as more accurate inlet conditions.[9] 1D flows include the patient-specific variation of velocity normal to inlet. 3D flows include patient-specific velocities in the inlet plane in addition to the velocities normal to the inlet. The more accurate inlet conditions are oftentimes not used because the high acquisition time and low spatial resolution of the PC-MRIs needed.[1]
In each patient-specific model there are multiple outlets. The most common outlet boundary conditions used are the two and three-element Windkessel models.[2] The two-element Windkessel model replicates viscous resistance immediately downstream from the outlet and the three-element Windkessel model accounts for the resistance of capillaries and venous circulation.[2] Comparing the results of the two outlet conditions there is no significant difference in wall shear stress.[2] It has been found that the outlet boundary condition has an effect on less of the total flow than the inlet boundary condition.[2] Because of this, the inlet boundary condition has been of higher focus in most CFD studies.
Limitations in modeling
Because there is no standard for setting inlet boundary condition in CFD models, they have to be verified with experimental results. Those results can be obtained either through in vivo measurements using 4D PC-MRIs. 4D PC-MRI's are also limited because the acquisition time is high, the spatial and temporal resolution is low, and the signal to noise ratio is also low.[10]
Particle tracking velocimetry
Particle tracking velocimetry (PTV) allows for researchers to create an experimental setup to evaluate aortic flow patterns.
Methods
A CT scan or MRI is taken of a patient to obtain the geometry of the aorta. The information from that scan is then used to create a physical model made from a transparent silicone material.[11] The material used can either be compliant, in order to replicate the dilation of the valve, or rigid.[10][1] The working fluid within the model should have a refractive index to match that of the material used to create the model. Fluorescent tracers in the working fluid are illuminated by a laser in the volume of interest. A single high speed camera can be used to capture four separate images of the same illuminated volume using and image splitter and four mirrors.[12]
The pulsating flow of the aorta is replicated by a ventricular assist device (VAD). The VAD is driven by a pump with a waveform to replicate the systole and diastole of the flow. When the pump is running, the high speed camera collects images of the tracers within the volume under investigation. A 3D velocity profile of the volume under investigation can be created from the particle movement from frame to frame. Focusing on different control volumes within the model allow for the creation of velocity profiles at different locations within the aorta.
Application of results
PTV velocity information can be used in place of 4D PC-MRI information in a CFD model.[10] The 3D velocity information from the inlet of the PTV model can be applied as the inlet boundary condition in a CFD model. That CFD model can then solve for wall shear stresses. The velocity information from PTV can also be used to create an MRI simulation.[1] The MRI simulation can then be used to assess progression of cardiovascular diseases.
References
- Gülan, Utku; Calen, Christelle; Duru, Firat; Holzner, Markus (July 2018). "Blood flow patterns and pressure loss in the ascending aorta: A comparative study on physiological and aneurysmal conditions". Journal of Biomechanics. 76: 152–159. doi:10.1016/j.jbiomech.2018.05.033. ISSN 0021-9290. PMID 29907330.
- Madhavan, Sudharsan; Kemmerling, Erica M. Cherry (2018-05-30). "The effect of inlet and outlet boundary conditions in image-based CFD modeling of aortic flow". BioMedical Engineering OnLine. 17 (1): 66. doi:10.1186/s12938-018-0497-1. ISSN 1475-925X. PMC 5975715. PMID 29843730.
- Pirola, S.; Jarral, O. A.; O'Regan, D. P.; Asimakopoulos, G.; Anderson, J. R.; Pepper, J. R.; Athanasiou, T.; Xu, X. Y. (June 2018). "Computational study of aortic hemodynamics for patients with an abnormal aortic valve: The importance of secondary flow at the ascending aorta inlet". APL Bioengineering. 2 (2): 026101. doi:10.1063/1.5011960. hdl:10044/1/57522. ISSN 2473-2877. PMC 6481743. PMID 31069298.
- Liu, Xiao; Fan, Yubo; Deng, Xiaoyan (2009-12-24). "Effect of Spiral Flow on the Transport of Oxygen in the Aorta: A Numerical Study". Annals of Biomedical Engineering. 38 (3): 917–926. doi:10.1007/s10439-009-9878-8. ISSN 0090-6964. PMID 20033776. S2CID 11591750.
- Cecchi, Emanuele; Giglioli, Cristina; Valente, Serafina; Lazzeri, Chiara; Gensini, Gian Franco; Abbate, Rosanna; Mannini, Lucia (February 2011). "Role of hemodynamic shear stress in cardiovascular disease". Atherosclerosis. 214 (2): 249–256. doi:10.1016/j.atherosclerosis.2010.09.008. ISSN 0021-9150. PMID 20970139.
- Garcia, Julio; van der Palen, Roel L.F.; Bollache, Emilie; Jarvis, Kelly; Rose, Michael J.; Barker, Alex J.; Collins, Jeremy D.; Carr, James C.; Robinson, Joshua (2017-05-26). "Distribution of blood flow velocity in the normal aorta: Effect of age and gender". Journal of Magnetic Resonance Imaging. 47 (2): 487–498. doi:10.1002/jmri.25773. ISSN 1053-1807. PMC 5702593. PMID 28556277.
- Uwe Janka, Hans (January 1996). "Increased cardiovascular morbidity and mortality in diabetes mellitus: identification of the high risk patient". Diabetes Research and Clinical Practice. 30: S85–S88. doi:10.1016/s0168-8227(96)80043-x. ISSN 0168-8227. PMID 8964198.
- Shin, Eunji; Kim, Jung Joo; Lee, Seonjoong; Ko, Kyung Soo; Rhee, Byoung Doo; Han, Jin; Kim, Nari (2018-08-23). "Hemodynamics in diabetic human aorta using computational fluid dynamics". PLOS ONE. 13 (8): e0202671. Bibcode:2018PLoSO..1302671S. doi:10.1371/journal.pone.0202671. ISSN 1932-6203. PMC 6107202. PMID 30138473.
- Morbiducci, Umberto; Ponzini, Raffaele; Gallo, Diego; Bignardi, Cristina; Rizzo, Giovanna (January 2013). "Inflow boundary conditions for image-based computational hemodynamics: Impact of idealized versus measured velocity profiles in the human aorta". Journal of Biomechanics. 46 (1): 102–109. doi:10.1016/j.jbiomech.2012.10.012. ISSN 0021-9290. PMID 23159094.
- Gallo, Diego; Gülan, Utku; Di Stefano, Antonietta; Ponzini, Raffaele; Lüthi, Beat; Holzner, Markus; Morbiducci, Umberto (September 2014). "Analysis of thoracic aorta hemodynamics using 3D particle tracking velocimetry and computational fluid dynamics". Journal of Biomechanics. 47 (12): 3149–3155. doi:10.1016/j.jbiomech.2014.06.017. ISSN 0021-9290. PMID 25017300.
- Doyle, B. J.; Morris, L. G.; Callanan, A.; Kelly, P.; Vorp, D. A.; McGloughlin, T. M. (2008). "3D Reconstruction and Manufacture of Real Abdominal Aortic Aneurysms: From CT Scan to Silicone Model". Journal of Biomechanical Engineering. 130 (3): 034501. doi:10.1115/1.2907765. ISSN 0148-0731. PMID 18532870.
- Gülan, Utku; Lüthi, Beat; Holzner, Markus; Liberzon, Alex; Tsinober, Arkady; Kinzelbach, Wolfgang (2012-09-02). "Experimental study of aortic flow in the ascending aorta via Particle Tracking Velocimetry" (PDF). Experiments in Fluids. 53 (5): 1469–1485. Bibcode:2012ExFl...53.1469G. doi:10.1007/s00348-012-1371-8. hdl:20.500.11850/58740. ISSN 0723-4864. S2CID 53508056.