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

Development and characterization of a Dextran/CaCl2-based blood-mimicking fluid: a comparative study of rheological and mechanical properties in artificial erythrocyte suspensions

The development of accurate blood-mimicking fluids (BMFs) is essential for in vitro studies of blood contacting medical devices. Experimental data typically relies on single-phase glycerin/water solutions as substitutes to visualize simplified blood flow. These models are accurate only at high shear rates, limiting their applicability at lower shear rates. In this study, we investigated three potential BMFs, each composed of poly(sodium acrylate-co-acrylamide) hydrogel microparticles (beads) as artificial erythrocytes. Microbeads were produced using microfluidic systems (MFS) and were suspended in three plasma-like solutions: 10% and 50% (v/v) glycerol/water solutions and a Dextran40/CaCl2 solution. The BMFs were evaluated for their rheological and mechanical properties, including particle elasticity, sedimentation behavior, and shear flow analysis, to assess their suitability for mimicking blood. Rheometric measurements were performed at room temperature using a plate-plate configuration, measuring viscosity and shear stress for shear rates of 5-500 s-1. Atomic force microscopy (AFM) measurements were conducted to assess their mechanical response. The Dextran40/CaCl2-based BMF was identified as the most promising, demonstrating rheological and mechanical properties that closely align with those of human blood. This research offers a refined approach to developing blood analogs that better simulate the mechanical response and flow characteristics of blood for the validation and development of blood contacting medical devices.

Fluid Dynamic and in Vitro Blood Study to Understand Catheter-Related Thrombosis

PURPOSE

Central venous catheters (CVCs) provide a direct route to the venous circulation but are prone to catheter-related thrombosis (CRT). A known CRT risk factor is a high catheter-to-vein ratio (CVR), or a large catheter diameter with respect to the indwelling vein size. In this study, the CVR's effect on CVC hemodynamics and its impact on CRT is investigated with in vitro and in silico experiments.

METHODS

An in vitro flow loop is used to characterize the hemodynamics around CVCs using particle image velocimetry. In addition, CRT is investigated using an in vitro flow loop with human blood and clinical catheters. The wall shear rate of flow around the CVC is computed numerically. CVRs of 0.20, 0.33, and 0.49 and Reynolds numbers of 200, 800, and 1300 are evaluated. No flow is used through CVC lumens to model chronic indwelling catheters.

RESULTS

Results show CVR ≥ 0.33 promotes platelet-rich clot growth at the device tip and at an increased rate compared to lower CVR cases. A high wall shear rate gradient on the CVC tip and an extended wake distal to the tip exists for higher CVR cases, promoting the aggregation of platelets and subsequent stagnation for clot formation. Further, the combination of the CVR and Reynolds number are crucial to CRT potential, not the CVR alone. Specifically, thrombosis risk is increased with low (stasis driven) and/or high (platelet activation driven) flow conditions, with the CVR and CVC's geometry playing an additional role in promoting fluid mechanic driven thrombus development. A high CVR (≥ 0.33) and high flow condition (≥ 1300) results in the highest risk for clot growth at the tip of the device; other locations of the device are at risk for thrombus development in lower flow conditions, regardless of the CVR. The importance of the device geometry and flow in promoting thrombus and fibrin sheath formation is also shown for the device investigated.

CONCLUSIONS

This work demonstrates that the CVR, flow, and device geometry affect CRT. For clinical cases with CVR ≥ 0.33 and/or Re ≥ 1300, the device tip may be monitored more consistently for clot formation. Thrombosis risks remain on the entire catheter, regardless of the flow condition, for a CVR = 0.49. Device placement should be chosen carefully with respect to the combination of the Reynolds number and CVR. Further study is needed on the effect of catheterization to confirm these findings.

Numerical simulation and in vitro experimental study of the hemodynamic performance of vena cava filters with helical forms

Inferior vena cava filter (IVCF) implantation is a common method of thrombus capture. By implanting a filter in the inferior vena cava (IVC), microemboli can be effectively blocked from entering the pulmonary circulation, thereby avoiding acute pulmonary embolism (PE). Inspired by the helical flow effect in the human arterial system, we propose a helical retrievable IVCF, which, due to the presence of a helical structure inducing a helical flow pattern of blood in the region near the IVCF, can effectively avoid the deposition of microemboli in the vicinity of the IVCF while promoting the cleavage of the captured thrombus clot. It also reduces the risk of IVCF dislodging and slipping in the vessel because its shape expands in the radial direction, allowing its distal end to fit closely to the IVC wall, and because its contact structure with the inner IVC wall is curved, increasing the contact area and reducing the risk of the vessel wall being punctured by the IVCF support structure. We used ANSYS 2023 software to conduct unidirectional fluid-structure coupling simulation of four different forms of IVCF, combined with microthrombus capture experiments in vitro, to explore the impact of these four forms of IVCF on blood flow patterns and to evaluate the risk of IVCF perforation and IVCF dislocation. It can be seen from the numerical simulation results that the helical structure does have the function of inducing blood flow to undergo helical flow dynamics, and the increase in wall shear stress (WSS) brought about by this function can improve the situation of thrombosis accumulation to a certain extent. Meanwhile, the placement of IVCF will change the flow state of blood flow and lead to the deformation of blood vessels. In in vitro experiments, we found that the density of the helical support rod is a key factor affecting the thrombus trapping efficiency, and in addition, the contact area between the IVCF and the vessel wall has a major influence on the risk of IVCF displacement.

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.