Effect of microchannel structure and fluid properties on non-inertial particle migration.

In this work, we investigate the influence of channel structure and fluid rheology on non-inertial migration of non-Brownian polystyrene beads. Particle migration in this regime can be found in biomedical, chemical, environmental and geological applications. However, the effect of fluid rheology on particle migration in porous media remains to be clearly understood. Here, we isolate the effects of elasticity and shear thinning by comparing a Newtonian fluid, a purely elastic (Boger) fluid, and a shear-thinning elastic fluid. To mimic the complexity of geometries in real-world application, a random porous structure is created through a disordered arrangement of cylindrical pillars in the microchannel. Experiments are repeated in an empty channel and in channels with an ordered arrangement of pillars, and the similarities and differences in the observed particle focusing are analyzed. It is found that elasticity drives the particles away from the channel walls in an empty microchannel. Notably, particle focusing is unaffected by curved streamlines in an ordered porous microchannel and particles stay away from pillars in elastic fluids. Shear-thinning is found to reduce the effect of focusing and a broader region of particle concentration is observed. It is also noteworthy that the rheological characteristics of the fluid are not important for the particle distribution in a randomly arranged pillared microchannel and particles have a uniform distribution for all suspending fluids. Moreover, discussion on the current discrepancy in the literature about the equilibrium positions of the particles in a channel is extended by analyzing the results obtained in the current experiments.

Viscoelastic flow past mono- and bidisperse random arrays of cylinders: flow resistance, topology and normal stress distribution.

We investigate creeping viscoelastic fluid flow through two-dimensional porous media consisting of random arrangements of monodisperse and bidisperse cylinders, using our finite volume-immersed boundary method introduced in S. De, et al., J. Non-Newtonian Fluid Mech., 2016, 232, 67-76. The viscoelastic fluid is modeled with a FENE-P model. The simulations show an increased flow resistance with increase in flow rate, even though the bulk response of the fluid to shear flow is shear thinning. We show that if the square root of the permeability is chosen as the characteristic length scale in the determination of the dimensionless Deborah number (De), then all flow resistance curves collapse to a single master curve, irrespective of the pore geometry. Our study reveals how viscoelastic stresses and flow topologies (rotation, shear and extension) are distributed through the porous media, and how they evolve with increasing De. We correlate the local viscoelastic first normal stress differences with the local flow topology and show that the largest normal stress differences are located in shear flow dominated regions and not in extensional flow dominated regions at higher viscoelasticity. The study shows that normal stress differences in shear flow regions may play a crucial role in the increase of flow resistance for viscoelastic flow through such porous media.

Momentum conserving Brownian dynamics propagator for complex soft matter fluids.

We present a Galilean invariant, momentum conserving first order Brownian dynamics scheme for coarse-grained simulations of highly frictional soft matter systems. Friction forces are taken to be with respect to moving background material. The motion of the background material is described by locally averaged velocities in the neighborhood of the dissolved coarse coordinates. The velocity variables are updated by a momentum conserving scheme. The properties of the stochastic updates are derived through the Chapman-Kolmogorov and Fokker-Planck equations for the evolution of the probability distribution of coarse-grained position and velocity variables, by requiring the equilibrium distribution to be a stationary solution. We test our new scheme on concentrated star polymer solutions and find that the transverse current and velocity time auto-correlation functions behave as expected from hydrodynamics. In particular, the velocity auto-correlation functions display a long time tail in complete agreement with hydrodynamics.

Coarse-grained simulations of moderately entangled star polyethylene melts.

In this paper, a previous coarse-grain model [J. T. Padding and W. J. Briels, J. Chem. Phys. 117, 925 (2002)] to simulate melts of linear polymers has been adapted to simulate polymers with more complex hierarchies. Bond crossings between highly coarse-grained soft particles are prevented by applying an entanglement algorithm. We first test our method on a virtual branch point inside a linear chain to make sure it works effectively when linking two linear arms. Next, we apply our method to study the diffusive and rheological behaviors of a melt of three-armed stars. We find that the diffusive behavior of the three-armed star is very close to that of a linear polymer with the same molecular weight, while its rheological properties are close to those of a linear chain with molecular mass equal to that of the longest linear sub-chain in the star.

