Steady Flow in a Patient-Averaged Inferior Vena Cava-Part I: Particle Image Velocimetry Measurements at Rest and Exercise Conditions

Clinical situations for which 3D Printing is considered an appropriate representation or extension of data contained in a medical imaging examination: vascular conditions

BACKGROUND

Medical three-dimensional (3D) printing has demonstrated utility and value in anatomic models for vascular conditions. A writing group composed of the Radiological Society of North America (RSNA) Special Interest Group on 3D Printing (3DPSIG) provides appropriateness recommendations for vascular 3D printing indications.

METHODS

A structured literature search was conducted to identify all relevant articles using 3D printing technology associated with vascular indications. Each study was vetted by the authors and strength of evidence was assessed according to published appropriateness ratings.

RESULTS

Evidence-based recommendations for when 3D printing is appropriate are provided for the following areas: aneurysm, dissection, extremity vascular disease, other arterial diseases, acute venous thromboembolic disease, venous disorders, lymphedema, congenital vascular malformations, vascular trauma, vascular tumors, visceral vasculature for surgical planning, dialysis access, vascular research/development and modeling, and other vasculopathy. Recommendations are provided in accordance with strength of evidence of publications corresponding to each vascular condition combined with expert opinion from members of the 3DPSIG.

CONCLUSION

This consensus appropriateness ratings document, created by the members of the 3DPSIG, provides an updated reference for clinical standards of 3D printing for the care of patients with vascular conditions.

Numerical simulation and in vitro experimental study of thrombus capture efficiency of a new retrievable vena cava filter

The vena cava filter is a filtering device to prevent pulmonary embolism caused by thrombosis from lower limbs and pelvis. A new retrievable vena cava filter was evaluated in this paper. To evaluate the hemodynamic performance and thrombus capture efficiency after transplantation, numerical simulation of computational fluid dynamics was performed. In this paper, the two-phase flow model of computational fluid dynamics software was used to analyze the outlet blood flow velocity, inlet-outlet pressure difference, filter wall shear stress, the ratio of area with wall shear stress, and the thrombus capture efficiency with the thrombus diameter of 5 mm, 10 mm, 15 mm and the thrombus content of 10%, 20%, 30%, respectively. Additionally, in vitro experimental test was performed to compare its thrombus capture efficiency with Denali and Aegisy Filters. The Denali Filter showed the least interference with the blood flow, followed by the new filter and the Aegisy Filter. The results indicated that the new filter had a higher capture rate in capturing 5mm small-diameter thrombus. This research certain theoretical significance and reference value for the research and development of the new vena cava filters as well as the evaluation of the thrombus capture efficiency of the filters.

Modeling Flow in an In Vitro Anatomical Cerebrovascular Model with Experimental Validation

Acute ischemic stroke (AIS) is a leading cause of mortality that occurs when an embolus becomes lodged in the cerebral vasculature and obstructs blood flow in the brain. The severity of AIS is determined by the location and how extensively emboli become lodged, which are dictated in large part by the cerebral flow and the dynamics of embolus migration which are difficult to measure in vivo in AIS patients. Computational fluid dynamics (CFD) can be used to predict the patient-specific hemodynamics and embolus migration and lodging in the cerebral vasculature to better understand the underlying mechanics of AIS. To be relied upon, however, the computational simulations must be verified and validated. In this study, a realistic in vitro experimental model and a corresponding computational model of the cerebral vasculature are established that can be used to investigate flow and embolus migration and lodging in the brain. First, the in vitro anatomical model is described, including how the flow distribution in the model is tuned to match physiological measurements from the literature. Measurements of pressure and flow rate for both normal and stroke conditions were acquired and corresponding CFD simulations were performed and compared with the experiments to validate the flow predictions. Overall, the CFD simulations were in relatively close agreement with the experiments, to within ±7% of the mean experimental data with many of the CFD predictions within the uncertainty of the experimental measurement. This work provides an in vitro benchmark data set for flow in a realistic cerebrovascular model and is a first step towards validating a computational model of AIS.

