Fast computation of many-particle hydrodynamic and electrostatic interactions in a confined geometry

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a "dry diffusion" approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

A fast spectral method for electrostatics in doubly periodic slit channels

We develop a fast method for computing the electrostatic energy and forces for a collection of charges in doubly periodic slabs with jumps in the dielectric permittivity at the slab boundaries. Our method achieves spectral accuracy by using Ewald splitting to replace the original Poisson equation for nearly singular sources with a smooth far-field Poisson equation, combined with a localized near-field correction. Unlike existing spectral Ewald methods, which make use of the Fourier transform in the aperiodic direction, we recast the problem as a two-point boundary value problem in the aperiodic direction for each transverse Fourier mode for which exact analytic boundary conditions are available. We solve each of these boundary value problems using a fast, well-conditioned Chebyshev method. In the presence of dielectric jumps, combining Ewald splitting with the classical method of images results in smoothed charge distributions, which overlap the dielectric boundaries themselves. We show how to preserve the spectral accuracy in this case through the use of a harmonic correction, which involves solving a simple Laplace equation with smooth boundary data. We implement our method on graphical processing units and combine our doubly periodic Poisson solver with Brownian dynamics to study the equilibrium structure of double layers in binary electrolytes confined by dielectric boundaries. Consistent with prior studies, we find strong charge depletion near the interfaces due to repulsive interactions with image charges, which points to the need for incorporating polarization effects in understanding confined electrolytes, both theoretically and computationally.

Coil-stretch-like transition of elastic sheets in extensional flows

The conformation of a long linear polymer dissolved in fluid and exposed to an extensional flow is well-known to exhibit a "coil-stretch" transition, which for sufficiently long chains can lead to bistability. The present work reports computations indicating that an analogous "compact-stretched" transition arises in the dynamics of a thin elastic sheet. Sheets of nominally circular, square or rectangular shape are simulated in planar and biaxial flows using a finite element method for the sheet conformations and a regularized Stokeslet method for the fluid flow. If a neo-Hookean constitutive model is used for the sheet elasticity, the sheets will stretch without bound once a critical extension rate, as characterized nondimensionally by a capillary number, is exceeded. Nonlinear elasticity, represented with the Yeoh model, arrests the stretching, leading to a highly-stretched steady state once the critical capillary number is exceeded. For all shapes and in both planar and biaxial extension, a parameter regime exists in which both weakly stretched (compact) and strongly stretched states can be found, depending on initial conditions. I.e. this parameter regime displays bistability. As in the long-chain polymer case, the bistable behavior arises from the hydrodynamic interaction between distant elements of the sheet, and vanishes if these interactions are artificially screened by use of a Brinkman model for the fluid motion. While the sheets can transiently display wrinkled shapes, all final shapes in planar and biaxial extension are planar.

A simple microswimmer model inspired by the general equation for nonequilibrium reversible-irreversible coupling

A simple mean-field microswimmer model is presented. The model is inspired by the nonequilibrium thermodynamics of multi-component fluids that undergo chemical reactions. These thermodynamics can be rigorously described in the context of the GENERIC (general equation for the nonequilibrium reversible-irreversible coupling) framework. More specifically, this approach was recently applied to non-ideal polymer solutions [T. Indei and J. D. Schieber, J. Chem. Phys. 146, 184902 (2017)]. One of the species of the solution is an unreactive polymer chain represented by the bead-spring model. Using this detailed description as inspiration, we then make several simplifying assumptions to obtain a mean-field model for a Janus microswimmer. The swimmer model considered here consists of a polymer dumbbell in a sea of reactants. One of the beads of the dumbbell is allowed to act as a catalyst for a chemical reaction between the reactants. We show that the mean-squared displacement (MSD) of the center of mass of this Janus dumbbell exhibits ballistic behavior at time scales at which the concentration of the reactant is large. The time scales at which the ballistic behavior is observed in the MSD coincide with the time scales at which the cross-correlation between the swimmer's orientation and the direction of its displacement exhibits a maximum. Since the swimmer model was inspired by the GENERIC framework, it is possible to ensure that the entropy generation is always positive, and therefore, the second law of thermodynamics is obeyed.

Flow-induced segregation and dynamics of red blood cells in sickle cell disease

