Two-particle friction in a mesoscopic solvent

The effects of hydrodynamic interactions on the friction tensors for two particles in solution are studied. The particles have linear dimensions on nanometer scales and are either simple spherical particles interacting with the solvent through repulsive Lennard-Jones forces or are composite cluster particles whose atomic components interact with the solvent through repulsive Lennard-Jones forces. The solvent dynamics is modeled at a mesoscopic level through multiparticle collisions that conserve mass, momentum, and energy. The dependence of the two-particle relative friction tensors on the interparticle separation indicates the importance of hydrodynamic interactions for these nanoparticles.     

Creeping motion of an assemblage of composite spheres relative to a fluid

The sedimentation of a homogeneous distribution of spherical composite particles and the fluid flow through a bed of these particles are investigated theoretically. Each composite particle is composed of a spherical solid core and a surrounding porous shell. In the fluid-permeable porous shell, idealized hydrodynamic frictional segments are assumed to distribute uniformly. The effect of interactions among the particles is taken into explicit account by employing a fundamental cell-model representation which is known to provide good predictions for the motion of a swarm of nonporous spheres within a fluid. In the limit of a small Reynolds number, the Stokes and Brinkman equations are solved for the flow field in a unit cell, and the drag force exerted by the fluid on the particle is obtained in a closed form. For a distribution of composite spheres, the normalized mobility of the particles decreases or the particle interactions increase monotonically with a decrease in the permeability of their porous shells. The effect of particle interactions on the creeping motion of composite spheres relative to a fluid can be quite significant in some situations. In the limiting cases, the analytical solutions describing the drag force or mobility for a suspension of composite spheres reduce to those for suspensions of solid spheres and of porous spheres. The hydrodynamic behavior for composite spheres may be approximated by that for permeable spheres when the porous layer is sufficiently thick, depending on the permeability.     

Short-time transport properties in dense suspensions: from neutral to charge-stabilized colloidal spheres

We present a detailed study of short-time dynamic properties in concentrated suspensions of charge-stabilized and of neutral colloidal spheres. The particles in many of these systems are subject to significant many-body hydrodynamic interactions. A recently developed accelerated Stokesian dynamics (ASD) simulation method is used to calculate hydrodynamic functions, wave-number-dependent collective diffusion coefficients, self-diffusion and sedimentation coefficients, and high-frequency limiting viscosities. The dynamic properties are discussed in dependence on the particle concentration and salt content. Our ASD simulation results are compared with existing theoretical predictions, notably those of the renormalized density fluctuation expansion method of Beenakker and Mazur [Physica A 126, 349 (1984)], and earlier simulation data on hard spheres. The range of applicability and the accuracy of various theoretical expressions for short-time properties are explored through comparison with the simulation data. We analyze, in particular, the validity of generalized Stokes-Einstein relations relating short-time diffusion properties to the high-frequency limiting viscosity, and we point to the distinctly different behavior of de-ionized charge-stabilized systems in comparison to hard spheres.     

Mesoscopic model for diffusion-influenced reaction dynamics

A hybrid mesoscopic multiparticle collision model is used to study diffusion-influenced reaction kinetics. The mesoscopic particle dynamics conserves mass, momentum, and energy so that hydrodynamic effects are fully taken into account. Reactive and nonreactive interactions with catalytic solute particles are described by full molecular dynamics. Results are presented for large-scale, three-dimensional simulations to study the influence of diffusion on the rate constants of the A + C <==> B + C reaction. In the limit of a dilute solution of catalytic C particles, the simulation results are compared with diffusion equation approaches for both the irreversible and reversible reaction cases. Simulation results for systems where the volume fraction phi of catalytic spheres is high are also presented, and collective interactions among reactions on catalytic spheres that introduce volume fraction dependence in the rate constants are studied.     

Long-time self-diffusion of charged colloidal particles: electrokinetic and hydrodynamic interaction effects

