Improving Blood Flow Simulations by Incorporating Measured Subject-Specific Wall Motion

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.

Analysis of Geometric and Hemodynamic Profiles in Rat Arteriovenous Fistula Following PDE5A Inhibition

Arteriovenous fistula (AVF) is essential for chronic kidney disease (CKD) patients on hemodialysis, but treatment for AVF maturation failure remains an unmet clinical need. Successful AVF remodeling occurs through sufficient lumen expansion to increase AVF blood flow and lumen area. Aberrant blood flow is thought to impair AVF remodeling, but previous literature has largely focused on hemodynamics averaged over the entire AVF or at a single location. We hypothesized that hemodynamics is heterogeneous, and thus any treatment's effect size is heterogeneous in the AVF. To test our hypothesis, we used the PDE5A inhibitor sildenafil to treat AVFs in a rat model and performed magnetic resonance imaging (MRI) based computational fluid dynamics (CFD) to generate a detailed spatial profile of hemodynamics in AVFs. 90 mg/kg of sildenafil was administered to rats in their drinking water for 14 days. On day 14 femoral AVFs were created in rats and sildenafil treatment continued for another 21 days. 21 days post-AVF creation, rats underwent non-contrast MRI for CFD and geometrical analysis. Lumen cross-sectional area (CSA) and flow rate were used to quantify AVF remodeling. Parameters used to describe aberrant blood flow include velocity magnitude, wall shear stress (WSS), oscillatory shear index (OSI), and vorticity. Geometrical parameters include arterial-venous (A-V) distance, anastomosis angle, tortuosity, and nonplanarity angle magnitude. When averaged across the entire AVF, sildenafil treated rats had significantly higher CSA, flow rate, velocity, WSS, OSI, and vorticity than control rats. To analyze heterogeneity, the vein was separated into zones: 0-5, 5-10, 10-15, and 15-20 mm from the anastomosis. In both groups: 1) CSA increased from the 0-5 to 15-20 zone; 2) velocity, WSS, and vorticity were highest in the 0-5 zone and dropped significantly thereafter; and 3) OSI increased at the 5-10 zone and then decreased gradually. Thus, the effect size of sildenafil on AVF remodeling and the relationship between hemodynamics and AVF remodeling depend on location. There was no significant difference between control and sildenafil groups for the other geometric parameters. Rats tolerated sildenafil treatment well, and our results suggest that sildenafil may be a safe and effective therapy for AVF maturation.

Numerical analysis of the hemodynamics of rat aorta based on magnetic resonance imaging and fluid-structure interaction

Murine models have been widely used to investigate the mechanobiology of aortic atherosclerosis and dissections, which develop preferably at different anatomic locations of aorta. Based MRI and finite element analysis with fluid-structure interaction, we numerically investigated factors that may affect the blood flow and structural mechanics of rat aorta. The results indicated that aortic root motion greatly increases time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), displacement of the aorta, and enhances helical flow pattern but has limited influence on effective stress, which is highly modulated by blood pressure. Moreover, the influence of the motion component on these indicators is different with axial motion more obvious than planar motion. Surrounding fixation of the intercostal arteries and the branch vessels on aortic arch would reduce the influence of aortic root motion. The compliance of the aorta has different influences at different regions, leading to decrease in TAWSS and helical flow, increase in OSI, RRT at the aortic arch, but has reversed effects on the branch vessels. When compared with the steady flow, the pulsatile blood flow would obviously increase the WSS, the displacement, and the effective stress in most regions. In conclusion, to accurately quantify the blood flow and structural mechanics of rat aorta, the motion of the aortic root, the compliance of aortic wall, and the pulsation of blood flow should be considered. However, when only focusing on the effective stress in rat aorta, the motion of the aortic root may be neglected.

A novel formulation for the study of the ascending aortic fluid dynamics with in vivo data

Numerical simulations to evaluate thoracic aortic hemodynamics include a computational fluid dynamic (CFD) approach or fluid-structure interaction (FSI) approach. While CFD neglects the arterial deformation along the cardiac cycle by applying a rigid wall simplification, on the other side the FSI simulation requires a lot of assumptions for the material properties definition and high computational costs. The aim of this study is to investigate the feasibility of a new strategy, based on Radial Basis Functions (RBF) mesh morphing technique and transient simulations, able to introduce the patient-specific changes in aortic geometry during the cardiac cycle. Starting from medical images, aorta models at different phases of cardiac cycle were reconstructed and a transient shape deformation was obtained by proper activating incremental RBF solutions during the simulation process. The results, in terms of main hemodynamic parameters, were compared with two performed CFD simulations for the aortic model at minimum and maximum volume. Our implemented strategy copes the actual arterial variation during cardiac cycle with high accuracy, capturing the impact of geometrical variations on fluid dynamics, overcoming the complexity of a standard FSI approach.

Comparison of respiratory motion suppression techniques for 4D flow MRI

PURPOSE

The purpose of this work was to assess the impact of respiratory motion and to compare methods for suppression of respiratory motion artifacts in 4D Flow MRI.

METHODS

A numerical 3D aorta phantom was designed based on an aorta velocity field obtained by computational fluid mechanics. Motion-distorted 4D Flow MRI measurements were simulated and several different motion-suppression techniques were evaluated: Gating with fixed acceptance window size, gating with different window sizes in inner and outer k-space, and k-space reordering. Additionally, different spatial resolutions were simulated.