Blood flow in sickle cell disease (SCD) can substantially differ from normal blood flow due to significant alterations in the physical properties of the red blood cells (RBCs). Chronic complications, such as inflammation of the endothelial cells lining blood vessel walls, are associated with SCD, for reasons that are unclear. Here, detailed boundary integral simulations are performed to investigate an idealized model flow flow in SCD, a binary suspension of flexible biconcave discoidal fluid-filled capsules and stiff curved prolate capsules that represent healthy and sickle RBCs, respectively, subjected to pressure-driven flow in a planar slit. The stiff component is dilute. The key observation is that, unlike healthy RBCs that concentrate around the center of the channel and form an RBC-depleted layer (i.e. cell-free layer) next to the walls, sickle cells are largely drained from the bulk of the suspension and aggregate inside the cell-free layer, displaying strong margination. These cells are found to undergo a rigid-body-like rolling orbit near the walls. A binary suspension of flexible biconcave discoidal capsules and stiff straight (non-curved) prolate capsules is also considered for comparison, and the curvature of the stiff component is found to play a minor role in the behavior. Additionally, by considering a mixture of flexible and stiff biconcave discoids, we reveal that rigidity difference by itself is sufficient to induce the segregation behavior in a binary suspension. Furthermore, the additional shear stress on the walls induced by the presence of cells is computed for the various cases. Compared to the small fluctuations in wall shear stress for a suspension of healthy RBCs, large local peaks in wall shear stress are observed for the binary suspensions, due to the proximity of the marginated stiff cells to the walls. This effect is most marked for the straight prolate capsules. As endothelial cells are known to mechanotransduce physical forces such as aberrations in shear stress and convert them to physiological processes such as activation of inflammatory signals, these results may aid in understanding mechanisms for endothelial dysfunction associated with SCD.

Hydrodynamic effects on the motility of crawling eukaryotic cells

Eukaryotic cell motility is crucial during development, wound healing, the immune response, and cancer metastasis. Some eukaryotic cells can swim, but cells more commonly adhere to and crawl along the extracellular matrix. We study the relationship between hydrodynamics and adhesion that describe whether a cell is swimming, crawling, or combining these motions. Our simple model of a cell, based on the three-sphere swimmer, is capable of both swimming and crawling. As cell-matrix adhesion strength increases, the influence of hydrodynamics on migration diminishes. Cells with significant adhesion can crawl with speeds much larger than their nonadherent, swimming counterparts. We predict that, while most eukaryotic cells are in the strong-adhesion limit, increasing environment viscosity or decreasing cell-matrix adhesion could lead to significant hydrodynamic effects even in crawling cells. Signatures of hydrodynamic effects include a dependence of cell speed on the presence of a nearby substrate or interactions between noncontacting cells. These signatures will be suppressed at large adhesion strengths, but even strongly adherent cells will generate relevant fluid flows that will advect nearby passive particles and swimmers.

Temporal evolution of concentration and microstructure of colloidal films during vertical drying: a lattice Boltzmann simulation study

We study the temporal and structural development of colloid films during vertical drying using the lattice Boltzmann (LB) simulation. The dispersed particles moving in Brownian motion have excluded volume and hydrodynamic interactions in the film. The concentrated colloidal film formed by solvent evaporation is modeled as an uniaxial compression of colloids with a planar moving interface. The simulation studies are carried out over a wide range of Péclet number (Pe), the relative ratio between the evaporation rate and the diffusion rate of colloids. The results clearly demonstrate a temporal variation of colloid concentration as the evaporation rate increases. In the case of high Pe, the increase of colloid concentration in the top layer creates structural features that can be distinguished along the height of the film, and eventually can induce a large tensile stress in the layer. However, surprisingly, the colloids are maximally crystallized in the case of moderate Pe. The LB simulation results are further compared with those from previous studies of the Brownian Dynamics (BD) simulation and the continuum model for the evaporation film. The LB and BD results match well both at low and high Pe limits. The qualitatively significant differences between LB and BD simulations at a moderate Pe indicate that hydrodynamic interactions (HIs) play an important role in this Pe. The presence of HIs induces a greater reduction of diffusion than under geometrical restriction alone, and the effect is conspicuous when particles are driven both by diffusion and by advection.

Structure and proton conduction in sulfonated poly(ether ether ketone) semi-permeable membranes: a multi-scale computational approach

The design of polymeric membranes for proton or ionic exchange highly depends on the fundamental understanding of the physical and molecular mechanisms that control the formation of the conduction channels. There is an inherent relation between the dynamical structure of the polymeric membrane and the electrostatic forces that drive membrane segregation and proton transport. Here, we used a multi-scale computational approach to analyze the morphology of sulfonated poly(ether ether ketone) membranes at the mesoscale. A self-consistent description of the electrostatic phenomenon was adopted, where discrete polymer chains and a continuum proton field were embedded in a continuum fluid. Brownian dynamics was used for the evolution of the suspended polymer molecules, while a convection-diffusion transport equation, including the Nernst-Planck diffusion mechanism, accounted for the dynamics of the proton concentration field. We varied the polymer concentration, the degree of sulfonation and the level of confinement to find relationships between membrane structure and proton conduction. Our results indicate that the reduced mobility of polymer chains, at concentrations above overlap, and a moderate degree of sulfonation - i.e., 30% - are essential elements for membrane segregation and proton domain connectivity. These conditions also ensure that the membrane structure is not affected by size or by potential gradients. Importantly, our analysis shows that membrane conductivity and current are linearly dependent on polymer concentration and quadratically dependent on the degree of sulfonation. We found that the optimal polymeric membrane design requires a polymer concentration above overlap and a degree of sulfonation around 50%. These conditions promote a dynamical membrane morphology with a constant density of proton channels. Our results and measurements agree with previous experimental works, thereby validating our model and observations.