The authors analyze the long-time self-diffusion of charge-stabilized colloidal macroions in nondilute suspensions using a mode-coupling scheme developed for multicomponent suspensions of interacting Brownian spheres. In this scheme, all ionic species, including counterions and electrolyte ions, are treated on an equal footing as charged hard spheres undergoing overdamped Brownian motion. Hydrodynamic interactions between all ions are accounted for on the far-field level. We show that the influence on the colloidal long-time self-diffusion coefficient arising from the relaxation of the microionic atmosphere surrounding the colloids, the so-called electrolyte friction effect, is usually insignificant in comparison with the friction contributions arising from direct and hydrodynamic interactions between the colloidal particles. This finding is true even for small colloid concentrations unless the mobility difference between colloidal particles and microions is not large. Furthermore, we observe an interesting nonmonotonic density dependence of the colloidal long-time self-diffusion coefficient in suspensions with low amount of added salt. We show that this unusual density dependence is due to colloid-colloid hydrodynamic interactions.     

Velocity domain and volume fraction distribution of heavy microparticles in low reynolds number flow in microchannel

In this paper, a two-dimensional Stokes flow having particle size distribution ranging from 10 to 50 µm in a rectangular channel is simulated numerically with focus on the hydrodynamic forces. The results show that due to the disparity between the density of the fluid and particles, velocity domain of particles deviates from the fluid velocity domain and this phenomenon occurs significantly for the larger particles. Also, with the increase of Reynolds number, the volume fraction of dispersed phase near the bottom wall of the channel increases either. Compared to similar studies, this investigation employs numerical simulation of microparticulate flow and interparticle hydrodynamic forces with emphasis on the dispersed phase volume fraction in order to present the microchannel flow properties.     

Magnetophoresis, sedimentation, and diffusion of particles in concentrated magnetic fluids

A dynamic mass transfer equation for describing magnetophoresis, sedimentation, and gradient diffusion of colloidal particles in concentrated magnetic fluids has been derived. This equation takes into account steric, magnetodipole, and hydrodynamic interparticle interactions. Steric interactions have been investigated using the Carnahan-Starling approximation for a hard-sphere system. In order to study the effective interparticle attraction, the free energy of the dipolar hard-sphere system is represented as a virial expansion with accuracy to the terms quadratic in particle concentration. The virial expansion gives an interpolation formula that fits well the results of computer simulation in a wide range of particle concentrations and interparticle interaction energies. The diffusion coefficient of colloidal particles is written with regard to steric, magnetodipole and hydrodynamic interactions. We thereby laid the foundation for the formulation of boundary-value problems and for calculation of concentration and magnetic fields in the devices (for example, magnetic fluid seals and acceleration sensors), which use a concentrated magnetic fluid as a working fluid. The Monte-Carlo methods and the analytical approach are employed to study the magnetic fluid stratification generated by the gravitational field in a cylinder of finite height. The coefficient of concentration stratification of the magnetic fluid is calculated in relation to the average concentration of particles and the dipolar coupling constant. It is shown that the effective particle attraction causes a many-fold increase in the concentration inhomogeneity of the fluid if the average volume fraction of particles does not exceed 30%. At high volume concentrations steric interactions play a crucial role.     

Theory of relaxation spectra of polymer networks with interchain friction and long-range hydrodynamic interactions

The dynamic network model taken into account the friction relative to incompressible effective viscous medium (EVM), the hydrodynamic interactions of a network with the EVM, and the interchain friction is considered. Two collective network relaxation spectra (RS) arise. The first RS corresponds to the network motion relative to immobile EVM, the EVM does not take part in this motion. The second RS includes the combined viscous motion of the network and EVM because of long-range hydrodynamic interactions. This RS is finite and narrow for infinitely large network. The symmetry of this motion is in concordance with that of incompressible EVM. The existence of interchain friction leads to the narrowing of both types of relaxation spectra.     

Hydrodynamic interactions between spheres in a viscous fluid with a flat free surface or hard wall

