A resolved two-way coupled CFD/6-DOF approach for predicting embolus transport and the embolus-trapping efficiency of IVC filters

Simulative Investigation of Different DLD Microsystem Designs with Increased Reynolds Numbers Using a Two-Way Coupled IBM-CFD/6-DOF Approach

Deterministic lateral displacement (DLD) microsystems are suitable for the size fractionation of particle suspensions in the size range of 0.1 to 10 µm. To be able to fractionate real particles beyond a laboratory scale, these systems have to be designed for higher throughputs. High flow resistances and increasing the clogging of the systems impose substantial challenges for industrial operation. Simulative parameter studies are suitable for improving the design of the systems; for example, the position and shape of the posts. A high-resolution, two-way coupled 6-DOF CFD-DEM approach was used to study the flow and particle behavior of different post shapes (circular and triangular) and post sizes at different Reynolds numbers. The results were compared with the classical first streamline width theory. It was shown that the streamline theory does not account for all effects responsible for the separation. Furthermore, a shift in the critical particle diameter to smaller values could be obtained when increasing the Reynolds number and also when using triangular posts with reduced post sizes compared to the post spacing. These findings can help to improve the efficiency of the systems as the post spacing could be extended, thus reducing the flow resistance and the probability of clogging.

Investigation of vessel occlusion during cell seeding process

The seeding of cells into an organ is an important step in cell therapy because the final functional properties of the organ are related to the initial cell distribution throughout the organ. However, vessel occlusion is a serious problem that prevents uniform distribution of the cells in the entire organ. Understanding the mechanism of vessel occlusion can help optimize the seeding process. In this study, the vessel occlusion phenomenon under perfusion conditions during cell seeding was investigated. First, we applied a microfluidic system that enabled the observation of the occlusion events during injection. Second, we applied a multiphase numerical model that can describe the cell-cell interactions and cell-fluid interactions to investigate the vessel occlusion phenomenon during the seeding process. In particular, the effects of cell concentration and flow rate were investigated. The results indicate the importance of cell-cell interactions and cell-vessel interactions for the occurrence of vessel occlusion. In addition, it is found that the probability of occurrence of vessel occlusion increases with the increase in cell concentration and decrease in flow rate. The simulation model can help determine the optimum parameters to enhance cell seeding efficiency.

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.

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 Analysis of VenaTech Convertible Vena Cava Filter Using Computational Fluid Dynamics

The VenaTech convertible filter (VTCF) has been widely used as an inferior vena cava (IVC) filter to prevent fatal pulmonary embolism in patients. However, its hemodynamics that greatly affect the filter efficacy and IVC patency are still unclear. This paper uses computational fluid dynamics with the Carreau model to simulate the non-Newtonian blood flows around the VTCF respectively deployed in the normal, reverse and three converted states in an IVC model. The results show that the prothrombotic stagnation zones are observed downstream from the normal, reverse and small open VTCFs, with the streamwise length is nearly eight times the IVC diameter. The no-slip boundary conditions of the thin-wire VTCF arms lead to the “viscous block” effect. The viscous block accelerates the blood flow by 5–15% inside the IVC and enhances the filter wall shear stress up to nearly 20 times that of the IVC only, which contributes to clot capture and thrombus lysis. The relative flow resistance is defined to evaluate the filter-induced resistance on the IVC blood flow that can be regarded as an index of IVC patency with the filter deployment. The flow resistance of the normal VTCF deployment increases dramatically by more than 60% compared with that of the IVC only and is a little higher (6%) than that of the reverse case. As the VTCF converts to a fully open configuration, the flow resistance gradually decreases to that of no filter. This work shows that even very thin VTCF arms can result in the viscous block effect and may cause significant hemodynamic impacts on clot capture, potential thrombosis and flow impedance inside the IVC. The present study also shows that CFD is a valuable and feasible in silico tool for analyzing the IVC filter hemodynamics to complement in vivo clinical and in vitro experimental studies.

Numerical Simulation of the Rheological Behavior of Nanoparticulate Suspensions