Discrete and Continuum Models for the Salt in Crowded Environments of Suspended Charged Particles

Electrostatic forces greatly affect the overall dynamics and diffusional activities of suspended charged particles in crowded environments. Accordingly, the concentration of counter- or co-ions in a fluid-''the salt"-determines the range, strength, and order of electrostatic interactions between particles. This environment fosters engineering routes for controlling directed assembly of particles at both the micro- and nanoscale. Here, we analyzed two computational modeling schemes that considered salt within suspensions of charged particles, or polyelectrolytes: discrete and continuum. Electrostatic interactions were included through a Green's function formalism, where the confined fundamental solution for Poisson's equation is resolved by the general geometry Ewald-like method. For the discrete model, the salt was considered as regularized point-charges with a specific valence and size, while concentration fields were defined for each ionic species for the continuum model. These considerations were evolved using Brownian dynamics of the suspended charged particles and the discrete salt ions, while a convection-diffusion transport equation, including the Nernst-Planck diffusion mechanism, accounted for the dynamics of the concentration fields. The salt/particle models were considered as suspensions under slit-confinement conditions for creating crowded "macro-ions", where density distributions and radial distribution functions were used to compare and differentiate computational models. Importantly, our analysis shows that disparate length scales or increased system size presented by the salt and suspended particles are best dealt with using concentration fields to model the ions. These findings were then validated by novel simulations of a semipermeable polyelectrolyte membrane, at the mesoscale, from which ionic channels emerged and enable ion conduction.

Brownian dynamics of confined suspensions of active microrollers

We develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M. Youssef, S. Sacanna, A. Donev, and P. Chaikin, Nat. Phys. (2016), preprint arXiv:1609.08673. We present a stochastic Adams-Bashforth integrator for the equations of Brownian dynamics, which has the same cost but is more accurate than the widely used Euler-Maruyama scheme, and use a random finite difference to capture the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. We generate the Brownian increments using a Krylov method and show that for particles confined to remain in the vicinity of a no-slip wall by gravity or active flows, the number of iterations is independent of the number of particles. Our numerical experiments with active rollers show that the thermal fluctuations set the characteristic height of the colloids above the wall, both in the initial condition and the subsequent evolution dominated by active flows. The characteristic height in turn controls the time scale and wavelength for the development of the fingering instability.

Multi-resolution polymer Brownian dynamics with hydrodynamic interactions

A polymer model given in terms of beads, interacting through Hookean springs and hydrodynamic forces, is studied. A Brownian dynamics description of this bead-spring polymer model is extended to multiple resolutions. Using this multiscale approach, a modeller can efficiently look at different regions of the polymer in different spatial and temporal resolutions with scalings given for the number of beads, statistical segment length, and bead radius in order to maintain macro-scale properties of the polymer filament. The Boltzmann distribution of a Gaussian chain for differing statistical segment lengths gives a diffusive displacement equation for the multi-resolution model with a mobility tensor for different bead sizes. Using the pre-averaging approximation, the translational diffusion coefficient is obtained as a function of the inverse of a matrix and then in closed form in the long-chain limit. This is then confirmed with numerical experiments.

Numerical simulation of particulate flows using a hybrid of finite difference and boundary integral methods

A combination of finite difference (FD) and boundary integral (BI) methods is used to formulate an efficient solver for simulating unsteady Stokes flow around particles. The two-dimensional (2D) unsteady Stokes equation is being solved on a Cartesian grid using a second order FD method, while the 2D steady Stokes equation is being solved near the particle using BI method. The two methods are coupled within the viscous boundary layer, a few FD grid cells away from the particle, where solutions from both FD and BI methods are valid. We demonstrate that this hybrid method can be used to accurately solve for the flow around particles with irregular shapes, even though radius of curvature of the particle surface is not resolved by the FD grid. For dilute particle concentrations, we construct a virtual envelope around each particle and solve the BI problem for the flow field located between the envelope and the particle. The BI solver provides velocity boundary condition to the FD solver at "boundary" nodes located on the FD grid, adjacent to the particles, while the FD solver provides the velocity boundary condition to the BI solver at points located on the envelope. The coupling between FD method and BI method is implicit at every time step. This method allows us to formulate an O(N) scheme for dilute suspensions, where N is the number of particles. For semidilute suspensions, where particles may cluster, an envelope formation method has been formulated and implemented, which enables solving the BI problem for each individual particle cluster, allowing efficient simulation of hydrodynamic interaction between particles even when they are in close proximity. The method has been validated against analytical results for flow around a periodic array of cylinders and for Jeffrey orbit of a moving ellipse in shear flow. Simulation of multiple force-free irregular shaped particles in the presence of shear in a 2D slit flow has been conducted to demonstrate the robustness of the FD-BI scheme. This method can also be used to solve the full Navier-Stokes equations around particles. In this case, it is shown that, below a threshold particle Reynolds number, which in turn depends on the required resolution of the particle surface, the FD-BI scheme will be more computationally efficient than a pure FD scheme solved on a multiblock grid.