Hydrodynamic interactions between spheres immersed in a low-Reynolds-number fluid flow close to a flat free surface or hard wall are investigated. The spheres may have different or equal radii, and may be separated from the boundary or at contact with the free surface. A simple and useful expression is derived for the propagator (Green operator) connecting centers of two spheres. In the derivation, the method of images and the displacement theorems are used. Symmetry of the displacement operators is explicitly shown. The significance of these results in efficient Stokesian and Brownian dynamics simulations is outlined. An example of an application is shown.     

Size segregation in a fluid-like or gel-like suspension settling under gravity or in a centrifuge

We investigate size segregation effects in a bidisperse concentrated suspension when slowly settling under gravity or when submitted to a centrifugal field. Experiments are carried out with PMMA spheres of two different mean diameters (190 and 25 microm) suspended in a hydrophobic index-matched fluid. Spatial repartitions of both small and large spheres and velocity fluctuations of particles are measured using fluorescently labeled PMMA spheres and a particle-image-velocimetry method. Large particles behave as hard spheres in purely hydrodynamic interactions, while small spheres interact through weakly attractive forces. For a small amount of small spheres among large ones, the suspension remains fluid during settling and the organization of the velocity field of particles into finite-sized structures also called "blobs" promotes size segregation. A larger proportion of weakly attractive small spheres in the bidisperse suspension causes a considerable slowdown of the settling process under gravity and the occurrence of a large-scale collective behavior together with a loss of size segregation. When centrifuging the gel-like bidisperse suspension, a shear-induced melting of the particle network induces a spectacular segregation of species. As a consequence, aging tests of soft yielding materials using centrifugation methods are not representative of the shelf-life stability of the products. A tentative model based on the competition between viscous stresses acting upon particles and adhesive stresses gives a correct estimate of the critical stationary acceleration for the destabilization of the particle network and the onset of size segregation in a gel-like suspension.     

Low Reynolds Number Interactions between Colloidal Particles near the Entrance to a Cylindrical Pore

The interaction between stable colloidal particles arriving at a pore entrance was studied using a numerical method for the case where the particle size is smaller than but of the same order as the pore size. The numerical method was adapted from a front-tracking technique developed for studying incompressible, multifluid flow by S. O. Unverdi and G. Tryggvason (J. Comp. Phys. 100, 25, 1992). The method is based on the finite difference solution of Navier-Stokes equation on a stationary, structured, Cartesian grid and the explicit representation of the particle-liquid interface using an unstructured grid that moves through the stationary grid. The simulations are in two dimensions, considering both deformable and nondeformable particles, and include interparticle colloidal interactions. The interparticle and particle-pore hydrodynamic interactions, which are very difficult to determine using existing analytical and semi-numerical, semi-analytical techniques in microhydrodynamics, are naturally accounted for in our numerical method and need not be explicity determined. Two- and three-particle motion toward a pore has been considered in our simulations. The simulations demonstrate how the competition between hydrodynamic forces and colloidal forces acting on particles dictate their flow behavior near the pore entrance. The predicted dependence of the particle flow behavior on the flow velocity and the ratio of pore size to particle size are qualitatively consistent with the experimental observations of V. Ramachandran and H. S. Fogler (J. Fluid Mech. 385, 129, 1999). Copyright 2000 Academic Press.     