How Péclet number affects microstructure and transient cluster aggregation in sedimenting colloidal suspensions.

We study how varying the Péclet number (Pe) affects the steady state sedimentation of colloidal particles that interact through short-ranged attractions. By employing a hybrid molecular dynamics simulation method we demonstrate that the average sedimentation velocity changes from a non-monotonic dependence on packing fraction φ at low Pe numbers, to a monotonic decrease with φ at higher Pe numbers. At low Pe number the pair correlation functions are close to their equilibrium values, but as the Pe number increases, important deviations from equilibrium forms are observed. Although the attractive forces we employ are not strong enough to form permanent clusters, they do induce transient clusters whose behaviour is also affected by Pe number. In particular, clusters are more likely to fragment and less likely to aggregate at larger Pe numbers, and the probability of finding larger clusters decreases with increasing Pe number. Interestingly, the lifetime of the clusters is more or less independent of Pe number in the range we study. Instead, the change in cluster distribution occurs because larger clusters are less likely to form with increasing Pe number. These results illustrate some of the subtleties that occur in the crossover from equilibrium like to purely non-equilibrium behaviour as the balance between convective and thermal forces changes.

Diffusion of spherical particles in microcavities.

The diffusive motion of a colloidal particle trapped inside a small cavity filled with fluid is reduced by hydrodynamic interactions with the confining walls. In this work, we study these wall effects on a spherical particle entrapped in a closed cylinder. We calculate the diffusion coefficient along the radial, azimuthal, and axial direction for different particle positions. At all locations the diffusion is smaller than in a bulk fluid and it becomes anisotropic near the container's walls. We present a simple model which reasonably well describes the simulation results for the given dimensions of the cylinder, which are taken from the recent experimental work.

Alignment of particles in sheared viscoelastic fluids.

We investigate the shear-induced structure formation of colloidal particles dissolved in non-Newtonian fluids by means of computer simulations. The two investigated visco-elastic fluids are a semi-dilute polymer solution and a worm-like micellar solution. Both shear-thinning fluids contain long flexible chains whose entanglements appear and disappear continually as a result of Brownian motion and the applied shear flow. To reach sufficiently large time and length scales in three-dimensional simulations with up to 96 spherical colloids, we employ the responsive particle dynamics simulation method of modeling each chain as a single soft Brownian particle with slowly evolving inter-particle degrees of freedom accounting for the entanglements. Parameters in the model are chosen such that the simulated rheological properties of the fluids, i.e., the storage and loss moduli and the shear viscosities, are in reasonable agreement with experimental values. Spherical colloids dispersed in both quiescent fluids mix homogeneously. Under shear flow, however, the colloids in the micellar solution align to form strings in the flow direction, whereas the colloids in the polymer solution remain randomly distributed. These observations agree with recent experimental studies of colloids in the bulk of these two liquids.

Effects of interparticle attractions on colloidal sedimentation.

We use a mesoscopic simulation technique to study the effect of short-ranged interparticle attractions on the steady-state sedimentation of colloidal suspensions. Attractions increase the average sedimentation velocity v(s) compared to the pure hard-sphere case, and for strong enough attractions, a nonmonotonic dependence on the packing fraction phi with a maximum velocity at intermediate phi is observed. Attractions also strongly enhance hydrodynamic velocity fluctuations, which show a pronounced maximum size as a function of phi. These phenomena arise from a complex interplay between nonequilibrium hydrodynamic effects and the thermodynamics of transient cluster formation.

Translational and rotational friction on a colloidal rod near a wall.

We present particulate simulation results for translational and rotational friction components of a shish-kebab model of a colloidal rod with aspect ratio (length over diameter) L/D=10 in the presence of a planar hard wall. Hydrodynamic interactions between rod and wall cause an overall enhancement of the friction tensor components. We find that the friction enhancements to reasonable approximation scale inversely linear with the closest distance d between the rod surface and the wall, for d in the range between D/8 and L. The dependence of the wall-induced friction on the angle theta between the long axis of the rod and the normal to the wall is studied and fitted with simple polynomials in cos theta.

Chain dynamics of poly(ethylene-alt-propylene) melts by means of coarse-grained simulations based on atomistic molecular dynamics.