Modeling the relaxation of internal DNA segments during genome mapping in nanochannels

We have developed a multi-scale model describing the dynamics of internal segments of DNA in nanochannels used for genome mapping. In addition to the channel geometry, the model takes as its inputs the DNA properties in free solution (persistence length, effective width, molecular weight, and segmental hydrodynamic radius) and buffer properties (temperature and viscosity). Using pruned-enriched Rosenbluth simulations of a discrete wormlike chain model with circa 10 base pair resolution and a numerical solution for the hydrodynamic interactions in confinement, we convert these experimentally available inputs into the necessary parameters for a one-dimensional, Rouse-like model of the confined chain. The resulting coarse-grained model resolves the DNA at a length scale of approximately 6 kilobase pairs in the absence of any global hairpin folds, and is readily studied using a normal-mode analysis or Brownian dynamics simulations. The Rouse-like model successfully reproduces both the trends and order of magnitude of the relaxation time of the distance between labeled segments of DNA obtained in experiments. The model also provides insights that are not readily accessible from experiments, such as the role of the molecular weight of the DNA and location of the labeled segments that impact the statistical models used to construct genome maps from data acquired in nanochannels. The multi-scale approach used here, while focused towards a technologically relevant scenario, is readily adapted to other channel sizes and polymers.

Computational immuno-biology for organ transplantation and regenerative medicine

Organ transplantation and regenerative medicine are adopted platforms that provide replacement tissues and organs from natural or engineered sources. Acceptance, tolerance and rejection depend greatly on the proper control of the immune response against graft antigens, motivating the development of immunological and genetical therapies that prevent organ failure. They rely on a complete, or partial, understanding of the immune system. Ultimately, they are innovative technologies that ensure permanent graft tolerance and indefinite graft survival through the modulation of the immune system. Computational immunology has arisen as a tool towards a mechanistic understanding of the biological and physicochemical processes surrounding an immune response. It comprehends theoretical and computational frameworks that simulate immuno-biological systems. The challenge is centered on the multi-scale character of the immune system that spans from atomistic scales, during peptide-epitope and protein interactions, to macroscopic scales, for lymph transport and organ-organ reactions. In this paper, we discuss, from an engineering perspective, the biological processes that are involved during the immune response of organ transplantation. Previous computational efforts, including their characteristics and visible limitations, are described. Finally, future perspectives and challenges are listed to motivate further developments.

Shape-mediated margination and demargination in flowing multicomponent suspensions of deformable capsules

We present detailed simulations and theory for flow-induced segregation in suspensions of deformable fluid-filled capsules with different shapes during simple shear flow in a planar slit. This system is an idealized model for transport for blood cells and/or drug carriers in the microcirculation or in microfluidic devices. For the simulations, an accelerated implementation of the boundary integral method was employed. We studied the binary mixtures of spherical and ellipsoidal capsules, varying the aspect ratio κ of the ellipsoid while keeping constant either (a) equatorial radius or (b) volume. Effects of a variety of parameters was studied, including κ, volume fraction and number fraction of the spherical capsules in the mixture. In suspensions where the ellipsoids have the same equatorial radius as the spheres, capsules with lower κ marginate. In suspension where the ellipsoids have the same volume as the spheres, ellipsoidal (both oblate and prolate) capsules are seen to demarginate in a mixture of primarily spherical capsules. To understand these results, a mechanistic framework based on the competition between wall-induced migration and shear-induced collisions is presented. A simplified drift-diffusion theory based on this framework shows excellent qualitative agreement with simulation results.

Getting drowned in a swirl: Deformable bead-spring model microswimmers in external flow fields

Deformability is a central feature of many types of microswimmers, e.g., for artificially generated self-propelled droplets. Here, we analyze deformable bead-spring microswimmers in an externally imposed solvent flow field as simple theoretical model systems. We focus on their behavior in a circular swirl flow in two spatial dimensions. Linear (straight) two-bead swimmers are found to circle around the swirl with a slight drift to the outside with increasing activity. In contrast to that, we observe for triangular three-bead or squarelike four-bead swimmers a tendency of being drawn into the swirl and finally getting drowned, although a radial inward component is absent in the flow field. During one cycle around the swirl, the self-propulsion direction of an active triangular or squarelike swimmer remains almost constant, while their orbits become deformed exhibiting an "egglike" shape. Over time, the swirl flow induces slight net rotations of these swimmer types, which leads to net rotations of the egg-shaped orbits. Interestingly, in certain cases, the orbital rotation changes sense when the swimmer approaches the flow singularity. Our predictions can be verified in real-space experiments on artificial microswimmers.

