Jonas Lantz

WALL SHEAR STRESS IN A SUBJECT SPECIFIC HUMAN AORTA — INFLUENCE OF FLUID-STRUCTURE INTERACTION

Vascular wall shear stress (WSS) has been correlated to the development of atherosclerosis in arteries. As WSS depends on the blood flow dynamics, it is sensitive to pulsatile effects and local changes in geometry. The aim of this study is therefore to investigate if the effect of wall motion changes the WSS or if a rigid wall assumption is sufficient. Magnetic resonance imaging (MRI) was used to acquire subject specific geometry and flow rates in a human aorta, which were used as inputs in numerical models. Both rigid wall models and fluid-structure interaction (FSI) models were considered, and used to calculate the WSS on the aortic wall. A physiological range of different wall stiffnesses in the FSI simulations was used in order to investigate its effect on the flow dynamics. MRI measurements of velocity in the descending aorta were used as validation of the numerical models, and good agreement was achieved. It was found that the influence of wall motion was low on time-averaged WSS and oscillating shear index, but when regarding instantaneous WSS values the effect from the wall motion was clearly visible. Therefore, if instantaneous WSS is to be investigated, a FSI simulation should be considered.

Large eddy simulation of LDL surface concentration in a subject specific human aorta

The development of atherosclerosis is correlated to the accumulation of lipids in the arterial wall, which, in turn, may be caused by the build-up of low-density lipoproteins (LDL) on the arterial surface. The goal of this study was to model blood flow within a subject specific human aorta, and to study how the LDL surface concentration changed during a cardiac cycle. With measured velocity profiles as boundary conditions, a scale-resolving technique (large eddy simulation, LES) was used to compute the pulsatile blood flow that was in the transitional regime. The relationship between wall shear stress (WSS) and LDL surface concentration was investigated, and it was found that the accumulation of LDL correlated well with WSS. In general, regions of low WSS corresponded to regions of increased LDL concentration and vice versa. The instantaneous LDL values changed significantly during a cardiac cycle; during systole the surface concentration was low due to increased convective fluid transport, while in diastole there was an increased accumulation of LDL on the surface. Therefore, the near-wall velocity was investigated at four representative locations, and it was concluded that in regions with disturbed flow the LDL concentration had significant temporal changes, indicating that LDL accumulation is sensitive to not only the WSS but also near-wall flow.

Numerical and experimental assessment of turbulent kinetic energy in an aortic coarctation

The turbulent blood flow through an aortic coarctation in a 63-year old female patient was studied experimentally using magnetic resonance imaging (MRI), and numerically using computational fluid dynamics (CFD), before and after catheter intervention. Turbulent kinetic energy (TKE) was computed in the numerical model using large eddy simulation and compared with direct in vivo MRI measurements. Despite the two totally different methods to obtain TKE values, both quantitative and qualitative results agreed very well. The results showed that even though both blood flow rate and Reynolds number increased after intervention, total turbulent kinetic energy levels decreased in the coarctation. Therefore, the use of the Reynolds number alone as a measure of turbulence in cardiovascular flows should be used with caution. Furthermore, the change in flow field and kinetic energy were assessed, and it was found that before intervention a jet formed in the throat of the coarctation, which impacted the arterial wall just downstream the constriction. After intervention the jet was significantly weaker and broke up almost immediately, presumably resulting in less stress on the wall. As there was a good agreement between measurements and numerical results (the increase and decrease of integrated TKE matched measurements almost perfectly while peak values differed by approximately 1mJ), the CFD results confirmed the MRI measurements while at the same time providing high-resolution details about the flow. Thus, this preliminary study indicates that MR-based TKE measurements might be useful as a diagnostic tool when evaluating intervention outcome, while the detailed numerical results might be useful for further understanding of the flow for treatment planning.

Quantifying turbulent wall shear stress in a subject specific human aorta using large eddy simulation

In this study, large-eddy simulation (LES) is employed to calculate the disturbed flow field and the wall shear stress (WSS) in a subject specific human aorta. Velocity and geometry measurements using magnetic resonance imaging (MRI) are taken as input to the model to provide accurate boundary conditions and to assure the physiological relevance. In total, 50 consecutive cardiac cycles were simulated from which a phase average was computed to get a statistically reliable result. A decomposition similar to Reynolds decomposition is introduced, where the WSS signal is divided into a pulsating part (due to the mass flow rate) and a fluctuating part (originating from the disturbed flow). Oscillatory shear index (OSI) is plotted against time-averaged WSS in a novel way, and locations on the aortic wall where elevated values existed could easily be found. In general, high and oscillating WSS values were found in the vicinity of the branches in the aortic arch, while low and oscillating WSS were present in the inner curvature of the descending aorta. The decomposition of WSS into a pulsating and a fluctuating part increases the understanding of how WSS affects the aortic wall, which enables both qualitative and quantitative comparisons.

Quantifying turbulent wall shear stress in a stenosed pipe using large eddy simulation.

