Stimulated by a recent controversy regarding pressure drops predicted in a giant aneurysm with a proximal stenosis, the present study sought to assess variability in the prediction of pressures and flow by a wide variety of research groups. In phase I, lumen geometry, flow rates, and fluid properties were specified, leaving each research group to choose their solver, discretization, and solution strategies. Variability was assessed by having each group interpolate their results onto a standardized mesh and centerline. For phase II, a physical model of the geometry was constructed, from which pressure and flow rates were measured. Groups repeated their simulations using a geometry reconstructed from a micro-computed tomography (CT) scan of the physical model with the measured flow rates and fluid properties. Phase I results from 25 groups demonstrated remarkable consistency in the pressure patterns, with the majority predicting peak systolic pressure drops within 8% of each other. Aneurysm sac flow patterns were more variable with only a few groups reporting peak systolic flow instabilities owing to their use of high temporal resolutions. Variability for phase II was comparable, and the median predicted pressure drops were within a few millimeters of mercury of the measured values but only after accounting for submillimeter errors in the reconstruction of the life-sized flow model from micro-CT. In summary, pressure can be predicted with consistency by CFD across a wide range of solvers and solution strategies, but this may not hold true for specific flow patterns or derived quantities. Future challenges are needed and should focus on hemodynamic quantities thought to be of clinical interest.

Heterozygous loss of the arterial-specific TGFβ type I receptor, activin receptor-like kinase 1 (ALK1; ACVRL1), causes hereditary hemorrhagic telangiectasia (HHT). HHT is characterized by development of fragile, direct connections between arteries and veins, or arteriovenous malformations (AVMs). However, how decreased ALK1 signaling leads to AVMs is unknown. To understand the cellular mis-steps that cause AVMs, we assessed endothelial cell behavior in alk1-deficient zebrafish embryos, which develop cranial AVMs. Our data demonstrate that alk1 loss has no effect on arterial endothelial cell proliferation but alters arterial endothelial cell migration within lumenized vessels. In wild-type embryos, alk1-positive cranial arterial endothelial cells generally migrate towards the heart, against the direction of blood flow, with some cells incorporating into endocardium. In alk1-deficient embryos, migration against flow is dampened and migration in the direction of flow is enhanced. Altered migration results in decreased endothelial cell number in arterial segments proximal to the heart and increased endothelial cell number in arterial segments distal to the heart. We speculate that the consequent increase in distal arterial caliber and hemodynamic load precipitates the flow-dependent development of downstream AVMs.

Wide variability exists in the prediction of intracranial aneurysm WSS. While segmentation and CFD solver techniques may be difficult to standardize across groups, our findings suggest that some of the variability in image-based CFD could be reduced by establishing guidelines for model extents, inflow rates, and blood properties, and by encouraging the reporting of normalized hemodynamic parameters.

With the increased availability of computational resources, the past decade has seen a rise in the use of computational fluid dynamics (CFD) for medical applications. There has been an increase in the application of CFD to attempt to predict the rupture of intracranial aneurysms, however, while many hemodynamic parameters can be obtained from these computations, to date, no consistent methodology for the prediction of the rupture has been identified. One particular challenge to CFD is that many factors contribute to its accuracy; the mesh resolution and spatial/temporal discretization can alone contribute to a variation in accuracy. This failure to identify the importance of these factors and identify a methodology for the prediction of ruptures has limited the acceptance of CFD among physicians for rupture prediction. The International CFD Rupture Challenge 2013 seeks to comment on the sensitivity of these various CFD assumptions to predict the rupture by undertaking a comparison of the rupture and blood-flow predictions from a wide range of independent participants utilizing a range of CFD approaches. Twenty-six groups from 15 countries took part in the challenge. Participants were provided with surface models of two intracranial aneurysms and asked to carry out the corresponding hemodynamics simulations, free to choose their own mesh, solver, and temporal discretization. They were requested to submit velocity and pressure predictions along the centerline and on specified planes. The first phase of the challenge, described in a separate paper, was aimed at predicting which of the two aneurysms had previously ruptured and where the rupture site was located. The second phase, described in this paper, aims to assess the variability of the solutions and the sensitivity to the modeling assumptions. Participants were free to choose boundary conditions in the first phase, whereas they were prescribed in the second phase but all other CFD modeling parameters were not prescribed. In order to compare the computational results of one representative group with experimental results, steady-flow measurements using particle image velocimetry (PIV) were carried out in a silicone model of one of the provided aneurysms. Approximately 80% of the participating groups generated similar results. Both velocity and pressure computations were in good agreement with each other for cycle-averaged and peak-systolic predictions. Most apparent "outliers" (results that stand out of the collective) were observed to have underestimated velocity levels compared to the majority of solutions, but nevertheless identified comparable flow structures. In only two cases, the results deviate by over 35% from the mean solution of all the participants. Results of steady CFD simulations of the representative group and PIV experiments were in good agreement. The study demonstrated that while a range of numerical schemes, mesh resolution, and solvers was used, similar flow predictions were observed in the majority of cases. To further validate the computati...