Effects of confinement on models of intracellular macromolecular dynamics

The motions of particles in a viscous fluid confined within a spherical cell have been simulated using Brownian and Stokesian dynamics simulations. High volume fractions mimicking the crowded interior of biological cells were used. Importantly, although confinement yields an overall slowdown in motion, the qualitative effects of motion in the interior of the cell can be effectively modeled as if the system were an infinite periodic system. However, we observe layering of particles at the cell wall due to steric interactions in the confined space. Motions of nearby particles are also strongly correlated at the cell wall, and these correlations increase when hydrodynamic interactions are modeled. Further, particles near the cell wall have a tendency to remain near the cell wall. A consequence of these effects is that the mean contact time between particles is longer at the cell wall than in the interior of the cell. These findings identify a specific way that confinement affects the interactions between particles and points to a previously unidentified mechanism that may play a role in signal transduction and other processes near the membrane of biological cells.

Dynamics of a single red blood cell in simple shear flow

This work describes simulations of a red blood cell (RBC) in simple shear flow, focusing on the dependence of the cell dynamics on the spontaneous curvature of the membrane. The results show that an oblate spheroidal spontaneous curvature maintains the dimple of the RBC during tank-treading dynamics as well as exhibits off-shear-plane tumbling consistent with the experimental observations of Dupire et al. [J. Dupire, M. Socol, and A. Viallat, Proc. Natl. Acad. Sci. USA 109, 20808 (2012)] and their hypothesis of an inhomogeneous spontaneous shape. As the flow strength (capillary number Ca) is increased at a particular viscosity ratio between inner and outer fluid, the dynamics undergo transitions in the following sequence: tumbling, kayaking or rolling, tilted tank-treading, oscillating-swinging, swinging, and tank-treading. The tilted tank-treading (or spinning frisbee) regime has been previously observed in experiments but not in simulations. Two distinct classes of regime are identified: a membrane reorientation regime, where the part of membrane that is at the dimple at rest moves to the rim and vice versa, is observed in motions at high Ca such as tilted tank-treading, oscillating-swinging, swinging, and tank-treading, and a nonreorientation regime, where the part of the membrane starting from the dimple stays at the dimple, is observed in motions at low Ca such as rolling, tumbling, kayaking, and flip-flopping.

Matrix-free Brownian dynamics simulation technique for semidilute polymeric solutions

Evaluating the concentration dependence of static and dynamic properties of macromolecules in semidilute polymer solutions requires accurate calculation of long-range hydrodynamic interactions (HI) and short range excluded volume (EV) forces. In conventional Brownian dynamics simulations (BDS), computation of HI necessitates construction of a dense diffusion tensor commonly performed via Ewald summation. Krylov subspace techniques allow efficient decomposition of this tensor [computational cost scales as O(N^{2}), where N is the total number of beads in bead-spring representation of macromolecules in a simulation box] and computation of Brownian displacements in the box. In this paper, a matrix-free approach for calculation of HI is implemented which leads to O(NlogN) scaling of computational expense. The fidelity of the algorithm is demonstrated by evaluating the asymptotic value of center-of-mass diffusivity of polymer molecules at very low concentrations and their radius of gyration scaling as a function of number of beads for dilute and semidilute solutions (with concentrations up to 5 times the overlap concentration). In turn, a favorable comparison between our results and the blob theory is shown.

Self-consistent description of electrokinetic phenomena in particle-based simulations

A new computational method is presented for study suspensions of charged particles undergoing fluctuating hydrodynamic and electrostatic interactions. The proposed model is appropriate for polymers, proteins, and porous particles embedded in a continuum electrolyte. A self-consistent Langevin description of the particles is adopted in which hydrodynamic and electrostatic interactions are included through a Green's function formalism. An Ewald-like split is adopted in order to satisfy arbitrary boundary conditions for the Stokeslet and Poisson Green functions, thereby providing a formalism that is applicable to any geometry and that can be extended to deformable objects. The convection-diffusion equation for the continuum ions is solved simultaneously considering Nernst-Planck diffusion. The method can be applied to systems at equilibrium and far from equilibrium. Its applicability is demonstrated in the context of electrokinetic motion, where it is shown that the ionic clouds associated with individual particles can be severely altered by the flow and concentration, leading to intriguing cooperative effects.

Motion of an elastic capsule in a constricted microchannel