We present coarse-grained molecular dynamics simulations of poly(ethylene-alt-propylene) (PEP) melts, ranging in chain length from about N(e) (the entanglement length) to N=6N(e). The coarse-grained parameters, potential of mean force and bare friction, were determined from fully atomistic molecular dynamics simulations carried out on a PEP cell containing 12 chains of 80 monomers each and subjected to periodic boundary conditions. These atomistic simulations were previously validated by means of extensive neutron scattering measurements. Uncrossability constrains were also introduced in the coarse-grained model to prevent unphysical bond crossing. The coarse-grained simulations were carried out at 492 K and focus on chain dynamics. The results obtained were analyzed in terms of Rouse coordinates and Rouse correlators. We observe deviations from Rouse behavior for all chain lengths investigated, even when the chain stiffness is incorporated in the Rouse model. These deviations become more important as the chain length increases. The general scenario emerging from the results obtained is that the deviations from Rouse-like behavior are due to correlations among the forces acting upon a chain bead, which seem to be related with the constraint of uncrossability among the chains. As consequence, nonexponentiality of the Rouse correlators and mode- and time-dependent friction are observed. It seems that, in the molecular weight explored, these effects still give not raise to reptation behavior but to a crossover regime between Rouse and reptation. On the other hand, the results obtained are in qualitative agreement with those expected from the so-called generalized Rouse models, based on memory function formalisms.

Efficient simulation of noncrossing fibers and chains in a hydrodynamic solvent.

An efficient simulation method is presented for Brownian fiber suspensions, which includes both uncrossability of the fibers and hydrodynamic interactions between the fibers mediated by a mesoscopic solvent. To conserve hydrodynamics, collisions between the fibers are treated such that momentum and energy are conserved locally. The choice of simulation parameters is rationalized on the basis of dimensionless numbers expressing the relative strength of different physical processes. The method is applied to suspensions of semiflexible fibers with a contour length equal to the persistence length, and a mesh size to contour length ratio ranging from 0.055 to 0.32. For such fibers the effects of hydrodynamic interactions are observable, but relatively small. The noncrossing constraint, on the other hand, is very important and leads to hindered displacements of the fibers, with an effective tube diameter in agreement with recent theoretical predictions. The simulation technique opens the way to study the effect of viscous effects and hydrodynamic interactions in microrheology experiments where the response of an actively driven probe bead in a fiber suspension is measured.

Flow around fishlike shapes studied using multiparticle collision dynamics.

Empirical measurements of hydrodynamics of swimming fish are very difficult. Therefore, modeling studies may be of great benefit. Here, we investigate the suitability for such a study of a recently developed mesoscale method, namely, multiparticle collision dynamics. As a first step, we confine ourselves to investigations at intermediate Reynolds numbers of objects that are stiff. Due to the lack of empirical data on the hydrodynamics of stiff fishlike shapes we use a previously published numerical simulation of the shapes of a fish and a tadpole for comparison. Because the shape of a tadpole resembles that of a circle with an attached splitter plate, we exploit the knowledge on hydrodynamic consequences of such an attachment to test the model further and study the effects of splitter plates for objects of several shapes at several Reynolds numbers. Further, we measure the angles of separation of flow around a circular cylinder and make small adjustments to the boundary condition and the method to drive the flow. Our results correspond with empirical data and with results from other models.

Dynamics and rheology of wormlike micelles emerging from particulate computer simulations.

We perform coarse-grained computer simulations of solutions of semidilute wormlike micelles and study their dynamic and rheological properties, both in equilibrium and under shear flow. The simulation model is tailored to the study of relatively large time and length scales (micrometers and several milliseconds), while it still retains the specific mechanical properties of the individual wormlike micelles. The majority of the mechanical properties (persistence length, diameter, and elastic modulus of a single worm) is determined from more detailed atomistic molecular dynamics simulations, providing the link with the chemistry of the surfactants. The method is applied to the case of a solution containing 8% (by weight) erucyl bis(hydroxymethyl)methylammonium chloride (EHAC). Different scission energies ranging from 15.5k(b)T to 19.1k(B)T are studied, leading to both unentangled and entangled wormlike micelles. We find a decrease in the average contour length and an increase in the average breaking rate with increasing shear rate. In equilibrium, the decay of the shear relaxation modulus of the unentangled samples agrees with predictions based on a theory of breakable Rouse chains. Under shear flow, transient over- and undershoots are measured in the stress tensor components. At high shear rates we observe a steady-state shear stress proportional to gamma(1/3), where gamma is the shear rate. This is confirmed by our high shear rate experiments of real EHAC in a parallel-plate geometry.