Modeling flow in an in vitro anatomical cerebrovascular model with experimental validation

Acute ischemic stroke (AIS) is a leading cause of mortality that occurs when an embolus becomes lodged in the cerebral vasculature and obstructs blood flow in the brain. The severity of AIS is determined by the location and how extensively emboli become lodged, which are dictated in large part by the cerebral flow and the dynamics of embolus migration which are difficult to measure in vivo in AIS patients. Computational fluid dynamics (CFD) can be used to predict the patient-specific hemodynamics and embolus migration and lodging in the cerebral vasculature to better understand the underlying mechanics of AIS. To be relied upon, however, the computational simulations must be verified and validated. In this study, a realistic in vitro experimental model and a corresponding computational model of the cerebral vasculature are established that can be used to investigate flow and embolus migration and lodging in the brain. First, the in vitro anatomical model is described, including how the flow distribution in the model is tuned to match physiological measurements from the literature. Measurements of pressure and flow rate for both normal and stroke conditions were acquired and corresponding CFD simulations were performed and compared with the experiments to validate the flow predictions. Overall, the CFD simulations were in relatively close agreement with the experiments, to within ±7% of the mean experimental data with many of the CFD predictions within the uncertainty of the experimental measurement. This work provides an in vitro benchmark data set for flow in a realistic cerebrovascular model and is a first step towards validating a computational model of AIS.

Experimental Hemodynamics within the Penn State Fontan Circulatory Assist Device

For children born with a single functional ventricle, the Fontan operation bypasses the right ventricle by forming a four-way total cavopulmonary connection adapting the existing ventricle for the systemic circulation. However, upon adulthood, many Fontan patients exhibit low cardiac output and elevated venous pressure, eventually requiring a heart transplantation. Despite efforts to develop a Fontan pump or use an existing ventricular assist device for failing Fontan support, there is still no device designed or tested for subpulmonary support. Penn State University is developing a hydrodynamically levitated Fontan circulatory assist device (FCAD) for bridge-to-transplant or destination therapy. The FCAD hemodynamics, at both steady and pulsatile conditions for three pump operating conditions, were quantified using particle image velocimetry to determine the velocity magnitudes and Reynolds normal and shear stresses. Data were acquired at three planes (0 mm and ±25% of the radius) for the inferior and superior vena cavae inlets and the pulmonary artery outlet. The inlets had a blunt velocity profile that became skewed towards the collecting volute as fluid approached the rotor. At the outlet, regardless of the flow condition, a high-velocity jet exited the volute and moved downstream in a helical pattern. Turbulent stresses observed at the volute exit were influenced by the rotor's rotation. Regardless of inlet conditions, the pump demonstrated advantageous behavior for clinical use with a predictable flow field and a low risk of platelet adhesion and hemolysis based on calculated wall shear rates and turbulent stresses, respectively.

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

This paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but the computational tests indicate that the methodology introduced here remains stable for selected test cases across a range of solid-fluid density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

Computational evaluation of inferior vena cava filters through computational fluid dynamics methods

Numerical simulation is growing in its importance toward the design, testing and evaluation of medical devices. Computational fluid dynamics and finite element analysis allow improved calculation of stress, heat transfer, and flow to better understand the medical device environment. Current research focuses not only on improving medical devices, but also on improving the computational tools themselves. As methods and computer technology allow for faster simulation times, iterations and trials can be performed faster to collect more data. Given the adverse events associated with long-term inferior vena cava (IVC) filter placement, IVC filter design and device evaluation are of paramount importance. This work reviews computational methods used to develop, test, and improve IVC filters to ultimately serve the needs of the patient.

Time-Resolved PIV Measurements and Turbulence Characteristics of Flow Inside an Open-Cell Metal Foam