Effect of adsorbed polymers on the slow motion of an assemblage of spherical particles relative to a fluid

 The effects of adsorbed polymers on the sedimentation of a homogeneous distribution of colloidal spheres and on the fluid flow through a bed of particles are investigated theoretically. The Reynolds number is assumed to be small, and the surface polymer layers are assumed to be thin with respect to the radius of particles and to the surface-to-surface spacing between neighboring particles. The effects of interaction of the individual particles are taken into explicit account by employing a fundamental cell-model representation which is known to provide good predictions for the motion of a swarm of spheres within a fluid in the absence of adsorbed polymers. To solve the Stokes flow equations within and outside the polymer layer a method of matched asymptotic expansions in a small parameter λ is used, where λ is the ratio of the length scale of the polymer layer to the particle radius. The results for the sedimentation rate and the pressure drop are expressed in terms of an effective hydrodynamic thickness (L) of the polymer layer, which are accurate to O(λ²). When the concentration of particles in a suspension or a bed is increased, L becomes larger, meaning the settling velocity decreases or the pressure drop increases. The O(λ) term for L normalized by its value in the limit λ→0 is found to be independent of the polymer segment distribution, the hydrodynamic inter-actions among the segments, and the volume fraction of the segments. The O(λ²) term for L, however, is a sensitive function of the polymer segment distribution and the volume fraction of the segments. In general, the particle-interaction effects on the motion of polymer-coated particles relative to a fluid can be quite significant. Received: 28 August 1996 Accepted: 23 January 1997     

Collective diffusion in charge-stabilized suspensions: concentration and salt effects

The authors present a joint experimental-theoretical study of collective diffusion properties in aqueous suspensions of charge-stabilized fluorinated latex spheres. Small-angle x-ray scattering and x-ray photon correlation spectroscopy have been used to explore the concentration and ionic-strength dependence of the static and short-time dynamic properties including the hydrodynamic function H(q), the wave-number-dependent collective diffusion coefficient D(q), and the intermediate scattering function over the entire accessible range. They show that all experimental data can be quantitatively described and explained by means of a recently developed accelerated Stokesian dynamics simulation method, in combination with a modified hydrodynamic many-body theory. In particular, the behavior of H(q) for de-ionized and dense suspensions can be attributed to the influence of many-body hydrodynamics, without any need for postulating hydrodynamic screening to be present, as it was done in earlier work. Upper and lower boundaries are provided for the peak height of the hydrodynamic function and for the short-time self-diffusion coefficient over the entire range of added salt concentrations.     

Strong magnetic field effects on solid-liquid and particle-particle interactions during the processing of a conducting liquid containing non-conducting particles

The behavior of micrometer-sized weak magnetic insulating particles migrating in a conductive liquid metal is of broad interest during strong magnetic field processing of materials. In the present paper, we develop a numerical method to investigate the solid-liquid and particle-particle interactions by using a computational fluid dynamics (CFDs) modeling. By applying a strong magnetic field, for example, 10 Tesla, the drag forces of a single spherical particle can be increased up to around 15% at a creeping flow limit. However, magnetic field effects are reduced when the Reynolds number becomes higher. For two identical particles migrating along their centerline in a conductive liquid, both the drag forces and the magnetic interaction will be influenced. Factors such as interparticle distance, Reynolds number and magnetic flux density are investigated. Shielding effects are found from the leading particle, which will subsequently induce a hydrodynamic interaction between two particles. Strong magnetic fields however do not appear to have a significant influence on the shielding effects. In addition, the magnetic interaction forces of magnetic dipole-dipole interaction and induced magneto-hydrodynamic interaction are considered. It can be found that the induced magneto-hydrodynamic interaction force highly depends on the flow field and magnetic flux density. Therefore, the interaction between insulating particles can be controlled by applying a strong magnetic field and modifying the flow field. The present research provides a better understanding of the magnetic field induced interaction during liquid metal processing, and a method of non-metallic particles manipulation for metal/ceramic based materials preparation may be proposed.     

Short-time rheology and diffusion in suspensions of Yukawa-type colloidal particles

