Effects of thermal noise on the transitional dynamics of an inextensible elastic filament in stagnation flow.

We investigate the dynamics of a single inextensible elastic filament subject to anisotropic friction in a viscous stagnation-point flow, by employing both a continuum model represented by Langevin type stochastic partial differential equations (SPDEs) and a dissipative particle dynamics (DPD) method. Unlike previous works, the filament is free to rotate and the tension along the filament is determined by the local inextensible constraint. The kinematics of the filament is recorded and studied with normal modes analysis.

Construction of dissipative particle dynamics models for complex fluids via the Mori-Zwanzig formulation.

We present a bottom-up coarse-graining procedure to construct mesoscopic force fields directly from microscopic dynamics. By grouping many bonded atoms in the molecular dynamics (MD) system into a single cluster, we compute both the conservative and non-conservative interactions between neighboring clusters. In particular, we perform MD simulations of polymer melts to provide microscopic trajectories for evaluating coarse-grained (CG) interactions.

Blood flow in small tubes: quantifying the transition to the non-continuum regime.

In small vessels blood is usually treated as a Newtonian fluid down to diameters of ~200 μm. We investigate the flow of red blood cell (RBC) suspensions driven through small tubes (diameters 10-150 μm) in the range marking the transition from arterioles and venules to the largest capillary vessels. The results of the simulations combined with previous simulations of uniform shear flow and experimental data show that for diameters less than ~100 μm the suspension's stress cannot be described as a continuum, even a heterogeneous one.

Simulation and modelling of slip flow over surfaces grafted with polymer brushes and glycocalyx fibres.

Fabrication of functionalized surfaces using polymer brushes is a relatively simple process and parallels the presence of glycocalyx filaments coating the luminal surface of our vasculature. In this paper, we perform atomistic-like simulations based on dissipative particle dynamics (DPD) to study both polymer brushes and glycocalyx filaments subject to shear flow, and we apply mean-field theory to extract useful scaling arguments on their response. For polymer brushes, a weak shear flow has no effect on the brush density profile or its height, while the slip length is independent of the shear rate and is of the order of the brush mesh size as a result of screening by hydrodynamic interactions.

Effect of chain chirality on the self-assembly of sickle hemoglobin.

We present simulation results on the self-assembly behavior of sickle hemoglobin (HbS). A coarse-grained HbS model, which contains hydrophilic and hydrophobic particles explicitly, is constructed to match the structural properties and physical description of HbS. The hydrophobic interactions are shown to be necessary with chirality being the main driver for the formation of HbS fibers.

Multiscale modeling of red blood cell mechanics and blood flow in malaria.

Red blood cells (RBCs) infected by a Plasmodium parasite in malaria may lose their membrane deformability with a relative membrane stiffening more than ten-fold in comparison with healthy RBCs leading to potential capillary occlusions. Moreover, infected RBCs are able to adhere to other healthy and parasitized cells and to the vascular endothelium resulting in a substantial disruption of normal blood circulation. In the present work, we simulate infected RBCs in malaria using a multiscale RBC model based on the dissipative particle dynamics method, coupling scales at the sub-cellular level with scales at the vessel size.

Predicting human blood viscosity in silico.

The viscosity of blood has long been used as an indicator in the understanding and treatment of disease, and the advent of modern viscometers allows its measurement with ever-improving clinical convenience. However, these advances have not been matched by theoretical developments that can yield a quantitative understanding of blood's microrheology and its possible connection to relevant biomolecules (e.g.

Predicting dynamics and rheology of blood flow: A comparative study of multiscale and low-dimensional models of red blood cells.

We compare the predictive capability of two mathematical models for red blood cells (RBCs) focusing on blood flow in capillaries and arterioles. Both RBC models as well as their corresponding blood flows are based on the dissipative particle dynamics (DPD) method, a coarse-grained molecular dynamics approach. The first model employs a multiscale description of the RBC (MS-RBC), with its membrane represented by hundreds or even thousands of DPD-particles connected by springs into a triangular network in combination with out-of-plane elastic bending resistance.

Wall shear stress-based model for adhesive dynamics of red blood cells in malaria.

Red blood cells (RBCs) infected by the Plasmodium falciparum (Pf-RBCs) parasite lose their membrane deformability and they also exhibit enhanced cytoadherence to vascular endothelium and other healthy and infected RBCs. The combined effect may lead to severe disruptions of normal blood circulation due to capillary occlusions. Here we extend the adhesion model to investigate the adhesive dynamics of Pf-RBCs as a function of wall shear stress (WSS) and other parameters using a three-dimensional, multiscale RBC model.

Blood flow and cell-free layer in microvessels.

Blood is modeled as a suspension of red blood cells using the dissipative particle dynamics method. The red blood cell membrane is coarse-grained for efficient simulations of multiple cells, yet accurately describes its viscoelastic properties. Blood flow in microtubes ranging from 10 to 40 μm in diameter is simulated in three dimensions for values of hematocrit in the range of 0.

A low-dimensional model for the red blood cell.