Open-cell metal foams are porous medium for thermo-fluidic systems. However, their complex geometry makes it difficult to perform time-resolved (TR) measurements inside them. In this study, a TR particle image velocimetry (PIV) method is introduced for use inside open-cell metal foam structures. Stereolithography 3D printing methods and conventional post-processing methods cannot be applied to metal foam structures; therefore, PolyJet 3D printing and post-processing methods were employed to fabricate a transparent metal foam replica. The key to obtaining acceptable transparency in this method is the complete removal of the support material from the printing surfaces. The flow characteristics inside a 10-pore-per-inch (PPI) metal foam were analyzed in which porosity is 0.92 while laminar flow condition is applied to inlet. The flow inside the foam replica is randomly divided and combined by the interconnected pore network. Robust crosswise motion occurs inside foam with approximately 23% bulk speed. Strong influence on transverse motion by metal foam is evident. In addition, span-wise vorticity evolution is similar to the integral time length scale of the stream-wise center plane. The span-wise vorticity fluctuation through the foam arrangement is presented. It is believed that this turbulent characteristic is caused by the interaction of jets that have different flow directions inside the metal foam structure. The finite-time Lyapunov exponent method is employed to visualize the vortex ridges. Fluctuating attracting and repelling material lines are expected to enhance the heat and mass transfer. The results presented in this study could be useful for understanding the flow characteristics inside metal foams.

In Vitro Clot Trapping Efficiency of the FDA Generic Inferior Vena Cava Filter in an Anatomical Model: An Experimental Fluid-Structure Interaction Benchmark

PURPOSE

Robust experimental data for performing validation of fluid-structure interaction (FSI) simulations of the transport of deformable solid bodies in internal flow are currently lacking. This in vitro experimental study characterizes the clot trapping efficiency of a new generic conical-type inferior vena cava (IVC) filter in a rigid anatomical model of the IVC with carefully characterized test conditions, fluid rheological properties, and clot mechanical properties.

METHODS

Various sizes of spherical and cylindrical clots made of synthetic materials (nylon and polyacrylamide gel) and bovine blood are serially injected into the anatomical IVC model under worst-case exercise flow conditions. Clot trapping efficiencies and their uncertainties are then quantified for each combination of clot shape, size, and material.

RESULTS

Experiments reveal the clot trapping efficiency increases with increasing clot diameter and length, with trapping efficiencies ranging from as low as approximately 42% for small 3.2 mm diameter spherical clots up to 100% for larger clot sizes. Because of the asymmetry of the anatomical IVC model, the data also reveal the iliac vein of clot origin influences the clot trapping efficiency, with the trapping efficiency for clots injected into the left iliac vein up to a factor of 7.5 times greater than that for clots injected into the right iliac (trapping efficiencies of approximately 10% versus 75%, respectively).

CONCLUSION

Overall, this data set provides a benchmark for validating simulations predicting IVC filter clot trapping efficiency and, more generally, low-Reynolds number FSI modeling.

Hemodynamic effects of blood clots trapped by an inferior vena cava filter

The alteration of blood flow around an OPTEASE inferior vena cava filter with one or two blood clots attached was investigated by means of computational fluid dynamics. We used a patient-specific vein wall geometry, and we generated different clot models with shapes adapted to the filter and vein wall geometries. A total of eight geometries, with one or two clots and a total clot volume of 0.5 or 1 cm³ , were considered. A non-Newtonian model for blood viscosity was adopted and the possible development of turbulence was accounted for by means of a three-equation model. Two blood flow rates were considered for each case, representative for rest and exercise conditions. In exercise conditions, flow unsteadiness and even turbulence was detected in some cases. Pressure and wall shear stress (WSS) distributions were modified in all cases. Clots attached to the filter downstream basket considerably increased averaged WSS values by up to almost 50%. In all the cases a flow recirculation region appeared downstream of the clot. The degree of flow stagnation in these regions, an indicator of propensity to thrombogenesis, was estimated in terms of mean residence times and mean blood viscosity. High levels of flow stagnation were detected in rest conditions in the wake of those clots that were placed upstream from the filter. Our results suggest that one downstream placed big clot, showing a higher tendency to induce flow instabilities and turbulence, might be more harmful than two small clots placed in tandem.