A comprehensive study is presented on the short-time dynamics in suspensions of charged colloidal spheres. The explored parameter space covers the major part of the fluid-state regime, with colloid concentrations extending up to the freezing transition. The particles are assumed to interact directly by a hard-core plus screened Coulomb potential, and indirectly by solvent-mediated hydrodynamic interactions. By comparison with accurate accelerated Stokesian Dynamics (ASD) simulations of the hydrodynamic function H(q), and the high-frequency viscosity η(∞), we investigate the accuracy of two fast and easy-to-implement analytical schemes. The first scheme, referred to as the pairwise additive (PA) scheme, uses exact two-body hydrodynamic mobility tensors. It is in good agreement with the ASD simulations of H(q) and η(∞), for smaller volume fractions up to about 10% and 20%, respectively. The second scheme is a hybrid method combining the virtues of the δγ scheme by Beenakker and Mazur with those of the PA scheme. It leads to predictions in good agreement with the simulation data, for all considered concentrations, combining thus precision with computational efficiency. The hybrid method is used to test the accuracy of a generalized Stokes-Einstein (GSE) relation proposed by Kholodenko and Douglas, showing its severe violation in low salinity systems. For hard spheres, however, this GSE relation applies decently well.     

Precise multipole method for calculating many-body hydrodynamic interactions in a microchannel

We introduce a novel and precise method for computing many-body hydrodynamic interactions in a cylindrical microchannel. The method is generic in the sense that we can easily change the radius and the character of particles (hard spheres, droplets, permeable spheres, etc.). These features are not available in any of the existing methods. Comparison with the available results validates our method. In particular we obtain excellent agreement with the analytically known expression for the single particle friction coefficient. Additionally we observe negative hydrodynamic coupling for finite particles which are consistent with the recently reported effect for point particles. As an example we compute the velocities of polymeric chains of particles in parabolic flow and compare them to unbounded space. The method will be helpful in the understanding of physical and physicochemical processes in a wide range of bio-, geophysical, and microfluidic systems.     

Hydrodynamic crossover in dynamic microparticle adhesion on surfaces of controlled nanoscale heterogeneity

This note documents the crossover from a regime where shear flow hinders microparticle adhesion on collecting surfaces to that where increased flow aids particle capture. Flow generally works against adhesion and successfully hinders particle capture when the net physicochemical attractions between the particles and collector are weak compared with hydrodynamic forces on the particle. Conversely, with strong attractions between particles and collector, flow aids particle capture by increasing the mass transport of particles to the interfacial region. Here, local hydrodynamics still generally oppose adhesion but are insufficient to pull particles off of the surface. Thus, flow actually increases the particle capture rate through the increased transport to the surface. These behaviors are demonstrated using 1 mum silica spheres flowing over electrostatically heterogeneous (length scales near 10 nm) collecting surfaces at shear rates from 22 to 795 s(-1). The net surface charge on the collector is varied systematically from strongly negative (pure silica) to strongly positive (a saturated polycationic overlayer), demonstrating the interplay between physicochemical and hydrodynamic contributions. These results clearly apply to situations where heterogeneous particle-surface interactions are electrostatic in nature; however, qualitatively similar behavior was previously reported for the effect receptor density on bacterial adhesion.     

Modeling of solvent flow effects in enzyme catalysis under physiological conditions

A stochastic model for the dynamics of enzymatic catalysis in explicit, effective solvents under physiological conditions is presented. Analytically-computed first passage time densities of a diffusing particle in a spherical shell with absorbing boundaries are combined with densities obtained from explicit simulation to obtain the overall probability density for the total reaction cycle time of the enzymatic system. The method is used to investigate the catalytic transfer of a phosphoryl group in a phosphoglycerate kinase-ADP-bis phosphoglycerate system, one of the steps of glycolysis. The direct simulation of the enzyme-substrate binding and reaction is carried out using an elastic network model for the protein, and the solvent motions are described by multiparticle collision dynamics which incorporates hydrodynamic flow effects. Systems where solvent-enzyme coupling occurs through explicit intermolecular interactions, as well as systems where this coupling is taken into account by including the protein and substrate in the multiparticle collision step, are investigated and compared with simulations where hydrodynamic coupling is absent. It is demonstrated that the flow of solvent particles around the enzyme facilitates the large-scale hinge motion of the enzyme with bound substrates, and has a significant impact on the shape of the probability densities and average time scales of substrate binding for substrates near the enzyme, the closure of the enzyme after binding, and the overall time of completion of the cycle.