A comparative study of altered hemodynamics in iliac vein compression syndrome.

Introduction: Iliac vein compression syndrome (IVCS) is present in over 20% of the population and is associated with left leg pain, swelling, and thrombosis. IVCS symptoms are thought to be induced by altered pelvic hemodynamics, however, there currently exists a knowledge gap on the hemodynamic differences between IVCS and healthy patients. To elucidate those differences, we carried out a patient-specific, computational modeling comparative study. Methods: Computed tomography and ultrasound velocity and area data were used to build and validate computational models for a cohort of IVCS (N = 4, Subject group) and control (N = 4, Control group) patients. Flow, cross-sectional area, and shear rate were compared between the right common iliac vein (RCIV) and left common iliac vein (LCIV) for each group and between the Subject and Control groups for the same vessel. Results: For the IVCS patients, LCIV mean shear rate was higher than RCIV mean shear rate (550 ± 103 s-1 vs. 113 ± 48 s-1, p = 0.0009). Furthermore, LCIV mean shear rate was higher in the Subject group than in the Control group (550 ± 103 s-1 vs. 75 ± 37 s-1, p = 0.0001). Lastly, the LCIV/RCIV shear rate ratio was 4.6 times greater in the Subject group than in the Control group (6.56 ± 0.9 vs. 1.43 ± 0.6, p = 0.00008). Discussion: Our analyses revealed that IVCS patients have elevated shear rates which may explain a higher thrombosis risk and suggest that their thrombus initiation process may share aspects of arterial thrombosis. We have identified hemodynamic metrics that revealed profound differences between IVCS patients and Controls, and between RCIV and LCIV in the IVCS patients. Based on these metrics, we propose that non-invasive measurement of shear rate may aid with stratification of patients with moderate compression in which treatment is highly variable. More investigation is needed to assess the prognostic value of shear rate and shear rate ratio as clinical metrics and to understand the mechanisms of thrombus formation in IVCS patients.

An ultrasound imaging and computational fluid dynamics protocol to assess hemodynamics in iliac vein compression syndrome.

OBJECTIVE: Elevated shear rates are known to play a role in arterial thrombosis; however, shear rates have not been thoroughly investigated in patients with iliac vein compression syndrome (IVCS) owing to imaging limitations and assumptions on the low shear nature of venous flows. This study was undertaken to develop a standardized protocol that quantifies IVCS shear rates and can aid in the diagnosis and treatment of patients with moderate yet symptomatic compression. METHODS: Study patients with and without IVCS had their iliac vein hemodynamics measured via duplex ultrasound (US) at two of the following three vessel locations: infrarenal inferior vena cava (IVC), right common iliac vein, and left common iliac vein, in addition to acquiring data at the right and left external iliac veins. US velocity spectra were multiplied by a weighted cross-sectional area calculated from US and computed tomography (CT) data to create flow waveforms. Flow waveforms were then scaled to enforce conservation of flow across the IVC and common iliac veins. A three-dimensional (3D), patient-specific model of the iliac vein anatomy was constructed from CT and US examination. Flow waveforms and the 3D model were used as a basis to run a computational fluid dynamics (CFD) simulation. Owing to collateral vessel flow and discrepancies between CT and US area measurements, flows in internal iliac veins and cross-sectional areas of the common iliac veins were calibrated iteratively against target common iliac flow. Simulation results on mean velocity were validated against US data at measurement locations. Simulation results were postprocessed to derive spatial and temporal values of quantities such as velocity and shear rate. RESULTS: Using our modeling protocol, we were able to build CFD models of the iliac veins that matched common iliac flow splits within 2% and measured US velocities within 10%. Proof-of-concept analyses (1 subject, 1 control) have revealed that patients with IVCS may experience elevated shear rates in the compressed left common iliac vein, more typical of the arterial rather than the venous circulation. These results encourage us to extend this protocol to a larger group of patients with IVCS and controls. CONCLUSIONS: We developed a protocol that obtains hemodynamic measurements of the IVC and iliac veins from US, creates patient-specific 3D reconstructions of the venous anatomy using CT and US examinations, and computes shear rates using calibrated CFD methods. Proof-of-concept results have indicated that patients with IVCS may experience elevated shear rates in the compressed left common iliac vein. Larger cohorts are needed to assess the relationship between venous compression and shear rates in patients with IVCS as compared with controls with noncompressed iliac veins. Further studies using this protocol may also give promising insights into whether or not to treat patients with moderate, yet symptomatic compression.

Patient-Specific Computational Modeling of Different Cannulation Strategies for Extracorporeal Membrane Oxygenation.