We study the motion of an elastic capsule through a microchannel characterized by a localized constriction. We consider a capsule with a stress-free spherical shape and impose its steady-state configuration in an infinitely long straight channel as the initial condition for our calculations. We report how the capsule deformation, velocity, retention time, and maximum stress of the membrane are affected by the capillary number, Ca , and the constriction shape. We estimate the deformation by measuring the variation of the three-dimensional surface area and a series of alternative quantities easier to extract from experiments. These are the Taylor parameter, the perimeter and the area of the capsule in the spanwise plane. We find that the perimeter is the quantity that best reproduces the behavior of the three-dimensional surface area. This is maximum at the centre of the constriction and shows a second peak after it, whose location depends on the Ca number. We observe that, in general, area-deformation-correlated quantities grow linearly with Ca , while velocity-correlated quantities saturate for large Ca but display a steeper increase for small Ca . The velocity of the capsule divided by the velocity of the flow displays, surprisingly, two different qualitative behaviors for small and large capillary numbers. Finally, we report that longer constrictions and spanwise wall bounded (versus spanwise periodic) domains cause larger deformations and velocities. If the deformation and velocity in the spanwise wall bounded domains are rescaled by the initial equilibrium deformation and velocity, their behavior is undistinguishable from that in a periodic domain. In contrast, a remarkably different behavior is reported in sinusoidally shaped and smoothed rectangular constrictions indicating that the capsule dynamics is particularly sensitive to abrupt changes in the cross section. In a smoothed rectangular constriction larger deformations and velocities occur over a larger distance.

Hydrodynamics of DNA confined in nanoslits and nanochannels

Modeling the dynamics of a confined, semi exible polymer is a challenging problem, owing to the complicated interplay between the configurations of the chain, which are strongly affected by the length scale for the confinement relative to the persistence length of the chain, and the polymer-wall hydrodynamic interactions. At the same time, understanding these dynamics are crucial to the advancement of emerging genomic technologies that use confinement to stretch out DNA and "read" a genomic signature. In this mini-review, we begin by considering what is known experimentally and theoretically about the friction of a wormlike chain such as DNA confined in a slit or a channel. We then discuss how to estimate the friction coefficient of such a chain, either with dynamic simulations or via Monte Carlo sampling and the Kirk-wood pre-averaging approximation. We then review our recent work on computing the diffusivity of DNA in nanoslits and nanochannels, and conclude with some promising avenues for future work and caveats about our approach.

Molecular-scale simulation of cross-flow migration in polymer melts

The first ever molecular-scale simulation of cross-flow migration effects in dense polymer melts is presented; simulations for both unentangled and entangled chains are presented. At quiescence a small depletion next to the wall for the segmental densities of longer chains is present, a corresponding excess exists about one-half a radii of gyration away from the wall, and uniform values are observed further from the wall. In shear flow the melts exhibit similar behavior as the quiescent case; a constant shear rate across the gap does not induce chain length based migration. In contradistinction, parabolic flow (where gradients in shear rate are present) causes profound migration for both unentangled and entangled melts. Mapping onto polyethylene and calculating stress shows the system is far below the stress required to break chains. Accordingly, our findings are consistent with flow induced migration mechanisms predominating over competing chain degradation mechanisms thus resolving a 40 year old controversy.

A microfluidic device to sort capsules by deformability: a numerical study

Guided by extensive numerical simulations, we propose a microfluidic device that can sort elastic capsules by their deformability. The device consists of a duct embedded with a semi-cylindrical obstacle, and a diffuser which further enhances the sorting capability. We demonstrate that the device can operate reasonably well under changes in the initial position of the capsule. The efficiency of the device remains essentially unaltered under small changes of the obstacle shape (from semi-circular to semi-elliptic cross-section). Confinement along the direction perpendicular to the plane of the device increases its efficiency. This work is the first numerical study of cell sorting by a realistic microfluidic device.

A minimally-resolved immersed boundary model for reaction-diffusion problems

We develop an immersed boundary approach to modeling reaction-diffusion processes in dispersions of reactive spherical particles, from the diffusion-limited to the reaction-limited setting. We represent each reactive particle with a minimally-resolved "blob" using many fewer degrees of freedom per particle than standard discretization approaches. More complicated or more highly resolved particle shapes can be built out of a collection of reactive blobs. We demonstrate numerically that the blob model can provide an accurate representation at low to moderate packing densities of the reactive particles, at a cost not much larger than solving a Poisson equation in the same domain. Unlike multipole expansion methods, our method does not require analytically computed Green's functions, but rather, computes regularized discrete Green's functions on the fly by using a standard grid-based discretization of the Poisson equation. This allows for great flexibility in implementing different boundary conditions, coupling to fluid flow or thermal transport, and the inclusion of other effects such as temporal evolution and even nonlinearities. We develop multigrid-based preconditioners for solving the linear systems that arise when using implicit temporal discretizations or studying steady states. In the diffusion-limited case the resulting linear system is a saddle-point problem, the efficient solution of which remains a challenge for suspensions of many particles. We validate our method by comparing to published results on reaction-diffusion in ordered and disordered suspensions of reactive spheres.

Dynamic simulation of concentrated macromolecular solutions with screened long-range hydrodynamic interactions: algorithm and limitations