Spinodal decomposition of asymmetric binary fluids in a micro-Couette geometry simulated with molecular dynamics.

The spinodal decomposition of quenched polymer/solvent and liquid-crystal/solvent mixtures in a miniature Taylor-Couette cell has been simulated by molecular dynamics. Three stacking motifs, each reflecting the geometry and symmetry of the cell, are most abundant among the fully phase separated stationary states. At zero or low angular velocity of the inner cylindrical drum, the two segregated domains have a clear preference for the stacking with the lowest free energy and hence the smallest total interfacial tension. For high shear rates, the steady state appears to be determined by a minimum dissipation mechanism, i.e., the mixtures are likely to evolve into the stacking demanding the least mechanical power by the rotating wall. The partial slip at the polymer-solvent interfaces then gives rise to a new pattern: A stack of three concentric cylindrical shells with the viscous polymer layer sandwiched between two solvent layers. Neither of these mechanisms can explain all simulation results, as the separating mixture easily becomes kinetically trapped in a long-lived suboptimal configuration. The phase separation process is observed to proceed faster under shear than in a quiescent mixture.

Finite system size effects in the interfacial dynamics of binary liquid films.

We study the relaxation dynamics of capillary waves in the interface between two confined liquid layers by means of molecular dynamics simulations. We measure the autocorrelations of the interfacial Fourier modes and find that the finite thickness of the liquid layers leads to a marked increase of the relaxation times as compared to the case of fluid layers of infinite depth. The simulation results are in good agreement with a theoretical first-order perturbation derivation, which starts from the overdamped Stokes' equation. The theory also takes into account an interfacial friction, but the difference with no-slip interfacial conditions is small. When the walls are sheared, it is found that the relaxation times of modes perpendicular to the flow are unaffected. Modes along the flow direction are relatively unaffected as long as the equilibrium relaxation time is sufficiently short compared to the rate of deformation. We discuss the consequences for experiments on thin layers and on ultralow surface tension fluids, as well as computer simulations.

Interplay between hydrodynamic and Brownian fluctuations in sedimenting colloidal suspensions.

We apply a hybrid molecular dynamics and mesoscopic simulation technique to study the steady-state sedimentation of hard sphere particles for Peclet number (Pe) ranging from 0.08 to 12. Hydrodynamic backflow causes a reduction of the average sedimentation velocity relative to the Stokes velocity. We find that this effect is independent of Pe number. Velocity fluctuations show the expected effects of thermal fluctuations at short correlation times. At longer times, nonequilibrium hydrodynamic fluctuations are visible, and their character appears to be independent of the thermal fluctuations. The hydrodynamic fluctuations dominate the diffusive behavior even for modest Pe number, while conversely the short-time fluctuations are dominated by thermal effects for surprisingly large Pe numbers. Inspired by recent experiments, we also study finite sedimentation in a horizontal planar slit. In our simulations distinct lateral patterns emerge, in agreement with observations in the experiments.

Viscoelastic effects on residual oil distribution in flows through pillared microchannels.

HYPOTHESIS: Multiphase flow through porous media is important in a number of industrial, natural and biological processes. One application is enhanced oil recovery (EOR), where a resident oil phase is displaced by a Newtonian or polymeric fluid. In EOR, the two-phase immiscible displacement through heterogonous porous media is usually governed by competing viscous and capillary forces, expressed through a Capillary number Ca, and viscosity ratio of the displacing and displaced fluid. However, when viscoelastic displacement fluids are used, elastic forces in the displacement fluid also become significant. It is hypothesized that elastic instabilities are responsible for enhanced oil recovery through an elastic microsweep mechanism. EXPERIMENTS: In this work, we use a simplified geometry in the form of a pillared microchannel. We analyze the trapped residual oil size distribution after displacement by a Newtonian fluid, a nearly inelastic shear thinning fluid, and viscoelastic polymers and surfactant solutions. FINDINGS: We find that viscoelastic polymers and surfactant solutions can displace more oil compared to Newtonian fluids and nearly inelastic shear thinning polymers at similar Ca numbers. Beyond a critical Ca number, the size of residual oil blobs decreases significantly for viscoelastic fluids. This critical Ca number directly corresponds to flow rates where elastic instabilities occur in single phase flow, suggesting a close link between enhancement of oil recovery and appearance of elastic instabilities.