Institution of extracorporeal membrane oxygenation (ECMO) results in unique blood flow characteristics to the end-organ vascular beds. We studied the interplay between cardiac-driven and extracorporeal membrane oxygenation (ECMO)-driven flow to vascular beds in different ECMO configurations using a patient-specific computational fluid dynamics (CFD) analysis. A computational ECMO model (femoral artery cannulation [FAC]) was constructed using patient-specific imaging and hemodynamic data. Following model calibration, we augmented the 3D geometrical model to represent alternative ECMO configurations (ascending aorta cannulation [AAC] and subclavian artery cannulation [SAC]). We performed CFD analyses, including a novel virtual color-dye analysis to compare global and regional blood flow and pressure characteristics as well as contributions of cardiac and ECMO-derived flow to the various vascular beds. Flow waveforms at all the aortic branch vessels were pulsatile, despite low cardiac output and predominant nonpulsatile ECMO-driven hemodynamics. Virtual color-dye analysis revealed differential contribution of cardiac and ECMO-derived flow to the end-organ vascular beds in the FAC model, while this was more evenly distributed in the AAC and SAC models. While global hemodynamics were relatively similar between various ECMO configurations, several distinct hemodynamic indices, in particular wall shear stress and oscillatory shear patterns, as well as differential contribution of ECMO-derived flow to various vascular beds, showed remarkable differences. The clinical impact of this study highlighting the relevance of CFD modeling in assessment of complex hemodynamics in ECMO warrants further evaluation.

Effects of non-Newtonian viscosity on arterial and venous flow and transport.

It is well known that blood exhibits non-Newtonian viscosity, but it is generally modeled as a Newtonian fluid. However, in situations of low shear rate, the validity of the Newtonian assumption is questionable. In this study, we investigated differences between Newtonian and non-Newtonian hemodynamic metrics such as velocity, vorticity, and wall shear stress. In addition, we investigated cardiovascular transport using two different approaches, Eulerian mass transport and Lagrangian particle tracking. Non-Newtonian solutions revealed important differences in both hemodynamic and transport metrics relative to the Newtonian model. Most notably for the hemodynamic metrics, in-plane velocity and vorticity were consistently larger in the Newtonian approximation for both arterial and venous flows. Conversely, wall shear stresses were larger for the non-Newtonian case for both the arterial and venous models. Our results also indicate that for the Lagrangian metrics, the history of accumulated shear was consistently larger for both arterial and venous flows in the Newtonian approximation. Lastly, our results also suggest that the Newtonian model produces larger near wall and luminal mass transport values compared to the non-Newtonian model, likely due to the increased vorticity and recirculation. These findings demonstrate the importance of accounting for non-Newtonian behavior in cardiovascular flows exhibiting significant regions of low shear rate and recirculation.

CRIMSON: An open-source software framework for cardiovascular integrated modelling and simulation.

In this work, we describe the CRIMSON (CardiovasculaR Integrated Modelling and SimulatiON) software environment. CRIMSON provides a powerful, customizable and user-friendly system for performing three-dimensional and reduced-order computational haemodynamics studies via a pipeline which involves: 1) segmenting vascular structures from medical images; 2) constructing analytic arterial and venous geometric models; 3) performing finite element mesh generation; 4) designing, and 5) applying boundary conditions; 6) running incompressible Navier-Stokes simulations of blood flow with fluid-structure interaction capabilities; and 7) post-processing and visualizing the results, including velocity, pressure and wall shear stress fields. A key aim of CRIMSON is to create a software environment that makes powerful computational haemodynamics tools accessible to a wide audience, including clinicians and students, both within our research laboratories and throughout the community. The overall philosophy is to leverage best-in-class open source standards for medical image processing, parallel flow computation, geometric solid modelling, data assimilation, and mesh generation. It is actively used by researchers in Europe, North and South America, Asia, and Australia. It has been applied to numerous clinical problems; we illustrate applications of CRIMSON to real-world problems using examples ranging from pre-operative surgical planning to medical device design optimization.

Numerical considerations for advection-diffusion problems in cardiovascular hemodynamics.

Numerical simulations of cardiovascular mass transport pose significant challenges due to the wide range of Péclet numbers and backflow at Neumann boundaries. In this paper we present and discuss several numerical tools to address these challenges in the context of a stabilized finite element computational framework. To overcome numerical instabilities when backflow occurs at Neumann boundaries, we propose an approach based on the prescription of the total flux. In addition, we introduce a "consistent flux" outflow boundary condition and demonstrate its superior performance over the traditional zero diffusive flux boundary condition. Lastly, we discuss discontinuity capturing (DC) stabilization techniques to address the well-known oscillatory behavior of the solution near the concentration front in advection-dominated flows. We present numerical examples in both idealized and patient-specific geometries to demonstrate the efficacy of the proposed procedures. The three contributions discussed in this paper successfully address commonly found challenges when simulating mass transport processes in cardiovascular flows.