Search

Pathological changes in blood flow lead to altered hemodynamic forces, which are responsible for a number of conditions related to the remodeling and regeneration of the vasculature. More specifically, wall shear stress (WSS) has been shown to be a significant hemodynamic parameter with respect to aneurysm growth and rupture, as well as plaque activation leading to increased risk of stroke. In-vivo measurement of shear stress is difficult due to the stringent requirements on spatial resolution near the wall boundaries, as well as the deviation from the commonly assumed parabolic flow behavior at the wall. In this work, we propose an experimental method of in-vitro WSS calculations from high-temporal resolution velocity distributions, which are derived from 1000 fps high-speed angiography (HSA). The high-spatial and temporal resolution of our HSA detector makes such high-resolution velocity gradient measurements feasible. Presented here is the methodology for calculation of WSS in the imaging plane, as well as initial results for a variety of vascular geometries at physiologically realistic flow rates. Further, the effect of spatial resolution on the gradient calculation is explored using CFD-derived velocity data. Such angiographic-based analysis with HSA has the potential to provide critical hemodynamic feedback in an interventional setting, with the overarching objective of supporting clinical decision-making and improving patient outcomes.

Stroke has become the main cause of mortality and morbidity in patients treated with Left Ventricular Assist Devices (LVADs). The hemodynamics of the left ventricle are altered by the implantation of an LVAD, with the increase of thrombogenic flow patterns, such as stagnation regions. Time-resolved stereo particle image velocimetry (Stereo-PIV) measurements of the flow inside a patient-specific model of the left ventricle (LV) implanted with an LVAD were performed. The effects of LVAD speed, peripheral resistance and afterload were investigated. The impact of activating the LVAD pulsatility mode (periodic speed modulation) was also evaluated. Analysis of the velocity measurements in two orthogonal planes revealed stagnation zones which may be favorable to thrombus formation. Increasing LVAD speed, despite increasing the flow rate through the inflow cannula, does not automatically result in smaller stagnation regions. These results demonstrated the strong interdependence of peripheral resistance, afterload and flow through the LVAD. As a consequence, the pulsatility mode showed very limited effect on overall flow rate. However, it did reduce the size of high stagnation areas. This study showed how LVAD speed, peripheral resistance and afterload impact the complex intraventricular flow patterns in a ventricle implanted with an LVAD and quantify their thrombogenic risk.

With increasing use of left ventricular assist devices (LVAD) it is critical to devise strategies to optimize LVAD speed while controlling mean arterial pressure (MAP) and flow according to patient physiology. The complex interdependency between LVAD speed, MAP, and flow frequently makes optimization difficult under clinical conditions. We propose a method to guide this procedure in silico, narrowing the conditions to test clinically. A computational model of the circulatory network that simulates HF and LVAD support, incorporating LVAD pressure-flow curves was applied retrospectively to anonymized patient hemodynamics data from the University of Washington Medical Center. MAP management on 61 patient-specific computational models with a target of 70 mm Hg, resulting flow for a given LVAD speed was analyzed, and compared to a target output of 5 L/min. Before performing virtual MAP management, 51% had a MAP>70 mm Hg and CO>5 L/min, and 33% had a MAP>70 mm Hg and CO

This meeting abstract was removed due to the OA licensing requirements of this journal. The full abstract is listed here : https://www.svin.org/files/SVIN_2021_Abstracts_for_Web.pdf