RESULTS

Respiratory motion reduced the image quality. All motion-suppression techniques improved the data quality. Flow rate errors of up to 30% without gating could be reduced to less than 2.5% with the most successful motion suppression methods. Weighted gating and gating combined with k-space reordering were advantageous compared with conventional fixed-window gating. Spatial resolutions finer than the amount of accepted motion did not lead to improved results.

CONCLUSION

Respiratory motion affects 4D Flow MRI data. Several different motion suppression techniques exist that are capable of reducing the errors associated with respiratory motion. Spatial resolutions finer than the degree of accepted respiratory motion do not result in improved data quality. Magn Reson Med 78:1877-1882, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

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.

GPU-accelerated model for fast, three-dimensional fluid-structure interaction computations

In this paper we introduce a methodology for performing one-way Fluid-Structure interaction (FSI), i.e. where the motion of the wall boundaries is imposed. We use a Graphics Processing Unit (GPU) accelerated Lattice-Boltzmann Method (LBM) implementation and present an efficient workflow for embedding the moving geometry, given as a set of polygonal meshes, in the LBM computation. The proposed method is first validated in a synthetic experiment: a vessel which is periodically expanding and contracting. Next, the evaluation focuses on the 3D Peristaltic flow problem: a fluid flows inside a flexible tube, where a periodic wave-like deformation produces a fluid motion along the centerline of the tube. Different geometry configurations are used and results are compared against previously published solutions. The efficient approach leads to an average execution time of approx. one hour per computation, whereas 50% of it is required for the geometry update operations. Finally, we also analyse the effect of changing the Reynolds number on the flow streamlines: the flow regime is significantly affected by the Reynolds number.

Incidence of branching patterns variations of the arch in aortic dissection in Chinese patients

Several authors have described anatomic variations of the aortic arch in 13% to 20% of the patients who do not have aortic disease. However, few studies have evaluated these patterns in the thoracic aortic dissection (TAD). In the authors' knowledge, this is the first survey that specifically investigates the frequency of these variations in a broad, nonselected group of Chinese patients with aortic dissection. Furthermore, it compares this group with a group of patients without aortic disease.The objective of this study was to define the variation frequency of the aortic arch branches pattern using the tomographic studies of 525 Chinese patients with a diagnosis of TAD. The Stanford classification was used to set the site of the initial tear of the dissection. In addition, we performed an epidemiological analysis of the aortic arch anatomic variations in TAD, and its possible implications for surgical or endovascular treatment. The general hypothesis proposal asserted that Chinese patients with dissection of the aorta have a similar incidence of variations of the aortic arch to the patients without aortic disease.A retrospective study of cases and controls was carried out using the tomographic studies (CT) of all patients admitted to the First Affiliated Hospital of Zhengzhou University, located at Henan-China, with a confirmed diagnosis of aortic dissection from January 2012 until December 2014. The group of cases consisted of 525 patients: 374 men and 151 women, with a mean age of 52.27 years (range, 20-89). The average age of the patients with Stanford A and B aortic dissection was 49.46 and 53.67, respectively. The control group consisted of 525 unselected patients without TAD who underwent a CT scan of the chest due to other indications. This group consisted of 286 men and 239 women, with a mean age of 53.60 years (range, 18-89). All the patients with aneurysm or dissection were excluded from the control group. We performed a statistical analysis of demographic data.The study found 7 different patterns of the aortic arch on both groups of cases and controls. Within the 525 patients with TAD were observed 85 (16.19%) anatomical variations, while the control group showed 112 variations (21.33%); P = 0.033. The most common anatomical variant was the bovine arch, found in 62 (11.80%) cases of TAD compared with 77 (14.66%) in the control group; P = 0.172. Anatomical variations were observed in 14.32% of the patients with Stanford A dissection and 17.09% of the patients with Stanford B dissection; P = 0.425. Patients with Stanford A dissection showed the pattern of bovine arch in 23 (13.21%) of 174 cases. In contrast, the patients with Stanford B dissection showed it in 39 (11.11%) of 351 cases; P = 0.481. The anatomical variant defined as vertebral artery of direct origin of the aortic arch was more frequent in the patients with Stanford B dissection (5.12%). The patients with Stanford A dissection presented this pattern in 1.14% of the cases; P = 0.025. This study observed an increased frequency of aortic dissection in the subgroup from 41 to 60 years old. In the subgroup from 41 to 60 years old without TAD, a greater frequency of anatomical variations were found than in the patients with TAD (20.81% vs 14.23%; P = 0.050). The same fashion was seen in patients older than 80 years (27.27% vs 0%; P = 0.030). The anatomical variations of the aortic arch with TAD occurred in 14.97% of the male patients and 19.20% of the female patients compared to 21.67% to 20.92% in the control group; P = 0.026 and P = 0.681, respectively.The aortic arch variations were found less frequently in the TAD group than in the control group in the present Chinese series. The bovine arch was considered the variant pattern of the major frequency in the patients with TAD and the control group. The anatomical variant of 4 branches, defined as vertebral artery of direct origin of the aortic arch, was more frequent in patients with Stanford B aortic dissection than in the patients with Stanford A.This finding might show an association between the geometry of the aortic arch and the site of onset of first intimal tear of dissection.