Nanoparticles significantly alter the rheological properties of a polymer or monomeric resin with major effect on the further processing of the materials. In this matter, especially the influence of particle material and disperse properties on the viscosity is not yet understood fully, but can only be modelled to some extent empirically after extensive experimental effort. In this paper, a numerical study on an uncured monomeric epoxy resin, which is filled with boehmite nanoparticles, is presented to elucidate the working principles, which govern the rheological behavior of nanoparticulate suspensions and to simulate the suspension viscosity based on assessable material and system properties. To account for the effect of particle surface forces and hydrodynamic interactions on the rheological behavior, a resolved CFD is coupled with DEM. It can be shown that the particle interactions caused by surface forces induce velocity differences between the particles and their surrounding fluid, which result in increased drag forces and cause the additional energy dissipation during shearing. The paper points out the limits of the used simulation method and presents a correction technique with respect to the Péclet number, which broadens the range of applicability. Valuable information is gained for a future mechanistic modelling of nanoparticulate suspension viscosity by elucidating the interdependency between surface forces, shear rate and resulting drag forces on the particles.

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.

Embolus Transport Simulations with Fully Resolved Particle Surfaces

PURPOSE

There has been interest in recent work in using computational fluid dynamics with Lagrangian analysis to calculate the trajectory of emboli-like particles in the vasculature. While previous studies have provided an understanding of the hemodynamic factors determining the fates of such particles and their relationship to risk of stroke, most analyses have relied on a particle equation of motion that assumes the particle is "small" e.g., much less than the diameter of the vessel. This work quantifies the limit when a particle can no longer be considered "small".

METHODS

The motion of embolus-like particles are simulated using an overset mesh technique. This allows the fluid stresses on the particle surface to be fully resolved. Consequently, the particles can be of arbitrary size or shape. The trajectory of resolved particles and "small" particles are simulated through a patient-specific carotid artery bifurcation model with particles 500, 1000, and 2000 μm in diameter. The proportions of particles entering the internal carotid artery are treated as the outcome of the particle fate, and statistical comparisons are made to ascertain the importance of non-small particle effects.

RESULTS

For the 2000 μm embolus, the proportion of particles traveling to the internal carotid artery is 74.7 ± 1.3% (mean ± 95% confidence margin) for the "small" particle model and is 85.7 ± 5.4% for a resolved particle model. The difference is statistically significant, [Formula: see text], based on the binomial test for the particle outcomes. No statistically discernible differences are found for the smaller diameter particles.

CONCLUSIONS

Quantitative differences are observable for the 2000 μm trajectories between the "small" and resolved particle models which is a particle diameter 27% relative to the common carotid artery diameter. A fully resolved particle model ought to be considered for emboli trajectory simulations when the particle size ratio is ≳ 20%.

Embolus Analog Trajectory Paths Under Physiological Flowrates Through Patient-Specific Aortic Arch Models

Atrial fibrillation (AF) is the most common irregular heartbeat among the world's population and is a major contributor to cardiogenic embolisms and acute ischemic stroke (AIS). However, the role AF flow plays in the trajectory paths of cardiogenic emboli has not been experimentally investigated. A physiological simulation system was designed to analyze the trajectory patterns of bovine embolus analogs (EAs) (n = 720) through four patient-specific models, under three flow conditions: steady flow, normal pulsatile flow, and AF pulsatile flow. It was seen that EA trajectory paths were proportional to the percentage flowrate split of 25–31% along the branching vessels. Overall, AF flow conditions increased trajectories through the left- (LCCA) and right (RCCA)-common carotid artery by 25% with respect to normal pulsatile flow. There was no statistical difference in the distribution of clot trajectories when the clot was released from the right, left, or anterior positions. Significantly, more EAs traveled through the brachiocephalic trunk (BCT) than through the LCCA or the left subclavian. Yet of the EAs that traveled through the common carotid arteries, there was a greater affiliation toward the LCCA compared to the RCCA (p < 0.05).

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

PURPOSE

Although many previous computational fluid dynamics (CFD) studies have investigated the hemodynamics in the inferior vena cava (IVC), few studies have compared computational predictions to experimental data, and only qualitative comparisons have been made. Herein, we provide particle image velocimetry (PIV) measurements of flow in a patient-averaged IVC geometry under idealized conditions typical of those used in the preclinical evaluation of IVC filters.

METHODS