The flow in a spherical cavity on a curved round pipe is a canonical flow that describes well the flow inside a sidewall aneurysm on an intracranial artery. Intracranial aneurysms are often treated with a flow-diverting stent (FDS), a low-porosity metal mesh that covers the entrance to the cavity, to reduce blood flow into the aneurysm sac and exclude it from mechanical stresses imposed by the blood flow. Successful treatment is highly dependent on the degree of reduction of flow inside the cavity, and the resulting altered fluid mechanics inside the aneurysm following treatment. Using stereoscopic particle image velocimetry, we characterize the fluid mechanics in a canonical configuration representative of an intracranial aneurysm treated with a FDS: a spherical cavity on the side of a curved round pipe covered with a metal mesh formed by an actual medical FDS. This porous mesh coverage is the focus of Part 2 of the paper, characterizing the effects of parent vessel Re, De and pulsatility, Wo, on the fluid dynamics, compared with the canonical configuration with no impediments to flow into the cavity that is described in Part 1 (Chassagne et al., J. Fluid Mech., vol. 915, 2021, A123). Coverage with a FDS markedly reduces the flow Re in the aneurysmal cavity, creating a viscous-dominated flow environment despite the parent vessel Re > 100. Under steady flow conditions, the topology that forms inside the cavity is shown to be a function of the parent vessel De. At low values of De, flow enters the cavity at the leading edge and remains attached to the wall before exiting at the trailing edge, a novel behaviour that was not found under any conditions of the high-Re, unimpeded cavity flow described in Part 1. Under these conditions, flow in the cavity co-rotates with the direction of the free-stream flow, similar to Stokes flow in a cavity. As De increases, the flow along the leading edge begins to separate, and the recirculation zone grows with increasing De, until, above De ≈ 180, the flow inside the cavity is fully recirculating, counter-rotating with respect to the free-stream flow. Under pulsatile flow conditions, the vortex inside the cavity progresses through the same cycle - switching from attached and co-rotating with the free-stream flow at the beginning of the cycle (low velocity and positive acceleration) to separated and counter-rotating as De reaches a critical value. The location of separation within the harmonic cycle is shown to be a function of both De and Wo. The values of aneurysmal cavity Re based on both the average velocity and the circulation inside the cavity are shown to increase with increasing values of De, while Wo is shown to have little influence on the time-averaged metrics. As De increases, the strength of the secondary flow in the parent vessel grows, due to the inertial instability in the curved pipe, and the flow rate entering the cavity increases. Thus, the effectiveness of FDS treatment to exclude the aneurysmal cavity from the haemodynamic stresses is compromised for aneurysms located on high-curvature arteries, i.e. vessels with high De, and this can be a fluid mechanics criterion to guide treatment selection.

Flow-diverting stents (FDS) are used to treat cerebral aneurysms. They promote the formation of a stable thrombus within the aneurysmal sac and, if successful, isolate the aneurysmal dome from mechanical stresses to prevent rupture. Platelet activation, a mechanism necessary for thrombus formation, is known to respond to biomechanical stimuli, particularly to the platelets’ residence time and shear stress exposure. Currently, there is no reliable method for predicting FDS treatment outcomes, either a priori or after the procedure. Eulerian computational fluid dynamic (CFD) studies of aneurysmal flow have searched for predictors of endovascular treatment outcome; however, the hemodynamics of thrombus formation cannot be fully understood without considering the platelets’ trajectories and their mechanics-triggered activation. Lagrangian analysis of the fluid mechanics in the aneurysmal vasculature provides novel metrics by tracking the platelets’ residence time (RT) and shear history (SH). Eulerian and Lagrangian parameters are compared for 19 patient-specific cases, both pre- and post-treatment, to assess the degree of change caused by the FDS and subsequent treatment efficacy.