Large eddy simulation was applied for flow of Re=2000 in a stenosed pipe in order to undertake a thorough investigation of the wall shear stress (WSS) in turbulent flow. A decomposition of the WSS into time averaged and fluctuating components is proposed. It was concluded that a scale resolving technique is required to completely describe the WSS pattern in a subject specific vessel model, since the poststenotic region was dominated by large axial and circumferential fluctuations. Three poststenotic regions of different WSS characteristics were identified. The recirculation zone was subject to a time averaged WSS in the retrograde direction and large fluctuations. After reattachment there was an antegrade shear and smaller fluctuations than in the recirculation zone. At the reattachment the fluctuations were the largest, but no direction dominated over time. Due to symmetry the circumferential time average was always zero. Thus, in a blood vessel, the axial fluctuations would affect endothelial cells in a stretched state, whereas the circumferential fluctuations would act in a relaxed direction.

4D flow MRI-Based pressure loss estimation in stenotic flows: Evaluation using numerical simulations

To assess how 4D flow MRI‐based pressure and energy loss estimates correspond to net transstenotic pressure gradients (TPGnet) and their dependence on spatial resolution.

Patient-Specific Simulation of Cardiac Blood Flow From High-Resolution Computed Tomography

Cardiac hemodynamics can be computed from medical imaging data, and results could potentially aid in cardiac diagnosis and treatment optimization. However, simulations are often based on simplified geometries, ignoring features such as papillary muscles and trabeculae due to their complex shape, limitations in image acquisitions, and challenges in computational modeling. This severely hampers the use of computational fluid dynamics in clinical practice. The overall aim of this study was to develop a novel numerical framework that incorporated these geometrical features. The model included the left atrium, ventricle, ascending aorta, and heart valves. The framework used image registration to obtain patient-specific wall motion, automatic remeshing to handle topological changes due to the complex trabeculae motion, and a fast interpolation routine to obtain intermediate meshes during the simulations. Velocity fields and residence time were evaluated, and they indicated that papillary muscles and trabeculae strongly interacted with the blood, which could not be observed in a simplified model. The framework resulted in a model with outstanding geometrical detail, demonstrating the feasibility as well as the importance of a framework that is capable of simulating blood flow in physiologically realistic hearts.

Quantification of turbulence and velocity in stenotic flow using spiral three-dimensional phase-contrast MRI.

Evaluate spiral three‐dimensional (3D) phase contrast MRI for the assessment of turbulence and velocity in stenotic flow.

Erratum to: Quantitative Assessment of Turbulence and Flow Eccentricity in an Aortic Coarctation: Impact of Virtual Interventions

Turbulence and flow eccentricity can be measured by magnetic resonance imaging (MRI) and may play an important role in the pathogenesis of numerous cardiovascular diseases. In the present study, we propose quantitative techniques to assess turbulent kinetic energy (TKE) and flow eccentricity that could assist in the evaluation and treatment of stenotic severities. These hemodynamic parameters were studied in a pre-treated aortic coarctation (CoA) and after several virtual interventions using computational fluid dynamics (CFD), to demonstrate the effect of different dilatation options on the flow field. Patient-specific geometry and flow conditions were derived from MRI data. The unsteady pulsatile flow was resolved by large eddy simulation (LES) including non-Newtonian blood rheology. Results showed an inverse asymptotic relationship between the total amount of TKE and degree of dilatation of the stenosis, where the pre-stenotic hypoplastic segment may limit the possible improvement by treating the CoA alone. Spatiotemporal maps of TKE and flow eccentricity could be linked to the characteristics of the post-stenotic jet, showing a versatile response between the CoA dilatations. By including these flow markers into a combined MRI-CFD intervention framework, CoA therapy has not only the possibility to produce predictions via simulation, but can also be validated pre- and immediate post treatment, as well as during follow-up studies.

Estimating the irreversible pressure drop across a stenosis by quantifying turbulence production using 4D Flow MRI

The pressure drop across a stenotic vessel is an important parameter in medicine, providing a commonly used and intuitive metric for evaluating the severity of the stenosis. However, non-invasive estimation of the pressure drop under pathological conditions has remained difficult. This study demonstrates a novel method to quantify the irreversible pressure drop across a stenosis using 4D Flow MRI by calculating the total turbulence production of the flow. Simulation MRI acquisitions showed that the energy lost to turbulence production can be accurately quantified with 4D Flow MRI within a range of practical spatial resolutions (1–3 mm; regression slope = 0.91, R2 = 0.96). The quantification of the turbulence production was not substantially influenced by the signal-to-noise ratio (SNR), resulting in less than 2% mean bias at SNR > 10. Pressure drop estimation based on turbulence production robustly predicted the irreversible pressure drop, regardless of the stenosis severity and post-stenosis dilatation (regression slope = 0.956, R2 = 0.96). In vitro validation of the technique in a 75% stenosis channel confirmed that pressure drop prediction based on the turbulence production agreed with the measured pressure drop (regression slope = 1.15, R2 = 0.999, Bland-Altman agreement = 0.75 ± 3.93 mmHg).