Correction to Doi type models for suspensions
David Gérard-Varet, Richard M. Höfer
TL;DR
The paper derives a corrected Doi-type model for suspensions of non-spherical particles in Stokes flow by starting from a microscopic $N$-particle system and addressing the breakdown of the naive mean-field limit for orientation dynamics due to a singular $-3$-homogeneous interaction. A deterministic intermediate system is analyzed via the method of reflections and an energy-based regularization of the singular term, yielding a refined angular dynamics that includes a new term $\phi M B\xi f$ with a microstructure-dependent tensor $B$. When particle configurations are generated by a stationary ergodic marked point process, the singular interaction limit becomes almost surely deterministic and expressible through 1- and 2-point correlations; under mixing, this limit has an explicit representation in terms of the correlations and the flow maps $\Phi$, $\Xi$, enabling a concrete form for $B$. The main result, encapsulated in Theorem thm.f-f_phi, shows that the microscopic empirical measure $f_N$ is approximated to first order in $\phi$ by $f_\phi$, with an $o(\phi)$ error in a negative Sobolev topology, thereby providing a practically usable correction to Doi-type models in dilute suspensions. This links the orientation dynamics directly to particle correlations and two-scale homogenization, suggesting that microstructure information is essential for accurately predicting orientational phenomena in suspensions of nonspherical particles.
Abstract
Starting from microscopic $N$ particle systems, we study the derivation of Doi type models for suspensions of non-spherical particles in Stokes flows. While Doi models accurately describe the effective evolution of the spatial particle density to the first order in the particle volume fraction, this accuracy fails regarding the evolution of the particle orientations. We rigorously attribute this failure to the singular interaction of the particles via a $-3$-homogeneous kernel. In the situation that the particles are initially distributed according to a stationary ergodic point process, we identify the limit of this singular interaction term. It consists of two parts. The first corresponds to a classical term in Doi type models. The second new term depends on the (microscopic) $2$-point correlation of the point process. By including this term, we provide a modification of the Doi model that is accurate to first order in the particle volume fraction.