Aims Accurate blood pressure (BP) measurement in continuous‐flow ventricular assist device (CF‐VAD) patients is imperative to reduce stroke risk. This study assesses the accuracy of the Doppler opening pressure method compared with the gold standard arterial line method in CF‐VAD patients. Methods and results In a longitudinal cohort of HeartMate II and HVAD patients, arterial line BP and simultaneously measured Doppler opening pressure were obtained. Overall correlation, agreement between Doppler opening pressure and arterial line mean vs. systolic pressure, and the effect of arterial pulsatility on the accuracy of Doppler opening pressure were analysed. A total of 1933 pairs of Doppler opening pressure and arterial line pressure readings within 1 min of each other were identified in 154 patients (20% women, mean age 55 ± 15, 50% HeartMate II and 50% HVAD). Doppler opening pressure had good correlation with invasive mean arterial pressure (r = 0.742, P

International audience; Perianeurysmal hemodynamics play a vital role in the initiation, growth and rupture of intracranial aneurysms. In vitro investigations of aneurysmal hemodynamics are helpful to visualize and measure blood flow, and aiding surgical planning approaches. Improving in vitro model creation can improve the feasibility and accuracy of hemodynamic investigations and surgical planning, improving clinical value. In this study, in vitro models were created from three-dimensional rotational angiography (3DRA) of six patients harboring intracranial aneurysms using a multi-step process involving 3D printing, index of refraction matching and silicone casting that renders the models transparent for flow visualization. Each model was treated with the same commercially-available, patient-specific, endovascular devices (coils and/or stents). All models were scanned by synchrotron X-ray microtomography to obtain high-resolution imaging of the vessel lumen, aneurysmal sac and endovascular devices. Dimensional accuracy was compared by quantifying the differences between the microtomographic reconstructions of the fabricated phantoms and the original 3DRA obtained during patient treatment. True-scale in vitro flow phantoms were successfully created for all six patients. Optical transparency was verified by using an index of refraction matched working fluid that replicated the mechanical behavior of blood. Synchrotron imaging of vessel lumen, aneurysmal sac and endovascular devices was successfully obtained, and dimensional errors were found to be O(100 μm). The creation of dimensionally-accurate, optically-transparent flow phantoms of patient-specific intracranial aneurysms is feasible using 3D printing technology. Such models may enable in vitro investigations of aneurysmal hemodynamics to aid in treatment planning and outcome prediction to devise optimal patient-specific neurointerventional strategies.

Flow in side-wall cerebral aneurysms can be ideally modelled as the combination of flow over a spherical cavity and flow in a curved circular pipe, two canonical flows. Flow in a curved pipe is known to depend on the Dean number De, combining the effects of Reynolds number, Re, and of the curvature along the pipe centreline, κ. Pulsatility in the flow introduces a dependency on the Womersley number Wo. Using stereo PIV measurements, this study investigated the effect of these three key non-dimensional parameters, by modifying pipe curvature (De), flow-rate (Re), and pulsatility frequency (Wo), on the flow patterns in a spherical cavity. A single counter-rotating vortex was observed in the cavity for all values of pipe curvature κ and Re, for both steady and pulsatile inflow conditions. Increasing the pipe curvature impacted both the flow patterns in the pipe and the cavity, by shifting the velocity profile towards the cavity opening and increasing the flow rate into the cavity. The circulation in the cavity was found to collapse well with only the Dean number, for both steady and pulsatile inflows. For pulsatile inflow, the counter-rotating vortex was unstable and the location of its centre over time was impacted by the curvature of the pipe, as well as the Re and the Wo in the freestream. The circulation in the cavity was higher for steady inflow than for the equivalent average Reynolds and Dean number pulsatile inflow, with very limited impact of the Womersley in the range studied.

As frequency of endovascular treatments for intracranial aneurysms increases, there is a growing need to understand the mechanisms for coil embolization failure. Computational fluid dynamics (CFD) modeling often simplifies modeling the endovascular coils as a homogeneous porous medium (PM), and focuses on the vascular wall endothelium, not considering the biomechanical environment of platelets. These assumptions limit the accuracy of computations for treatment predictions. We present a rigorous analysis using X-ray microtomographic imaging of the coils and a combination of Lagrangian (platelet) and Eulerian (endothelium) metrics. Four patient-specific, anatomically accurate in vitro flow phantoms of aneurysms are treated with the same patient-specific endovascular coils. Synchrotron tomography scans of the coil mass morphology are obtained. Aneurysmal hemodynamics are computationally simulated before and after coiling, using patient-specific velocity/pressure measurements. For each patient, we analyze the trajectories of thousands of platelets during several cardiac cycles, and calculate residence times (RTs) and shear exposure, relevant to thrombus formation. We quantify the inconsistencies of the PM approach, comparing them with coil-resolved (CR) simulations, showing the under- or overestimation of key hemodynamic metrics used to predict treatment outcomes. We fully characterize aneurysmal hemodynamics with converged statistics of platelet RT and shear stress history (SH), to augment the traditional wall shear stress (WSS) on the vascular endothelium. Incorporating microtomographic scans of coil morphology into hemodynamic analysis of coiled intracranial aneurysms, and augmenting traditional analysis with Lagrangian platelet metrics improves CFD predictions, and raises the potential for understanding and clinical translation of computational hemodynamics for intracranial aneurysm treatment outcomes.

This study evaluates quantitatively the impact that intermittent Aortic Valve opening has on the thrombogenicity in the aortic arch region for patients under Left Ventricular Assist Device therapy. The influence of flow through the Aortic Valve, opening once every five cardiac cycles, on the flow patterns in the ascending aortic is measured in a patient-derived CT-image-based model, after LVAD implantation. The mechanical environment of flowing platelets is investigated, by statistical treatment of outliers in Lagrangian particle tracking, and thrombogenesis metrics (platelet residence times and activation state characterized by shear stress accumulation) are compared for the cases of no Aortic Valve opening and intermittent Aortic Valve opening. All hemodynamics metrics are improved by Aortic Valve opening, even at a reduced frequency and flow rate. Residence times of platelets or microthrombi are reduced significantly by transvalvular flow, as are the shear stress history experienced and the shear stress magnitude and gradients on the aortic root endothelium. The findings of this device-neutral study support the multiple advantages of management that enables Aortic Valve opening, providing a rationale for establishing this as a standard in long-term treatment and care for advanced heart failure patients.

This study quantifies thrombogenic potential (TP) of a wide range of left ventricular assist device (LVAD) outflow graft anastomosis angles through state-of-the-art techniques: 3D imaged-based patient-specific models created via virtual surgery and unsteady computational fluid dynamics with Lagrangian particle tracking. This study aims at clarifying the influence of a single parameter (outflow graft angle) on the thrombogenesis associated with flow patterns in the aortic root after LVAD implantation. This is an important and poorly-understood aspect of LVAD therapy, because several studies have shown strong inter and intrapatient thrombogenic variability and current LVAD implantation strategies do not incorporate outflow graft angle optimization. Accurate platelet-level investigation, enabled by statistical treatment of outliers in Lagrangian particle tracking, demonstrates a strong influence of outflow graft anastomoses angle on thrombogenicity (platelet residence times and activation state characterized by shear stress accumulation) with significantly reduced TP for acutely-angled anastomosed outflow grafts. The methodology presented in this study provides a device-neutral platform for conducting comprehensive thrombogenicity evaluation of LVAD surgical configurations, empowering optimal patient-focused surgical strategies for long-term treatment and care for advanced heart failure patients.

Left ventricular assist device (LVAD) use has continued to grow. Despite recent advances in technology, LVAD patients continue to suffer from devastating complications, including stroke and device thrombosis. Among several variables affecting thrombogenicity, we hypothesize that insertion depth of the inflow cannula into the left ventricle (LV) influences hemodynamics and thrombosis risk. Blood flow patterns were studied in a patient-derived computational model of the LV, mitral valve (MV), and LVAD inflow cannula using unsteady computational fluid dynamics (CFD). Hundreds of thousands of platelets were tracked individually, for two inflow cannula insertion depth configurations (12 mm-reduced and 27 mm-conventional) using platelet-level (Lagrangian) metrics to quantify thrombogenicity. Particularly in patients with small LV dimensions, the deeper inflow cannula insertion resulted in much higher platelet shear stress histories (SH), consistent with markedly abnormal intraventricular hemodynamics. A larger proportion of platelets in this deeper insertion configuration was found to linger in the domain for long residence times (RT) and also accumulated much higher SH. The reduced inflow depth configuration promoted LV washout and reduced platelet SH. The increase of both SH and RT in the LV demonstrates the impact of inflow cannula depth on platelet activation and increased stroke risk in these patients. Inflow cannula depth of insertion should be considered as an opportunity to optimize surgical planning of LVAD therapy.

ObjectivesTo study the correlation between wall shear stress and endothelial cell expression in a patient-specific, three-dimensional (3D)-printed model of a cerebral aneurysm.Materials and methodsA 3D-printed model of a cerebral aneurysm was created from a patient’s angiogram. After populating the model with human endothelial cells, it was exposed to media under flow for 24 hours. Endothelial cell morphology was characterized in five regions of the 3D-printed model using confocal microscopy. Endothelial cells were then harvested from distinct regions of the 3D-printed model for mRNA collection and gene analysis via quantitative polymerase chain reaction (qPCR.) Cell morphology and mRNA measurement were correlated with computational fluid dynamics simulations.ResultsThe model was successfully populated with endothelial cells, which survived under flow for 24 hours. Endothelial morphology showed alignment with flow in the proximal and distal parent vessel and aneurysm neck, but disorganization in the aneurysm dome. Genetic analysis of endothelial mRNA expression in the aneurysm dome and distal parent vessel was compared with the proximal parent vessels. ADAMTS-1 and NOS3 were downregulated in the aneurysm dome, while GJA4 was upregulated in the distal parent vessel. Disorganized morphology and decreased ADAMTS-1 and NOS3 expression correlated with areas of substantially lower wall shear stress and wall shear stress gradient in computational fluid dynamics simulations.ConclusionsCreating 3D-printed models of patient-specific cerebral aneurysms populated with human endothelial cells is feasible. Analysis of these cells after exposure to flow demonstrates differences in both cell morphology and genetic expression, which correlate with areas of differential hemodynamic stress.

Introduction: Despite significant technical advancements in the design and manufacture of Left Ventricular Assist Devices, post-implant thrombotic and thromboembolic complications continue to affect long-term outcomes. Previous efforts, aimed at optimizing pump design as a means of reducing supraphysiologic shear stresses generated within the pump and associated prothrombotic shear-mediated platelet injury, have only partially altered the device hemocompatibility. Methods: We examined hemodynamic mechanisms that synergize with hypershear within the pump to contribute to the thrombogenic potential of the overall Left Ventricular Assist Device system. Results: Numerical simulations of blood flow in differing regions of the Left Ventricular Assist Device system, that is the diseased native left ventricle, the pump inflow cannula, the impeller, the outflow graft and the anastomosed downstream aorta, reveal that prothrombotic hemodynamic conditions might occur at these specific sites. Furthermore, we show that beyond hypershear, additional hemodynamic abnormalities exist within the pump, which may elicit platelet activation, such as recirculation zones and stagnant platelet trajectories. We also provide evidences that particular Left Ventricular Assist Device implantation configurations and specific post-implant patient management strategies, such as those allowing aortic valve opening, are more hemodynamically favorable and reduce the thrombotic risk. Conclusion: We extend the perspective of pump thrombosis secondary to the supraphysiologic shear stress environment of the pump to one of Left Ventricular Assist Device system thrombosis, raising the importance of comprehensive characterization of the different prothrombotic risk factors of the total system as the target to achieve enhanced hemocompatibility and improved clinical outcomes.

The prevalence of ventricular assist device (VAD) therapy has continued to increase due to a stagnant donor supply and growing advanced heart failure (HF) population. We hypothesize that left ventricular (LV) size strongly influences biocompatibility and risk of thrombosis. Unsteady computational fluid dynamics (CFD) was used in conjunction with patient-derived computational modeling and virtual surgery with a standard, apically implanted inflow cannula. A dual-focus approach of evaluating thrombogenicity was employed: platelet-based metrics to characterize the platelet environment and flow-based metrics to investigate hemodynamics. Left ventricular end-diastolic dimensions (LVEDds) ranging from 4.5 to 6.5 cm were studied and ranked according to relative thrombogenic potential. Over 150,000 platelets were individually tracked in each LV model over 15 cardiac cycles. As LV size decreased, platelets experienced markedly increased shear stress histories (SHs), whereas platelet residence time (RT) in the LV increased with size. The complex interplay between increased SH and longer RT has profound implications on thrombogenicity, with a significantly higher proportion of platelets in small LVs having long RT times and being subjected to high SH, contributing to thrombus formation. Our data suggest that small LV size, rather than decreased VAD speed, is the primary pathologic mechanism responsible for the increased incidence of thrombosis observed in VAD patients with small LVs.

Congenital heart defects (CHDs) represent the most common form of human birth defects; approximately one-third of heart defects involve malformations of the outflow tract (OFT). Maternal diabetes increases the risk of CHD by 3–5 fold. During heart organogenesis, little is known about the effects of hyperglycemia on hemodynamics, which are critical to normal heart development. Heart development prior to septation in the chick embryo was studied under hyperglycemic conditions. Sustained hyperglycemic conditions were induced, raising the average plasma glucose concentration from 70 mg/dL to 180 mg/dL, akin to the fasting plasma glucose of a patient with diabetes. The OFTs were assessed for structural and hemodynamic alterations using optical coherence tomography (OCT), confocal microscopy, and microcomputed tomography. In hyperglycemic embryos, the endocardial cushions of the proximal OFT were asymmetric, and the OFTs curvature and torsion were significantly altered. The blood flow velocity through the OFT of hyperglycemic embryos was significantly decreased, including flow reversal in 30% of the cardiac cycle. Thus, hyperglycemia at the onset of gestation results in asymmetric proximal endocardial cushions, abnormal OFT curvature, and altered hemodynamics in the developing heart. If present in humans, these results may identify early developmental alterations that contribute to the increased risk for cardiac malformations in babies from diabetic mothers.

SUMMARY The current research presents a novel method in which blood particulates – biconcave red blood cells (RBCs) and spherical cells are modeled using isogeometric analysis, specifically Non-Uniform Rational B-Splines (NURBS) in 3-D. The use of NURBS ensures that even with a coarse representation, the geometry of the blood particulates maintains an accurate description when subjected to large deformations. The fundamental advantage of this method is the coupling of the geometrical description and the stress analysis of the cell membrane into a single, unified framework. Details on the modeling approach, implementation of boundary conditions and the membrane mechanics analysis using isogeometric modeling are presented, along with validation cases for spherical and biconcave cells. Using NURBS – based isogeometric analysis, the behavior of individual cells in fluid flow is presented and analyzed in different flow regimes using as few as 176 elements for a spherical cell and 220 elements for a biconcave RBC. This work provides a framework for modeling a large number of 3-D deformable biological cells, each with its own geometric description and membrane properties. To the best knowledge of the authors, this is the first application of the NURBS – based isogeometric analysis to model and simulate blood particulates in flow in 3D. Copyright © 2014 John Wiley & Sons, Ltd.

BackgroundComputational modeling of intracranial aneurysms provides insights into the influence of hemodynamics on aneurysm growth, rupture, and treatment outcome. Standard modeling of coiled aneurysms simplifies the complex geometry of the coil mass into a homogeneous porous medium that fills the aneurysmal sac. We compare hemodynamics of coiled aneurysms modeled from high-resolution imaging with those from the same aneurysms modeled following the standard technique, in an effort to characterize sources of error from the simplified model.MaterialsPhysical models of two unruptured aneurysms were created using three-dimensional printing. The models were treated with coil embolization using the same coils as those used in actual patient treatment and then scanned by synchrotron X-ray microtomography to obtain high-resolution imaging of the coil mass. Computational modeling of each aneurysm was performed using patient-specific boundary conditions. The coils were modeled using the simplified porous medium or by incorporating the X-ray imaged coil surface, and the differences in hemodynamic variables were assessed.ResultsX-ray microtomographic imaging of coils and incorporation into computational models were successful for both aneurysms. Porous medium calculations of coiled aneurysm hemodynamics overestimated intra-aneurysmal flow, underestimated oscillatory shear index and viscous dissipation, and over- or underpredicted wall shear stress (WSS) and WSS gradient compared with X-ray-based coiled computational fluid dynamics models.ConclusionsComputational modeling of coiled intracranial aneurysms using the porous medium approach may inaccurately estimate key hemodynamic variables compared with models incorporating high-resolution synchrotron X-ray microtomographic imaging of complex aneurysm coil geometry.

Background: As heart failure prevalence continues to increase in the setting of a static donor supply, left ventricular assist device (LVAD) therapy for end-stage heart failure continues to grow. Anecdotal evidence suggests that malalignment of the LVAD inflow cannula may increase thrombosis risk, but this effect has not been explored mechanistically or quantified statistically. Our objective is to elucidate the impact of surgical angulation of the inflow cannula on thrombogenicity. Methods and Results: Unsteady computational fluid dynamics is used in conjunction with computational modeling and virtual surgery to model flow through the left ventricle for 5 different inflow cannula angulations. We use a holistic approach to evaluate thrombogenicity: platelet-based (Lagrangian) metrics to evaluate the platelet mechanical environment, combined with flow-based (Eulerian) metrics to investigate intraventricular hemodynamics. The thrombogenic potential of each LVAD inflow cannula angulation is quantitatively evaluated based on platelet shear stress history and residence time. Intraventricular hemodynamics are strongly influenced by LVAD inflow cannula angulation. Platelet behavior indicates elevated thrombogenic potential for certain inflow cannula angles, potentially leading to platelet activation. Our analysis demonstrates that the optimal range of inflow angulation is within 0±7° of the left ventricular apical axis. Conclusions: Angulation of the inflow cannula >7° from the apical axis (axis connecting mitral valve and ventricular apex) leads to markedly unfavorable hemodynamics as determined by computational fluid dynamics. Computational hemodynamic simulations incorporating Lagrangian and Eulerian metrics are a powerful tool for studying optimization of LVAD implantation strategies, with the long-term potential of improving outcomes.

We analyzed heart wall motion and blood flow dynamics in chicken embryos using in vivo optical coherence tomography (OCT) imaging and computational fluid dynamics (CFD) embryo-specific modeling. We focused on the heart outflow tract (OFT) region of day 3 embryos, and compared normal (control) conditions to conditions after performing an OFT banding intervention, which alters hemodynamics in the embryonic heart and vasculature. We found that hemodynamics and cardiac wall motion in the OFT are affected by banding in ways that might not be intuitive a priori. In addition to the expected increase in ventricular blood pressure, and increase blood flow velocity and, thus, wall shear stress (WSS) at the band site, the characteristic peristaltic-like motion of the OFT was altered, further affecting flow and WSS. Myocardial contractility, however, was affected only close to the band site due to the physical restriction on wall motion imposed by the band. WSS were heterogeneously distributed in both normal and banded OFTs. Our results show how banding affects cardiac mechanics and can lead, in the future, to a better understanding of mechanisms by which altered blood flow conditions affect cardiac development leading to congenital heart disease.

The dynamics of red blood cells (RBCs) is one of the major aspects of the cardiovascular system that has been studied intensively in the past few decades. The dynamics of biconcave RBCs are thought to have major influences in cardiovascular diseases, the problems associated with cardiovascular assistive devices, and the determination of blood rheology and properties. This article provides an overview of the works that have been accomplished in the past few decades and aim to study the dynamics of RBCs under different flow conditions. While significant progress has been made in both experimental and numerical studies, a detailed understanding of the behavior of RBCs is still faced with many challenges. Experimentally, the size of RBCs is considered to be a major limitation that allows measurements to be performed under conditions similar to physiological conditions. In numerical computations, researchers still are working to develop a model that can cover the details of the RBC mechanics as it deforms and moves in the bloodstream. Moreover, most of reported computational models have been confined to the behavior of a single RBC in 2-dimensional domains. Advanced models are yet to be devel-oped for accurate description of RBC dynamics under physiological flow conditions in 3-dimensional regimes.

Current threshold recommendations for elective abdominal aortic aneurysm (AAA) repair are based solely on maximal AAA diameter. Peak wall stress (PWS) has been demonstrated to be a better predictor than AAA diameter of AAA rupture risk. However, PWS calculations are time-intensive, not widely available, and therefore not yet clinically practical. In addition, PWS analysis does not account for variations in wall strength between patients. We therefore sought to identify surrogate clinical markers of increased PWS and decreased aortic wall strength to better predict AAA rupture risk.Patients treated at our institution from 2001 to 2014 for ruptured AAA (rAAA) were retrospectively identified and grouped into patients with small rAAA (maximum diameter 6 cm) or large rAAA (6 cm). Patients with large (6 cm) non-rAAA were also identified sequentially from 2009 for comparison. Demographics, vascular risk factors, maximal aortic diameter, and aortic outflow occlusion (AOO) were recorded. AOO was defined as complete occlusion of the common, internal, or external iliac artery. Computational fluid dynamics and finite element analysis simulations were performed to calculate wall stress distributions and to extract PWS.We identified 61 patients with rAAA, of which 15 ruptured with AAA diameter 60 mm (small rAAA group). Patients with small rAAAs were more likely to have peripheral arterial disease (PAD) and chronic obstructive pulmonary disease (COPD) than were patients in the large non-rAAA group. Patients with small rAAAs were also more likely to have AOO compared with non-rAAAs60 mm (27% vs 8%; P = .047). Among all patients with rAAAs, those with AOO ruptured at smaller mean AAA diameters than in patients without AOO (62.1 ± 11.8 mm vs 72.5 ± 16.4 mm; P = .024). PWS calculations of a representative small rAAA and a large non-rAAA showed a substantial increase in PWS with AOO.We demonstrate that AOO, PAD, and COPD in AAA are associated with rAAAs at smaller diameters. AOO appears to increase PWS, whereas COPD and PAD may be surrogate markers of decreased aortic wall strength. We therefore recommend consideration of early, elective AAA repair in patients with AOO, PAD, or COPD to minimize risk of early rupture.

Blood flow is inherently linked to embryonic cardiac development, as haemodynamic forces exerted by flow stimulate mechanotransduction mechanisms that modulate cardiac growth and remodelling. This study evaluated blood flow in the embryonic heart outflow tract (OFT) during normal development at each stage between HH13 and HH18 in chicken embryos, in order to characterize changes in haemodynamic conditions during critical cardiac looping transformations. Two-dimensional optical coherence tomography was used to simultaneously acquire both structural and Doppler flow images, in order to extract blood flow velocity and structural information and estimate haemodynamic measures. From HH13 to HH18, peak blood flow rate increased by 2.4-fold and stroke volume increased by 2.1-fold. Wall shear rate (WSR) and lumen diameter data suggest that changes in blood flow during HH13–HH18 may induce a shear-mediated vasodilation response in the OFT. Embryo-specific four-dimensional computational fluid dynamics modelling at HH13 and HH18 complemented experimental observations and indicated heterogeneous WSR distributions over the OFT. Characterizing changes in haemodynamics during cardiac looping will help us better understand the way normal blood flow impacts proper cardiac development.

Blood flow plays a critical role in regulating embryonic cardiac growth and development, with altered flow leading to congenital heart disease. Progress in the field, however, is hindered by a lack of quantification of hemodynamic conditions in the developing heart. In this study, we present a methodology to quantify blood flow dynamics in the embryonic heart using subject-specific computational fluid dynamics (CFD) models. While the methodology is general, we focused on a model of the chick embryonic heart outflow tract (OFT), which distally connects the heart to the arterial system, and is the region of origin of many congenital cardiac defects. Using structural and Doppler velocity data collected from optical coherence tomography, we generated 4D ( $$\hbox {3D}\,+\,\hbox {time}$$ ) embryo-specific CFD models of the heart OFT. To replicate the blood flow dynamics over time during the cardiac cycle, we developed an iterative inverse-method optimization algorithm, which determines the CFD model boundary conditions such that differences between computed velocities and measured velocities at one point within the OFT lumen are minimized. Results from our developed CFD model agree with previously measured hemodynamics in the OFT. Further, computed velocities and measured velocities differ by $$

Venkat Keshav Chivukula,1 Benjamin L Krog,1,2 Jones T Nauseef,2 Michael D Henry,2 Sarah C Vigmostad1 1Department of Biomedical Engineering, 2Department of Molecular Physiology and Biophysics, Holden Comprehensive Cancer Center, University of Iowa, Seamans Center for the Engineering Arts and Sciences, Iowa City, IA, USA Abstract: Over 90% of cancer deaths result not from primary tumor development, but from metastatic tumors that arise after cancer cells circulate to distal sites via the circulatory system. While it is known that metastasis is an inefficient process, the effect of hemodynamic parameters such as fluid shear stress (FSS) on the viability and efficacy of metastasis is not well understood. Recent work has shown that select cancer cells may be able to survive and possibly even adapt to FSS in vitro. The current research seeks to characterize the effect of FSS on the mechanical properties of suspended cancer cells in vitro. Nontransformed prostate epithelial cells (PrEC LH) and transformed prostate cancer cells (PC-3) were used in this study. The Young's modulus was determined using micropipette aspiration. We examined cells in suspension but not exposed to FSS (unsheared) and immediately after exposure to high (6,400 dyn/cm2) and low (510 dyn/cm2) FSS. The PrEC LH cells were ~140% stiffer than the PC-3 cells not exposed to FSS. Post-FSS exposure, there was an increase of ~77% in Young's modulus after exposure to high FSS and a ~47% increase in Young's modulus after exposure to low FSS for the PC-3 cells. There was no significant change in the Young's modulus of PrEC LH cells post-FSS exposure. Our findings indicate that cancer cells adapt to FSS, with an increased Young's modulus being one of the adaptive responses, and that this adaptation is specific only to PC-3 cells and is not seen in PrEC LH cells. Moreover, this adaptation appears to be graded in response to the magnitude of FSS experienced by the cancer cells. This is the first study investigating the effect of FSS on the mechanical properties of cancer cells in suspension, and may provide significant insights into the mechanism by which some select cancer cells may survive in the circulation, ultimately leading to metastasis at distal sites. Our findings suggest that biomechanical analysis of cancer cells could aid in identifying and diagnosing cancer in the future. Keywords: cancer, metastasis, fluid shear stress, micropipette aspiration, elastic modulus