Measurements are acquired under rest and exercise flow rate conditions in an optically transparent model fabricated using 3D printing. To ensure that boundary conditions are well-defined and to make follow-on CFD validation studies more convenient, fully-developed flow is provided at the inlets (i.e., the iliac veins) by extending them with straight rigid tubing longer than the estimated entrance lengths. Velocity measurements are then obtained at the downstream end of the tubing to confirm Poiseuille inflow boundary conditions.

RESULTS

Measurements in the infrarenal IVC reveal that flow profiles are blunter in the sagittal plane (minor axis) than in the coronal plane (major axis). Peak in-plane velocity magnitudes are 4.9 cm/s and 27 cm/s under the rest and exercise conditions, respectively. Flow profiles are less parabolic and exhibit more inflection points at the higher flow rate. Bimodal velocity peaks are also observed in the sagittal plane at the elevated flow condition.

CONCLUSIONS

The IVC geometry, boundary conditions, and infrarenal velocity measurements are provided for download on a free and publicly accessible repository at https://doi.org/10.6084/m9.figshare.7198703 . These data will facilitate future CFD validation studies of idealized, in vitro IVC hemodynamics and of similar laminar flows in vascular geometries.

Reversed Auxiliary Flow to Reduce Embolism Risk During TAVI: A Computational Simulation and Experimental Study

INTRODUCTION

Endovascular treatments, such as transcatheter aortic valve implantation (TAVI), carry a risk of embolization due to debris dislodgement during various procedural steps. Although embolic filters are already available and marketed, mechanisms underlying cerebral embolism still need to be elucidated in order to further reduce cerebrovascular events.

METHODS

We propose an experimental framework with an in silico duplicate allowing release of particles at the level of the aortic valve and their subsequent capture in the supra-aortic branches, simulating embolization under constant inflow and controlled hemodynamic conditions. The effect of a simple flow modulation, consisting of an auxiliary constant flow via the right subclavian artery (RSA), on the amount of particle entering the brachiocephalic trunk was investigated. Preliminary computational fluid dynamics (CFD) simulations were performed in order to assess the minimum retrograde flow-rate from RSA required to deviate particles.

RESULTS

Our results show that a constant reversed auxiliary flow of 0.5 L/min from the RSA under a constant inflow of 4 L/min from the ascending aorta is able to protect the brachiocephalic trunk from particle embolisms. Both computational and experimental results also demonstrate that the distribution of the bulk flow dictates the distribution of the particles along the aortic branches. This effect has also shown to be independent of release location and flow rate.

CONCLUSIONS

The present study confirms that the integration of in vitro experiments and in silico analyses allows designing and benchmarking novel solutions for cerebral embolic protection during TAVI such as the proposed embo-deviation technique based on an auxiliary retrograde flow from the right subclavian artery.

Advancing Regulatory Science With Computational Modeling for Medical Devices at the FDA's Office of Science and Engineering Laboratories

Protecting and promoting public health is the mission of the U.S. Food and Drug Administration (FDA). FDA's Center for Devices and Radiological Health (CDRH), which regulates medical devices marketed in the U.S., envisions itself as the world's leader in medical device innovation and regulatory science-the development of new methods, standards, and approaches to assess the safety, efficacy, quality, and performance of medical devices. Traditionally, bench testing, animal studies, and clinical trials have been the main sources of evidence for getting medical devices on the market in the U.S. In recent years, however, computational modeling has become an increasingly powerful tool for evaluating medical devices, complementing bench, animal and clinical methods. Moreover, computational modeling methods are increasingly being used within software platforms, serving as clinical decision support tools, and are being embedded in medical devices. Because of its reach and huge potential, computational modeling has been identified as a priority by CDRH, and indeed by FDA's leadership. Therefore, the Office of Science and Engineering Laboratories (OSEL)-the research arm of CDRH-has committed significant resources to transforming computational modeling from a valuable scientific tool to a valuable regulatory tool, and developing mechanisms to rely more on digital evidence in place of other evidence. This article introduces the role of computational modeling for medical devices, describes OSEL's ongoing research, and overviews how evidence from computational modeling (i.e., digital evidence) has been used in regulatory submissions by industry to CDRH in recent years. It concludes by discussing the potential future role for computational modeling and digital evidence in medical devices.