skip to main content


Title: A POROUS VISCOELASTIC MODEL FOR THE CELL CYTOSKELETON
The immersed boundary method is a widely used mixed Eulerian/Lagrangian framework for simulating the motion of elastic structures immersed in viscous fluids. In this work, we consider a poroelastic immersed boundary method in which a fluid permeates a porous, elastic structure of negligible volume fraction, and extend this method to include stress relaxation of the material. The porous viscoelastic method presented here is validated for a prescribed oscillatory shear and for an expansion driven by the motion at the boundary of a circular material by comparing numerical solutions to an analytical solution of the Maxwell model for viscoelasticity. Finally, an application of the modelling framework to cell biology is provided: passage of a cell through a microfluidic channel. We demonstrate that the rheology of the cell cytoplasm is important for capturing the transit time through a narrow channel in the presence of a pressure drop in the extracellular fluid.  more » « less
Award ID(s):
1664679
NSF-PAR ID:
10073319
Author(s) / Creator(s):
;
Date Published:
Journal Name:
The ANZIAM Journal
Volume:
59
Issue:
04
ISSN:
1446-1811
Page Range / eLocation ID:
472 to 498
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Microorganisms often navigate a complex environment composed of a viscous fluid with suspended microstructures such as elastic polymers and filamentous networks. These microstructures can have similar length scales to the microorganisms, leading to complex swimming dynamics. Some microorganisms secrete enzymes that dynamically change the elastic properties of the viscoelastic networks through which they move. In addition to biological organisms, microrobots have been engineered with the goals of mucin gel penetration or dissolving blood clots. In order to gain insight into the coupling between swimming performance and network remodeling, we used a regularized Stokeslet boundary element method to compute the motion of a microswimmer consisting of a rotating spherical body and counter-rotating helical flagellum. The viscoelastic network is represented by a network of points connected by virtual elastic linkages immersed in a viscous fluid. Here, we model the enzymatic dissolution of the network by bacteria or microrobots by dynamically breaking elastic linkages when the cell body of the swimmer falls within a given distance from the link. We investigate the swimming performance of the microbes as they penetrate and move through networks of different material properties, and also examine the effect of network remodeling. 
    more » « less
  2. Cell migration is critical for many vital processes, such as wound healing, as well as harmful processes, such as cancer metastasis. Experiments have highlighted the diversity in migration strategies employed by cells in physiologically relevant environments. In 3D fibrous matrices and confinement between two surfaces, some cells migrate using round membrane protrusions, called blebs. In bleb-based migration, the role of substrate adhesion is thought to be minimal, and it remains unclear if a cell can migrate without any adhesion complexes. We present a 2D computational fluid-structure model of a cell using cycles of bleb expansion and retraction in a channel with several geometries. The cell model consists of a plasma membrane, an underlying actin cortex, and viscous cytoplasm. Cellular structures are immersed in viscous fluid which permeates them, and the fluid equations are solved using the method of regularized Stokeslets. Simulations show that the cell cannot effectively migrate when the actin cortex is modeled as a purely elastic material. We find that cells do migrate in rigid channels if actin turnover is included with a viscoelastic description for the cortex. Our study highlights the non-trivial relationship between cell rheology and its external environment during migration with cytoplasmic streaming. 
    more » « less
  3. This study discusses turbulent suspension flows of non-Brownian, non-colloidal, neutrally buoyant and rigid spherical particles in a Newtonian fluid over porous media with particles too large to penetrate and move through the porous layer. We consider suspension flows with the solid volume fraction ${{\varPhi _b}}$ ranging from 0 to 0.2, and different wall permeabilities, while porosity is constant at 0.6. Direct numerical simulations with an immersed boundary method are employed to resolve the particles and flow phase, with the volume-averaged Navier–Stokes equations modelling the flow within the porous layer. The results show that in the presence of particles in the free-flow region, the mean velocity and the concentration profiles are altered with increasing porous layer permeability because of the variations in the slip velocity and wall-normal fluctuations at the suspension-porous interface. Furthermore, we show that variations in the stress condition at the interface significantly affect the particle near-wall dynamics and migration toward the channel core, thereby inducing large modulations of the overall flow drag. At the highest volume fraction investigated here, ${{\varPhi _b}}= 0.2$ , the velocity fluctuations and the Reynolds shear stress are found to decrease, and the overall drag increases due to the increase in the particle-induced stresses. 
    more » « less
  4. Although microtubules in plant cells have been extensively studied, the mechanisms that regulate the spatial organization of microtubules are poorly understood. We hypothesize that the interaction between microtubules and cytoplasmic flow plays an important role in the assembly and orientation of microtubules. To test this hypothesis, we developed a new computational modeling framework for microtubules based on theory and methods from the fluid-structure interaction. We employed the immersed boundary method to track the movement of microtubules in cytoplasmic flow. We also incorporated details of the encounter dynamics when two microtubules collide with each other. We verified our computational model through several numerical tests before applying it to the simulation of the microtubule-cytoplasm interaction in a growing plant cell. Our computational investigation demonstrated that microtubules are primarily oriented in the direction orthogonal to the axis of cell elongation. We validated the simulation results through a comparison with the measurement from laboratory experiments. We found that our computational model, with further calibration, was capable of generating microtubule orientation patterns that were qualitatively and quantitatively consistent with the experimental results. The computational model proposed in this study can be naturally extended to many other cellular systems that involve the interaction between microstructures and the intracellular fluid. 
    more » « less
  5. Abstract

    Although microtubules in plant cells have been extensively studied, the mechanisms that regulate the spatial organization of microtubules are poorly understood. We hypothesize that the interaction between microtubules and cytoplasmic flow plays an important role in the assembly and orientation of microtubules. To test this hypothesis, we developed a new computational modeling framework for microtubules based on theory and methods from the fluid–structure interaction. We employed the immersed boundary method to track the movement of microtubules in cytoplasmic flow. We also incorporated details of the encounter dynamics when two microtubules collide with each other. We verified our computational model through several numerical tests before applying it to the simulation of the microtubule–cytoplasm interaction in a growing plant cell. Our computational investigation demonstrated that microtubules are primarily oriented in the direction orthogonal to the axis of cell elongation. We validated the simulation results through a comparison with the measurement from laboratory experiments. We found that our computational model, with further calibration, was capable of generating microtubule orientation patterns that were qualitatively and quantitatively consistent with the experimental results. The computational model proposed in this study can be naturally extended to many other cellular systems that involve the interaction between microstructures and the intracellular fluid.

     
    more » « less