Hydrodynamic interactions exert a critical effect on the dynamics of macromolecules. As the concentration of macromolecules increases, by analogy to the behavior of semidilute polymer solutions or the flow in porous media, one might expect hydrodynamic screening to occur. Hydrodynamic screening would have implications both for the understanding of macromolecular dynamics as well as practical implications for the simulation of concentrated macromolecular solutions, e.g., in cells. Stokesian dynamics (SD) is one of the most accurate methods for simulating the motions of N particles suspended in a viscous fluid at low Reynolds number, in that it considers both far-field and near-field hydrodynamic interactions. This algorithm traditionally involves an O(N(3)) operation to compute Brownian forces at each time step, although asymptotically faster but more complex SD methods are now available. Motivated by the idea of hydrodynamic screening, the far-field part of the hydrodynamic matrix in SD may be approximated by a diagonal matrix, which is equivalent to assuming that long range hydrodynamic interactions are completely screened. This approximation allows sparse matrix methods to be used, which can reduce the apparent computational scaling to O(N). Previously there were several simulation studies using this approximation for monodisperse suspensions. Here, we employ newly designed preconditioned iterative methods for both the computation of Brownian forces and the solution of linear systems, and consider the validity of this approximation in polydisperse suspensions. We evaluate the accuracy of the diagonal approximation method using an intracellular-like suspension. The diffusivities of particles obtained with this approximation are close to those with the original method. However, this approximation underestimates intermolecular correlated motions, which is a trade-off between accuracy and computing efficiency. The new method makes it possible to perform large-scale and long-time simulation with an approximate accounting of hydrodynamic interactions.

Presentation of large DNA molecules for analysis as nanoconfined dumbbells

The analysis of very large DNA molecules intrinsically supports long-range, phased sequence information, but requires new approaches for their effective presentation as part of any genome analysis platform. Using a multi-pronged approach that marshaled molecular confinement, ionic environment, and DNA elastic properties-but tressed by molecular simulations-we have developed an efficient and scalable approach for presentation of large DNA molecules within nanoscale slits. Our approach relies on the formation of DNA dumbbells, where large segments of the molecules remain outside the nanoslits used to confine them. The low ionic environment, synergizing other features of our approach, enables DNA molecules to adopt a fully stretched conformation, comparable to the contour length, thereby facilitating analysis by optical microscopy. Accordingly, a molecular model is proposed to describe the conformation and dynamics of the DNA molecules within the nanoslits; a Langevin description of the polymer dynamics is adopted in which hydrodynamic effects are included through a Green's function formalism. Our simulations reveal that a delicate balance between electrostatic and hydrodynamic interactions is responsible for the observed molecular conformations. We demonstrate and further confirm that the "Odijk regime" does indeed start when the confinement dimensions size are of the same order of magnitude as the persistence length of the molecule. We also summarize current theories concerning dumbbell dynamics.

Dynamics of solutes with hydrodynamic interactions: comparison between Brownian dynamics and stochastic rotation dynamics simulations

The dynamics of particles in solution or suspension is influenced by thermal fluctuations and hydrodynamic interactions. Several mesoscale methods exist to account for these solvent-induced effects such as Brownian dynamics with hydrodynamic interactions and hybrid molecular dynamics-stochastic rotation dynamics methods. Here we compare two ways of coupling solutes to the solvent with stochastic rotation dynamics (SRD) to Brownian dynamics with and without explicit hydrodynamic interactions. In the first SRD scheme [SRD with collisional coupling (CC)] the solutes participate in the collisional step with the solvent and in the second scheme [SRD with central force coupling (CFC)] the solutes interact through direct forces with the solvent, generating slip boundary conditions. We compare the transport coefficients of neutral and charged solutes in a model system obtained by these simulation schemes. Brownian dynamics without hydrodynamic interactions is used as a reference to quantify the influence of hydrodynamics on the transport coefficients as modeled by the different methods. We show that, in the dilute range, the SRD CFC method provides results similar to those of Brownian dynamics with hydrodynamic interactions for the diffusion coefficients and for the electrical conductivity. The SRD CC scheme predicts diffusion coefficients close to those obtained by Brownian dynamic simulations without hydrodynamic interactions, but accounts for part of the influence of hydrodynamics on the electrical conductivity.

Krylov subspace methods for computing hydrodynamic interactions in brownian dynamics simulations

Hydrodynamic interactions play an important role in the dynamics of macromolecules. The most common way to take into account hydrodynamic effects in molecular simulations is in the context of a brownian dynamics simulation. However, the calculation of correlated brownian noise vectors in these simulations is computationally very demanding and alternative methods are desirable. This paper studies methods based on Krylov subspaces for computing brownian noise vectors. These methods are related to Chebyshev polynomial approximations, but do not require eigenvalue estimates. We show that only low accuracy is required in the brownian noise vectors to accurately compute values of dynamic and static properties of polymer and monodisperse suspension models. With this level of accuracy, the computational time of Krylov subspace methods scales very nearly as O(N(2)) for the number of particles N up to 10 000, which was the limit tested. The performance of the Krylov subspace methods, especially the "block" version, is slightly better than that of the Chebyshev method, even without taking into account the additional cost of eigenvalue estimates required by the latter. Furthermore, at N = 10,000, the Krylov subspace method is 13 times faster than the exact Cholesky method. Thus, Krylov subspace methods are recommended for performing large-scale brownian dynamics simulations with hydrodynamic interactions.

