A modified method of computed fluid dynamics simulation in abdominal aorta and visceral arteries.

PURPOSE: The flow velocity of visceral arteries was measured by 2D PCMRI to produce the patient-specific flow BC imposed on the outlets of visceral arteries in CFD simulation. This modified method aimed to improve the CFD accuracy in the abdominal aorta and visceral arteries. METHODS: A volunteer underwent non-contrast-enhanced MRA to scan the abdominal aorta and visceral arteries, and 2D PCMRI to obtain the flow velocity of the aforementioned vessels. The three-dimensional geometric model was reconstructed using the MRI scan data of the abdominal aorta and visceral arteries. The flow waveforms measured by 2D PCMRI were processed and then imposed on the aortic inlet and the outlets of all visceral arteries as the flow BC. The RCR parameters of the three elements Windkessel model were modulated and imposed on the aortic outlet. CFD simulation was run in the open-source software: svSolver. The same volunteer underwent 4D flow MRI to compare the flow field with those extracted from CFD results. RESULTS: Four specific time points in a cardiac cycle and three cross-sectional planes of aorta were selected to analyze the flow field, pressure and wall shear stress (WSS) from CFD. The flow waveforms and streamlines of CFD agreed with those of 4D flow MRI. The pressure waveforms, pressure distribution and WSS distribution from CFD conformed with the physiological condition of human body. CONCLUSION: These results suggest this modified CFD method may yield reasonable flow field, pressure and WSS in the abdominal aorta and visceral arteries.

Elongation Index as a Sensitive Measure of Cell Deformation in High-Throughput Microfluidic Systems.

One of the promising approaches for high-throughput screening of cell mechanotype is microfluidic deformability cytometry (mDC), in which the apparent deformation index (DI) of the cells stretched by extensional flow at the stagnation point of a cross-slot microchannel is measured. The DI is subject to substantial measurement errors due to cell offset from the flow centerline and velocity fluctuations in inlet channels, leading to artificial widening of DI versus cell size plots. Here, we simulated an mDC experiment using a custom computational algorithm for viscoelastic cell migration. Cell motion and deformation in a cross-slot channel was modeled for fixed or randomized values of cellular mechanical properties (diameter, shear elasticity, cortical tension) and initial cell placement, with or without sinusoidal fluctuations between the inlet velocities. Our numerical simulation indicates that mDC loses sensitivity to changes in shear elasticity when the offset distance exceeds 5 µm, and just 1% velocity fluctuation causes an 11.7% drop in the DI. The obtained relationships between the cell diameter, shear elasticity, and offset distance were used to establish a new measure of cell deformation, referred to as the "elongation index" (EI). In the randomized study, the EI scatter plots were visibly separated for the low- and high-elasticity populations of cells, with a mean of 300 and 3500 Pa, whereas the standard DI output was unable to distinguish between these two groups of cells. The successful suppression of the offset artifacts with a narrower data distribution was shown for the EI output of MCF-7 cells.

Cell trapping in Y-junction microchannels: A numerical study of the bifurcation angle effect in inertial microfluidics.

The majority of microfluidic technologies for cell sorting and isolation involve bifurcating (e.g., Y- or T-shaped junction) microchannels to trap the cells of a specific type. However, the microfluidic trapping efficiency remains low, independently of whether the cells are separated by a passive or an active sorting method. Using a custom computational algorithm, we studied the migration of separated deformable cells in a Y-junction microchannel, with a bifurcation angle ranging from 30° to 180°. Single or two cells of initially spherical shape were considered under flow conditions corresponding to inertial microfluidics. Through the numerical simulation, we identified the effects of cell size, cytoplasmic viscoelasticity, cortical tension, flow rate, and bifurcation angle on the critical separation distance for cell trapping. The results of this study show that the trapping and isolation of blood cells, and circulating tumor cells in a Y-junction microchannel was most efficient and least dependent on the flow rate at the bifurcation angle of 120°. At this angle, the trapping efficiency for white blood cells and circulating tumor cells increased, respectively, by 46% and 43%, in comparison with the trapping efficiency at 60°. The efficiency to isolate invasive tumor cells from noninvasive ones increased by 32%. This numerical study provides important design criteria to optimize microfluidic technology for deformability-based cell sorting and isolation.