Particle-based modeling of heterogeneous chemical kinetics including mass transfer.

Connecting the macroscopic world of continuous fields to the microscopic world of discrete molecular events is important for understanding several phenomena occurring at physical boundaries of systems. An important example is heterogeneous catalysis, where reactions take place at active surfaces, but the effective reaction rates are determined by transport limitations in the bulk fluid and reaction limitations on the catalyst surface. In this work we study the macro-micro connection in a model heterogeneous catalytic reactor by means of stochastic rotation dynamics. The model is able to resolve the convective and diffusive interplay between participating species, while including adsorption, desorption, and reaction processes on the catalytic surface. Here we apply the simulation methodology to a simple straight microchannel with a catalytic strip. Dimensionless Damkohler numbers are used to comment on the spatial concentration profiles of reactants and products near the catalyst strip and in the bulk. We end the discussion with an outlook on more complicated geometries and increasingly complex reactions.

Hydrodynamic interactions and Brownian forces in colloidal suspensions: coarse-graining over time and length scales.

We describe in detail how to implement a coarse-grained hybrid molecular dynamics and stochastic rotation dynamics simulation technique that captures the combined effects of Brownian and hydrodynamic forces in colloidal suspensions. The importance of carefully tuning the simulation parameters to correctly resolve the multiple time and length scales of this problem is emphasized. We systematically analyze how our coarse-graining scheme resolves dimensionless hydrodynamic numbers such as the Reynolds number Re, which indicates the importance of inertial effects, the Schmidt number Sc, which indicates whether momentum transport is liquidlike or gaslike, the Mach number, which measures compressibility effects, the Knudsen number, which describes the importance of noncontinuum molecular effects, and the Peclet number, which describes the relative effects of convective and diffusive transport. With these dimensionless numbers in the correct regime the many Brownian and hydrodynamic time scales can be telescoped together to maximize computational efficiency while still correctly resolving the physically relevant processes. We also show how to control a number of numerical artifacts, such as finite-size effects and solvent-induced attractive depletion interactions. When all these considerations are properly taken into account, the measured colloidal velocity autocorrelation functions and related self-diffusion and friction coefficients compare quantitatively with theoretical calculations. By contrast, these calculations demonstrate that, notwithstanding its seductive simplicity, the basic Langevin equation does a remarkably poor job of capturing the decay rate of the velocity autocorrelation function in the colloidal regime, strongly underestimating it at short times and strongly overestimating it at long times. Finally, we discuss in detail how to map the parameters of our method onto physical systems and from this extract more general lessons-keeping in mind that there is no such thing as a free lunch-that may be relevant for other coarse-graining schemes such as lattice Boltzmann or dissipative particle dynamics.

Mechanical properties of surfactant bilayer membranes from atomistic and coarse-grained molecular dynamics simulations.

We use simulations to predict the stability and mechanical properties of two amphiphilic bilayer membranes. We carry out atomistic MD simulations and investigate whether it is possible to use an existing coarse-grained (CG) surfactant model to map the membrane properties. We find that certain membranes can be represented well by the CG model, whereas others cannot. Atomistic MD simulations of the erucate membrane yield a headgroup area per surfactant a(0) of 0.26 nm(2), an elastic modulus K(A) of 1.7 N/m, and a bending rigidity kappa of 5 k(B)T. We find that the CG model, with the right choice for the size and potential well depth of the head, correctly reproduces a(0), kappa, as well as the fluctuation spectrum over the whole range of q values. Atomistic MD simulations of EHAC, on the other hand, suggest that this membrane is unstable. This is indicated by the fact that kappa is of the order of k(B)T, which means that the interface is extremely flexible and diffuse, and K(A) is close to zero, which means that the surface tension is zero. We argue that the CG model can be used if the headgroups are uncharged, dipolar, or effectively dipolar due to headgroup charge screening induced by counterion condensation.