Source author record

Michael Shelley

Michael Shelley appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.soft Biological Physics Computational Engineering, Finance, and Science cond-mat.mtrl-sci math.NA physics.comp-ph physics.flu-dyn Subcellular Processes

Catalog footprint

What is connected

5works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Comparison of explicit and mean-field models of cytoskeletal filaments with crosslinking motors

In cells, cytoskeletal filament networks are responsible for cell movement, growth, and division. Filaments in the cytoskeleton are driven and organized by crosslinking molecular motors. In reconstituted cytoskeletal systems, motor activity is responsible for far-from-equilibrium phenomena such as active stress, self-organized flow, and spontaneous nematic defect generation. How microscopic interactions between motors and filaments lead to larger-scale dynamics remains incompletely understood. To build from motor-filament interactions to predict bulk behavior of cytoskeletal systems, more computationally efficient techniques for modeling motor-filament interactions are needed. Here we derive a coarse-graining hierarchy of explicit and continuum models for crosslinking motors that bind to and walk on filament pairs. We compare the steady-state motor distribution and motor-induced filament motion for the different models and analyze their computational cost. All three models agree well in the limit of fast motor binding kinetics. Evolving a truncated moment expansion of motor density speeds the computation by $10^3$--$10^6$ compared to the explicit or continuous-density simulations, suggesting an approach for more efficient simulation of large networks. These tools facilitate further study of motor-filament networks on micrometer to millimeter length scales.

preprint2020arXiv

A scalable computational platform for particulate Stokes suspensions

We describe a computational framework for simulating suspensions of rigid particles in Newtonian Stokes flow. One central building block is a collision-resolution algorithm that overcomes the numerical constraints arising from particle collisions. This algorithm extends the well-known complementarity method for non-smooth multi-body dynamics to resolve collisions in dense rigid body suspensions. This approach formulates the collision resolution problem as a linear complementarity problem with geometric `non-overlapping' constraints imposed at each timestep. It is then reformulated as a constrained quadratic programming problem and the Barzilai-Borwein projected gradient descent method is applied for its solution. This framework is designed to be applicable for any convex particle shape, e.g., spheres and spherocylinders, and applicable to any Stokes mobility solver, including the Rotne-Prager-Yamakawa approximation, Stokesian Dynamics, and PDE solvers (e.g., boundary integral and immersed boundary methods). In particular, this method imposes Newton's Third Law and records the entire contact network. Further, we describe a fast, parallel, and spectrally-accurate boundary integral method tailored for spherical particles, capable of resolving lubrication effects. We show weak and strong parallel scalings up to $8\times 10^4$ particles with approximately $4\times 10^7$ degrees of freedom on $1792$ cores. We demonstrate the versatility of this framework with several examples, including sedimentation of particle clusters, and active matter systems composed of ensembles of particles driven to rotate.

preprint2016arXiv

A fast platform for simulating flexible fiber suspensions applied to cell mechanics

We present a novel platform for the large-scale simulation of fibrous structures immersed in a Stokesian fluid and evolving under confinement or in free-space. One of the main motivations for this work is to study the dynamics of fiber assemblies within biological cells. For this, we also incorporate the key biophysical elements that determine the dynamics of these assemblies, which include the polymerization and depolymerization kinetics of fibers, their interactions with molecular motors and other objects, their flexibility, and hydrodynamic coupling. This work, to our knowledge, is the first technique to include many-body hydrodynamic interactions (HIs), and the resulting fluid flows, in cellular fiber assemblies. We use the non-local slender body theory to compute the fluid-structure interactions of the fibers and a second-kind boundary integral formulation for other rigid bodies and the confining boundary. A kernel-independent implementation of the fast multiple method is utilized for efficient evaluation of HIs. The deformation of the fibers is described by the nonlinear Euler--Bernoulli beam theory and their polymerization is modeled by the reparametrization of the dynamic equations in the appropriate non-Lagrangian frame. We use a pseudo-spectral representation of fiber positions and implicit HIs in the time-stepping to resolve large fiber deformations, and to allow time-steps not constrained by temporal stiffness or fiber-fiber interactions. The entire computational scheme is parallelized, which enables simulating assemblies of thousands of fibers. We use our method to investigate two important questions in the mechanics of cell division: (i) the effect of confinement on the hydrodynamic mobility of microtubule asters; and (ii) the dynamics of the positioning of mitotic spindle in complex cell geometries.

preprint2015arXiv

Cytoplasmic flows as signatures for the mechanics of mitotic positioning

The proper positioning of the mitotic spindle is crucial for asymmetric cell division and generating cell diversity during development. Proper position in the single-cell embryo of Caenorhabditis elegans is achieved initially by the migration and rotation of the pronuclear complex (PNC) and its two associated centrosomal arrays of microtubules (MTs). We present here the first systematic theoretical study of how these $O(1000)$ centrosomal microtubules (MTs) interact through the immersing cytoplasm, the cell periphery and PNC, and with each other, to achieve proper position. This study is made possible through our development of a highly efficient and parallelized computational framework that accounts explicitly for long-ranged hydrodynamic interactions (HIs) between the MTs, while also capturing their flexibility, dynamic instability, and interactions with molecular motors and boundaries. First, we show through direct simulation that previous estimates of the PNC drag coefficient, based on either ignoring or partially including HIs, lead to misprediction of the active forces and time-scales of migration. We then directly study the dynamics of PNC migration under various force-transduction models, including the pushing or pulling of MTs at the cortex, and the pulling of MTs by cytoplasmically-bound force generators. While achieving proper position and orientation on physiologically reasonable time-scales does not uniquely choose a model, we find that each model produces a different signature in its induced cytoplasmic flow and MT conformations. We suggest then that cytoplasmic flows and MT conformations can be used to differentiate between mechanisms and to determine their contribution to the migration process.

preprint2010arXiv

Modeling and Simulation of Liquid Crystal Elastomers

We consider a continuum model describing the dynamic behavior of nematic liquid crystal elastomers (LCEs) and implement a numerical scheme to solve the governing equations. In the model, the Helmholtz free energy and Rayleigh dissipation are used, within a Lagrangian framework, to obtain the equations of motion. The free energy consists of both elastic and liquid crystalline contributions, each of which is a function of the material displacement and the orientational order parameter. The model gives dynamics for the material displacement, the scalar order parameter and the nematic director, the latter two of which correspond to the orientational order parameter tensor. Our simulations are carried out by solving the governing equations using an implicit-explicit scheme and the Chebyshev polynomial method. The simulations show that the model can successfully capture the shape changing dynamics of LCEs that have been observed in experiments, and also track the evolution of the order parameter tensor.