The red blood cell (RBC) is an important determinant of the rheological properties of blood because of its predominant number density, special mechanical properties and dynamics. Here, we develop a new low-dimensional RBC model based on dissipative particle dynamics (DPD). The model is constructed as a closed-torus-like ring of 10 colloidal particles connected by wormlike chain springs combined with bending resistance.

Systematic coarse-graining of spectrin-level red blood cell models.

We present a rigorous procedure to derive coarse-grained red blood cell (RBC) models, which yield accurate mechanical response. Based on a semi-analytic theory the linear and nonlinear elastic properties of healthy and infected RBCs in malaria can be matched with those obtained in optical tweezers stretching experiments. The present analysis predicts correctly the membrane Young's modulus in contrast to about 50% error in predictions by previous models.

A multiscale red blood cell model with accurate mechanics, rheology, and dynamics.

Red blood cells (RBCs) have highly deformable viscoelastic membranes exhibiting complex rheological response and rich hydrodynamic behavior governed by special elastic and bending properties and by the external/internal fluid and membrane viscosities. We present a multiscale RBC model that is able to predict RBC mechanics, rheology, and dynamics in agreement with experiments. Based on an analytic theory, the modeled membrane properties can be uniquely related to the experimentally established RBC macroscopic properties without any adjustment of parameters.

Steady shear rheometry of dissipative particle dynamics models of polymer fluids in reverse Poiseuille flow.

Polymer fluids are modeled with dissipative particle dynamics (DPD) as undiluted bead-spring chains and their solutions. The models are assessed by investigating their steady shear-rate properties. Non-Newtonian viscosity and normal stress coefficients, for shear rates from the lower to the upper Newtonian regimes, are calculated from both plane Couette and plane Poiseuille flows.

Direct construction of mesoscopic models from microscopic simulations.

Starting from microscopic molecular-dynamics (MD) simulations of constrained Lennard-Jones (LJ) clusters (with constant radius of gyration R(g)), we construct two mesoscopic models [Langevin dynamics and dissipative particle dynamics (DPD)] by coarse graining the LJ clusters into single particles. Both static and dynamic properties of the coarse-grained models are investigated and compared with the MD results. The effective mean force field is computed as a function of the intercluster distance, and the corresponding potential scales linearly with the number of particles per cluster and the temperature.

Rheology, microstructure and migration in brownian colloidal suspensions.

We demonstrate that suspended spherical colloidal particles can be effectively modeled as single dissipative particle dynamics (DPD) particles provided that the conservative repulsive force is appropriately chosen. The suspension model is further improved with a new formulation, which augments standard DPD with noncentral dissipative shear forces between particles while preserving angular momentum. Using the new DPD formulation we investigate the rheology, microstructure and shear-induced migration of a monodisperse suspension of colloidal particles in plane shear flows (Couette and Poiseuille).

Coarse-grained red blood cell model with accurate mechanical properties, rheology and dynamics.

We present a coarse-grained red blood cell (RBC) model with accurate and realistic mechanical properties, rheology and dynamics. The modeled membrane is represented by a triangular mesh which incorporates shear inplane energy, bending energy, and area and volume conservation constraints. The macroscopic membrane elastic properties are imposed through semi-analytic theory, and are matched with those obtained in optical tweezers stretching experiments.

Hydrodynamic interactions for single dissipative-particle-dynamics particles and their clusters and filaments.

We have investigated low Reynolds number flow past single dissipative-particle-dynamics particles (point centers of repulsion), their clusters, and their filaments using dissipative-particle-dynamics (DPD) simulations. The objective of our study was to verify whether DPD particles immersed in a sea of DPD particles behave like Langevin particles suspended in a continuous Newtonian fluid solvent, the basis of Brownian dynamics. Our principal test is to compare two effective DPD radii calculated by independent means.

Dissipative particle dynamics simulation of depletion layer and polymer migration in micro- and nanochannels for dilute polymer solutions.

The flows of dilute polymer solutions in micro- and nanoscale channels are of both fundamental and practical importance in variety of applications in which the channel gap is of the same order as the size of the suspended particles or macromolecules. In such systems depletion layers are observed near solid-fluid interfaces, even in equilibrium, and the imposition of flow results in further cross-stream migration of the particles. In this work we employ dissipative particle dynamics to study depletion and migration in dilute polymer solutions in channels several times larger than the radius of gyration (Rg) of bead-spring chains.

Schmidt number effects in dissipative particle dynamics simulation of polymers.

Simulation studies for dilute polymeric systems are presented using the dissipative particle dynamics method. By employing two different thermostats, the velocity-Verlet and Lowe's scheme, we show that the Schmidt number (S(c)) of the solvent strongly affects nonequilibrium polymeric quantities. The fractional extension of wormlike chains subjected to steady shear is obtained as a function of S(c).

Dissipative particle dynamics simulations of polymer chains: scaling laws and shearing response compared to DNA experiments.

Dissipative particle dynamics simulations of several bead-spring representations of polymer chains in dilute solution are used to demonstrate the correct static scaling laws for the radius of gyration. Shear flow results for the wormlike chain simulating single DNA molecules compare well with average extensions from experiments, irrespective of the number of beads. However, coarse graining with more than a few beads degrades the agreement of the autocorrelation of the extension.