Treating inertia in passive microbead rheology

Phys Rev E Stat Nonlin Soft Matter Phys. 2012 Feb;85(2 Pt 1):021504. doi: 10.1103/PhysRevE.85.021504. Epub 2012 Feb 24.

Abstract

The dynamic modulus G(*) of a viscoelastic medium is often measured by following the trajectory of a small bead subject to Brownian motion in a method called "passive microbead rheology." This equivalence between the positional autocorrelation function of the tracer bead and G(*) is assumed via the generalized Stokes-Einstein relation (GSER). However, inertia of both bead and medium are neglected in the GSER so that the analysis based on the GSER is not valid at high frequency where inertia is important. In this paper we show how to treat both contributions to inertia properly in one-bead passive microrheological analysis. A Maxwell fluid is studied as the simplest example of a viscoelastic fluid to resolve some apparent paradoxes of eliminating inertia. In the original GSER, the mean-square displacement (MSD) of the tracer bead does not satisfy the correct initial condition. If bead inertia is considered, the proper initial condition is realized, thereby indicating an importance of including inertia, but the MSD oscillates at a time regime smaller than the relaxation time of the fluid. This behavior is rather different from the original result of the GSER and what is observed. What is more, the discrepancy from the GSER result becomes worse with decreasing bead mass, and there is an anomalous gap between the MSD derived by naïvely taking the zero-mass limit in the equation of motion and the MSD for finite bead mass as indicated by McKinley et al. [J. Rheol. 53, 1487 (2009)]. In this paper we show what is necessary to take the zero-mass limit of the bead safely and correctly without causing either the inertial oscillation or the anomalous gap, while obtaining the proper initial condition. The presence of a very small purely viscous element can be used to eliminate bead inertia safely once included in the GSER. We also show that if the medium contains relaxation times outside the window where the single-mode Maxwell behavior is observed, the oscillation can be attenuated inside the window. This attenuation is realized even in the absence of a purely viscous element. Finally, fluid inertia also affects the bead autocorrelation through the Basset force and the fluid dragged around with the bead. We show that the Basset force plays the same role as the purely viscous element in high-frequency regime, and the oscillation of MSD is suppressed if fluid density and bead density are comparable.

Publication types

  • Research Support, U.S. Gov't, Non-P.H.S.

MeSH terms

  • Computer Simulation
  • Diffusion*
  • Elastic Modulus
  • Microspheres*
  • Models, Chemical*
  • Models, Molecular*
  • Rheology
  • Viscosity