A Re-Engineered Software Interface and Workflow for the Open-Source SimVascular Cardiovascular Modeling Package.

Patient-specific simulation plays an important role in cardiovascular disease research, diagnosis, surgical planning and medical device design, as well as education in cardiovascular biomechanics. simvascular is an open-source software package encompassing an entire cardiovascular modeling and simulation pipeline from image segmentation, three-dimensional (3D) solid modeling, and mesh generation, to patient-specific simulation and analysis. SimVascular is widely used for cardiovascular basic science and clinical research as well as education, following increased adoption by users and development of a GATEWAY web portal to facilitate educational access. Initial efforts of the project focused on replacing commercial packages with open-source alternatives and adding increased functionality for multiscale modeling, fluid-structure interaction (FSI), and solid modeling operations. In this paper, we introduce a major SimVascular (SV) release that includes a new graphical user interface (GUI) designed to improve user experience. Additional improvements include enhanced data/project management, interactive tools to facilitate user interaction, new boundary condition (BC) functionality, plug-in mechanism to increase modularity, a new 3D segmentation tool, and new computer-aided design (CAD)-based solid modeling capabilities. Here, we focus on major changes to the software platform and outline features added in this new release. We also briefly describe our recent experiences using SimVascular in the classroom for bioengineering education.

SimVascular: An Open Source Pipeline for Cardiovascular Simulation.

Patient-specific cardiovascular simulation has become a paradigm in cardiovascular research and is emerging as a powerful tool in basic, translational and clinical research. In this paper we discuss the recent development of a fully open-source SimVascular software package, which provides a complete pipeline from medical image data segmentation to patient-specific blood flow simulation and analysis. This package serves as a research tool for cardiovascular modeling and simulation, and has contributed to numerous advances in personalized medicine, surgical planning and medical device design. The SimVascular software has recently been refactored and expanded to enhance functionality, usability, efficiency and accuracy of image-based patient-specific modeling tools. Moreover, SimVascular previously required several licensed components that hindered new user adoption and code management and our recent developments have replaced these commercial components to create a fully open source pipeline. These developments foster advances in cardiovascular modeling research, increased collaboration, standardization of methods, and a growing developer community.

Numerical simulation of the pairwise interaction of deformable cells during migration in a microchannel.

Leukocytes and other circulating cells deform and move relatively to the channel flow in the lateral and translational directions. Their migratory property is important in immune response, hemostasis, cancer progression, delivery of nutrients, and microfluidic technologies such as cell separation and enrichment, and flow cytometry. Using our three-dimensional computational algorithm for multiphase viscoelastic flow, we have investigated the effect of pairwise interaction on the lateral and translational migration of circulating cells in a microchannel. The numerical simulation data show that when two cells with the same size and small separation distance interact, repulsive interaction take place until they reach the same lateral equilibrium position. During this process, they undergo swapping or passing, depending on the initial separation distance between each other. The threshold value of this distance increases with cell deformation, indicating that the cells experiencing larger deformation are more likely to swap. When a series of closely spaced cells with the same size are considered, they generally undergo damped oscillation in both lateral and translational directions until they reach equilibrium positions where they become evenly distributed in the flow direction (self-assembly phenomenon). A series of cells with a large lateral separation distance could collide repeatedly with each other, eventually crossing the centerline and entering the other side of the channel. For a series of cells with different deformability, more deformable cells, upon impact with less deformable cells, move to an equilibrium position closer to the centerline. The results of our study show that the bulk deformation of circulating cells plays a key role in their migration in a microchannel.