The majority of severe clinically significant forms of congenital heart disease (CHD) is associated with great artery lesions, including hypoplastic, double, right or interrupted aortic arch morphologies. While fetal and neonatal interventions are advancing, their potential ability to restore cardiac function, optimal timing, location, and intensity required for intervention remain largely unknown. We here combine computational fluid dynamics (CFD) simulations with in vivo experiments to test how individual pharyngeal arch artery hemodynamics alters as a result of local interventions to obstruct individual arch artery flow. Simulated isolated occlusions within each pharyngeal arch artery were created with image derived three-dimensional (3D) reconstructions of normal chick pharyngeal arch anatomy at Hamburger-Hamilton (HH) developmental stages HH18 and HH24. Acute flow redistributions were then computed using in vivo measured subject-specific aortic sinus inflow velocity profiles. A kinematic vascular growth-rendering algorithm was then developed and implemented to test the role of changing local wall shear stress patterns in downstream 3D morphogenesis of arch arteries. CFD simulations predicted that altered pressure gradients and flow redistributions were most sensitive to occlusion of the IVth arches. To evaluate these simulations experimentally, a novel in vivo experimental model of pharyngeal arch occlusion was developed and implemented using two-photon microscopy guided femtosecond laser based photodisruption surgery. The right IVth arch was occluded at HH18, and resulting diameter changes were followed for up to 24 hours. Pharyngeal arch diameter responses to acute hemodynamic changes were predicted qualitatively but poorly quantitatively. Chronic growth and adaptation to hemodynamic changes however were predicted in a subset of arches. Our findings suggest that this complex biodynamic process is governed through more complex forms of mechanobiological vascular growth rules. Other factors in addition to wall shear stress, or more complex WSS rules are likely important in the long-term arterial growth and patterning. Combination in-silico/experimental platforms are essential for accelerating our understanding and prediction of consequences from embryonic/fetal cardiovascular occlusions, and lay the foundation for non-invasive methods to guide CHD diagnosis and fetal intervention.

Background Hypoplastic left heart syndrome (HLHS) is a major human congenital heart defect that results in single ventricle physiology and high mortality. Clinical data indicate that intracardiac blood flow patterns during cardiac morphogenesis are a significant etiology. We used the left atrial ligation (LAL) model in the chick embryo to test the hypothesis that LAL immediately alters intracardiac flow streams and the biomechanical environment, preceding morphologic and structural defects observed in HLHS. Results Using fluorescent dye injections, we found that intracardiac flow patterns from the right common cardinal vein, right vitelline vein, and left vitelline vein were altered immediately following LAL. Furthermore, we quantified a significant ventral shift of the right common cardinal and right vitelline vein flow streams. We developed an in silico model of LAL, which revealed that wall shear stress was reduced at the left atrioventricular canal and left side of the common ventricle. Conclusions Our results demonstrate that intracardiac flow patterns change immediately following LAL, supporting the role of hemodynamics in the progression of HLHS. Sites of reduced WSS revealed by computational modeling are commonly affected in HLHS, suggesting that changes in the biomechanical environment may lead to abnormal growth and remodeling of left heart structures.

The pathology of the human abdominal aortic aneurysm (AAA) and its relationship to the later complication of intraluminal thrombus (ILT) formation remains unclear. The hemodynamics in the diseased abdominal aorta are hypothesized to be a key contributor to the formation and growth of ILT. The objective of this investigation is to establish a reliable 3D flow visualization method with corresponding validation tests with high confidence in order to provide insight into the basic hemodynamic features for a better understanding of hemodynamics in AAA pathology and seek potential treatment for AAA diseases. A stereoscopic particle image velocity (PIV) experiment was conducted using transparent patient-specific experimental AAA models (with and without ILT) at three axial planes. Results show that before ILT formation, a 3D vortex was generated in the AAA phantom. This geometry-related vortex was not observed after the formation of ILT, indicating its possible role in the subsequent appearance of ILT in this patient. It may indicate that a longer residence time of recirculated blood flow in the aortic lumen due to this vortex caused sufficient shear-induced platelet activation to develop ILT and maintain uniform flow conditions. Additionally, two computational fluid dynamics (CFD) modeling codes (Fluent and an in-house cardiovascular CFD code) were compared with the two-dimensional, threecomponent velocity stereoscopic PIV data. Results showed that correlation coefficients of the out-of-plane velocity data between PIV and both CFD methods are greater than 0.85, demonstrating good quantitative agreement. The stereoscopic PIV study can be utilized as test case templates for ongoing efforts in cardiovascular CFD solver development. Likewise, it is envisaged that the patient-specific data may provide a benchmark for further studying hemodynamics of actual AAA, ILT, and their convolution effects under physiological conditions for clinical applications.

No abstract