skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Hybrid Continuum-Particle Approach for Fluid-Structure Interaction Simulation of Red Blood Cells in Fluid Flows
Transport of cells in fluid flow plays a critical role in many physiological processes of the human body. Recent developments of in vitro techniques have enabled the understanding of cellular dynamics in laboratory conditions. However, it is challenging to obtain precise characteristics of cellular dynamics using experimental method alone, especially under in vivo conditions. This challenge motivates new developments of computational methods to provide complementary data that experimental techniques are not able to provide. Since there exists a large disparity in spatial and temporal scales in this problem, which requires a large number of cells to be simulated, it is highly desirable to develop an efficient numerical method for the interaction of cells and fluid flows. In this work, a new Fluid-Structure Interaction formulation is proposed based on the use of hybrid continuum-particle approach, which can resolve local dynamics of cells while providing large-scale flow patterns in the vascular vessel. Here, the Dissipative Particle Dynamics (DPD) model for the cellular membrane is used in conjunction with the Immersed Boundary Method (IBM) for the fluid plasma. Our results show that the new formulation is highly efficient in computing the deformation of cells within fluid flow while satisfying the incompressibility constraints of the fluid. We demonstrate that it is possible to couple the DPD with the IBM to simulate the complex dynamics of Red Blood Cells (RBC) such as parachuting. Our key observation is that the proposed coupling enables the simulation of RBC dynamics in realistic arterioles while ensuring the incompressibility constraint for fluid plasma. Therefore, the proposed method allows an accurate estimation of fluid shear stresses on the surface of simulated RBC. Our results suggest that this hybrid methodology can be extended for a variety of cells in physiological conditions.  more » « less
Award ID(s):
1946202 2019077
PAR ID:
10318357
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Fluids
Volume:
6
Issue:
4
ISSN:
2311-5521
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We propose a novel solid-fluid interaction method for coupling elastic solids with impulse flow maps. Our key idea is to unify the representation of fluid and solid components as particle flow maps with different lengths and dynamics. The solid-fluid coupling is enabled by implementing two novel mechanisms: first, we developed an impulse-to-velocity transfer mechanism to unify the exchanged physical quantities; second, we devised a particle path integral mechanism to accumulate coupling forces along each flow-map trajectory. Our framework integrates these two mechanisms into an Eulerian-Lagrangian impulse fluid simulator to accommodate traditional coupling models, exemplified by the Material Point Method (MPM) and Immersed Boundary Method (IBM), within a particle flow map framework. We demonstrate our method's efficacy by simulating solid-fluid interactions exhibiting strong vortical dynamics, including various vortex shedding and interaction examples across swimming, falling, breezing, and combustion. 
    more » « less
  2. We propose a novel Particle Flow Map (PFM) method to enable accurate long-range advection for incompressible fluid simulation. The foundation of our method is the observation that a particle trajectory generated in a forward simulation naturally embodies a perfect flow map. Centered on this concept, we have developed an Eulerian-Lagrangian framework comprising four essential components: Lagrangian particles for a natural and precise representation of bidirectional flow maps; a dual-scale map representation to accommodate the mapping of various flow quantities; a particle-to-grid interpolation scheme for accurate quantity transfer from particles to grid nodes; and a hybrid impulse-based solver to enforce incompressibility on the grid. The efficacy of PFM has been demonstrated through various simulation scenarios, highlighting the evolution of complex vortical structures and the details of turbulent flows. Notably, compared to NFM, PFM reduces computing time by up to 49 times and memory consumption by up to 41%, while enhancing vorticity preservation as evidenced in various tests like leapfrog, vortex tube, and turbulent flow. 
    more » « less
  3. Abstract We performed the transport of a breast cancer cell (MB231-TGFb) in a microvessel using high-resolution simulations. Using open-source imaging software Slicer3D and Meshmixer, the 3D surface mesh forming the cell membrane was reconstructed from confocal microscopic images. The Dissipative Particle Dynamics method is used to model the cell membrane. The extracellular fluid flow is modeled with the Immersed Boundary Method to solve the governing equations of the blood plasma. The unsteady flow is applied at the inlet of the microchannel with an oscillatory pattern. Our results showed that the extracellular flow patterns are highly dependent on the waveform profile. The oscillatory flow showed the creation of vortices that influence the cellular deformations in the microchannel. These results could have implications on the destination of the cancer cells during transport in physiological flows. 
    more » « less
  4. We propose Coadjoint Orbit FLIP (CO-FLIP), a high order accurate, structure preserving fluid simulation method in the hybrid Eulerian-Lagrangian framework. We start with a Hamiltonian formulation of the incompressible Euler Equations, and then, using a local, explicit, and high order divergence free interpolation, construct a modified Hamiltonian system that governs our discrete Euler flow. The resulting discretization, when paired with a geometric time integration scheme, is energy and circulation preserving (formally the flow evolves on a coadjoint orbit) and is similar to the Fluid Implicit Particle (FLIP) method. CO-FLIP enjoys multiple additional properties including that the pressure projection is exact in the weak sense, and the particle-to-grid transfer is an exact inverse of the grid-to-particle interpolation. The method is demonstrated numerically with outstanding stability, energy, and Casimir preservation. We show that the method produces benchmarks and turbulent visual effects even at low grid resolutions. 
    more » « less
  5. Abstract This work presents the development of a novel approach to model cancer cell dynamics in microcirculation. The proposed numerical model is based on a hybrid continuum-particle approach. The cancer cell model includes the cell membrane, nucleus, cytoplasm and the cytoskeleton. The Dissipative Particle Dynamics method was employed to simulate the mechanical components. The blood plasma is modeled as a Newtonian incompressible fluid. A Fluid-Structure Interaction coupling, leveraging the Immersed Boundary Method is developed to simulate the cell's response to flow dynamics. The model is applied to resolve the transport of cancer cells with realistic morphologies in microcirculatory flows. Our results suggest that the controlling of oscillatory flows can be utilized to induce specific morphological shapes and the surrounding fluid patterns, allowing full manipulation and control of the cell. Furthermore, the intracellular and extracellular dynamics response of the cancer cell is intrinsically linked to their shape, in which certain morphologies displayed strong resistance to the fluid-induced forces and the ability to migrate in various directions. Our computational framework provides new capabilities for designing bioengineering devices for cell manipulation and separation. 
    more » « less