Optimization of a Brownian-dynamics algorithm for semidilute polymer solutions

Simulating the static and dynamic properties of semidilute polymer solutions with Brownian dynamics (BD) requires the computation of a large system of polymer chains coupled to one another through excluded-volume and hydrodynamic interactions. In the presence of periodic boundary conditions, long-ranged hydrodynamic interactions are frequently summed with the Ewald summation technique. By performing detailed simulations that shed light on the influence of several tuning parameters involved both in the Ewald summation method, and in the efficient treatment of Brownian forces, we develop a BD algorithm in which the computational cost scales as O(N(1.8)), where N is the number of monomers in the simulation box. We show that Beenakker's original implementation of the Ewald sum, which is only valid for systems without bead overlap, can be modified so that θ solutions can be simulated by switching off excluded-volume interactions. A comparison of the predictions of the radius of gyration, the end-to-end vector, and the self-diffusion coefficient by BD, at a range of concentrations, with the hybrid lattice Boltzmann-molecular dynamics (LB-MD) method shows excellent agreement between the two methods. In contrast to the situation for dilute solutions, the LB-MD method is shown to be significantly more computationally efficient than the current implementation of BD for simulating semidilute solutions. We argue, however, that further optimizations should be possible.

Segregation by membrane rigidity in flowing binary suspensions of elastic capsules

Spatial segregation in the wall normal direction is investigated in suspensions containing a binary mixture of neo-Hookean capsules subjected to pressure driven flow in a planar slit. The two components of the binary mixture have unequal membrane rigidities. The problem is studied numerically using an accelerated implementation of the boundary integral method. The effect of a variety of parameters was investigated, including the capillary number, rigidity ratio between the two species, volume fraction, confinement ratio, and number fraction of the more floppy particle Xf in the mixture. It was observed that in suspensions of pure species, the mean wall normal positions of the stiff and the floppy particles are comparable. In mixtures, however, the stiff particles were found to be increasingly displaced toward the walls with increasing Xf, while the floppy particles were found to increasingly accumulate near the centerline with decreasing Xf. This segregation behavior was universally observed independent of the parameters. The origin of this segregation is traced to the effect of the number fraction Xf on the localization of the stiff and the floppy particles in the near wall region--the probability of escape of a stiff particle from the near wall region to the interior is greatly reduced with increasing Xf, while the exact opposite trend is observed for a floppy particle with decreasing Xf. Simple model studies on heterogeneous pair collisions involving a stiff and a floppy particle mechanistically explain the contrasting effect of Xf on the near wall localization of the two species. The key observation in these studies is that the stiff particle experiences much larger cross-stream displacement in heterogeneous collisions than the floppy particle. A unified mechanism incorporating the wall-induced migration of deformable particles away from the wall and the particle fluxes associated with heterogeneous and homogeneous pair collisions is presented.

Coexistence of tight and loose bundled states in a model of bacterial flagellar dynamics

Many microorganisms propel themselves through their fluid environment by means of multiple rotating flagella that self-organize to form bundles, a process that is complex and poorly understood. In the present work, the bundling behavior of a pair of flexible flagella, each driven by a constant torque motor, is investigated, using a mathematical model incorporating the fluid motion generated by each flagellum as well as the finite flexibility of the flagella. The initial stage of bundling is driven purely by hydrodynamics but the final state of the bundle is determined by a nontrivial balance between hydrodynamics and elasticity. As the flexibility of the flagella increases a regime is found where, depending on initial conditions, one finds bundles that are either tight, with the flagella in mechanical contact, or loose, with the flagella intertwined but not touching. That is, multiple coexisting states of bundling are found. The parameter regime (in terms of flexibility and distance between motors) at which this multiplicity occurs is comparable to the parameters for a number of bacteria.

Designing microcapsule arrays that propagate chemical signals

Using analysis and simulation, we show how ordered arrays of microcapsules in solution can be harnessed to propagate chemical signals in directed and controllable ways, allowing the signals to be transmitted over macroscopic distances. The system encompasses two types of capsules that are localized on an adhesive surface. The "signaling" capsules release inducer molecules, which trigger "targets" to release nanoparticles. The released nanoparticles can bind to the underlying surface and thus, create adhesion gradients, which then propel the signaling capsules to shuttle between neighboring targets. This arrangement acts like a relay, so that triggering target capsules at a particular location in the array also triggers target capsules in adjacent locations. For an array containing two target columns, our simulations and analysis show that steady input signal leads to a sustained periodic output. For an array containing multiple target columns, we show that by introducing a prescribed ratio of nanoparticle release rates between successive target columns, a chemical signal can be propagated along the array without dissipation. We also demonstrate that similar signal transmission cannot be performed via diffusion alone.