Over-the-wire deployment techniques of option elite inferior vena cava filter: 3D printing vena cava phantom study

•

Bent stiff-wire technique with transfemoral access had lower filter tilt ratio at Option IVC filter deployment.

•

Original push wire with transjugular access had lower filter tilt ratio at Option IVC filter deployment.

•

Filter jumping was common using the original push wire with transjugular access.

Purpose

To compare filter tilt and filter jumping during Option inferior vena cava (IVC) filter deployment with 3 different wires techniques using a 3-dimensional (3D) printing vena cava phantom.

Materials and methods

An IVC 3D printed vena cava phantom was made from a healthy young male’s computed tomographic data. Option IVC filters were deployed with 3 different wires: i) original push wire, ii) hydrophilic stiff wire, and iii) bent stiff wire. Right internal jugular and right femoral access were used 5 times with each wire. Filter tilt angle, tilt ratio, jumping, and tip abutment to the IVC wall were analyzed.

Results

The transfemoral approach with original push wire had significantly higher tilt angle than did the transjugular approach (6.1˚ ± 1.9 vs. 3.5˚ ± 1.3, p = 0.04). Mean tilt ratio was significantly lower with the bent wire with transfemoral access (0.49 ± 0.13 vs. 0.78 ± 0.18 [original push-wire] and 0.67 ± 0.08 [stiff wire], p = 0.019). The ratio was lower also with original push wire with transjugular access (0.34 ± 0.19 vs. 0.57 ±0.11 [stiff wire] and 0.58 ±0.17 [bent wire], p = 0.045). Filter jumping occurred more often with the transjugular approach with original push wire than with stiff or bent-wire delivery. Filter tip abutment to the IVC wall occurred only with the transfemoral approach.

Conclusions

Bent wire with transfemoral access and original push wire with transjugular access had lower filter tilt ratio at Option IVC filter deployment. However, filter jumping was common using the original push wire with transjugular access.

An Immersed Interface Method for Discrete Surfaces

Fluid-structure systems occur in a range of scientific and engineering applications. The immersed boundary (IB) method is a widely recognized and effective modeling paradigm for simulating fluid-structure interaction (FSI) in such systems, but a difficulty of the IB formulation of these problems is that the pressure and viscous stress are generally discontinuous at fluid-solid interfaces. The conventional IB method regularizes these discontinuities, which typically yields low-order accuracy at these interfaces. The immersed interface method (IIM) is an IB-like approach to FSI that sharply imposes stress jump conditions, enabling higher-order accuracy, but prior applications of the IIM have been largely restricted to numerical methods that rely on smooth representations of the interface geometry. This paper introduces an immersed interface formulation that uses only a C 0 representation of the immersed interface, such as those provided by standard nodal Lagrangian finite element methods. Verification examples for models with prescribed interface motion demonstrate that the method sharply resolves stress discontinuities along immersed boundaries while avoiding the need for analytic information about the interface geometry. Our results also demonstrate that only the lowest-order jump conditions for the pressure and velocity gradient are required to realize global second-order accuracy. Specifically, we demonstrate second-order global convergence rates along with nearly second-order local convergence in the Eulerian velocity field, and between first- and second-order global convergence rates along with approximately first-order local convergence for the Eulerian pressure field. We also demonstrate approximately second-order local convergence in the interfacial displacement and velocity along with first-order local convergence in the fluid traction along the interface. As a demonstration of the method's ability to tackle more complex geometries, the present approach is also used to simulate flow in a patient-averaged anatomical model of the inferior vena cava, which is the large vein that carries deoxygenated blood from the lower extremities back to the heart. Comparisons of the general hemodynamics and wall shear stress obtained by the present IIM and a body-fitted discretization approach show that the present method yields results that are in good agreement with those obtained by the body-fitted approach.