pith. sign in

arxiv: 2508.13022 · v1 · submitted 2025-08-18 · ⚛️ physics.flu-dyn · cond-mat.soft

On the Rheology of Two-Dimensional Dilute Emulsions

Pith reviewed 2026-05-18 22:17 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn cond-mat.soft
keywords two-dimensional emulsionsdilute suspensionsrheologydroplet deformationStokes flowapparent viscositycapillary numberviscosity ratio
0
0 comments X

The pith

Two-dimensional dilute emulsions have apparent viscosity μ(1 + ((2λ + 1)/(λ + 1))φ) and steady droplet deformation equal to the capillary number.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives analytical expressions for the rheology and small-deformation behavior of droplets suspended in a two-dimensional fluid. It begins from the Lamb solution for two-dimensional Stokes flow to obtain the velocity and pressure fields around a droplet in extensional flow, then integrates these fields to find the effective viscosity and the steady Taylor deformation. The resulting viscosity increases linearly with droplet area fraction according to a factor that depends on the viscosity ratio λ, while the deformation parameter reaches a steady value exactly equal to the capillary number with no λ dependence. A reader would care because two-dimensional simulations are widely used for their computational speed, and these closed-form results supply exact benchmarks that three-dimensional analogs lack.

Core claim

Using the Lamb solution for two-dimensional Stokes flows, the flow fields around a droplet in a purely extensional flow are obtained. These yield the apparent viscosity μ* = μ(1 + f(λ)φ) + O(φ²) where f(λ) = (2λ + 1)/(λ + 1), and the steady-state Taylor deformation D_T^∞ = Ca in the capillarity-dominated regime, with the prefactor g(λ) equal to one. This contrasts with the three-dimensional case in which g(λ) depends on λ. The results are validated by direct numerical simulations for viscosity ratios between 0.01 and 100.

What carries the argument

The Lamb solution for two-dimensional Stokes flows applied to a viscous droplet in pure extensional flow, which supplies the velocity and pressure fields needed to compute effective viscosity and deformation.

If this is right

  • The apparent viscosity of the emulsion increases linearly with area fraction φ according to the exact prefactor (2λ + 1)/(λ + 1).
  • The steady Taylor deformation parameter equals the capillary number with no dependence on the viscosity ratio λ.
  • The derived expressions supply precise benchmarks that can be used to verify numerical codes for two-dimensional droplet problems.
  • The analytical results remain accurate across viscosity ratios from 0.01 to 100 when tested against direct simulations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same Lamb-solution approach could be applied to other two-dimensional flows such as simple shear to obtain analogous rheology expressions.
  • The lack of λ dependence in the deformation result may allow simpler reduced-order models for droplet dynamics in two-dimensional microfluidic geometries.
  • Extensions that relax the dilute assumption could quantify pairwise droplet interactions and their effect on the O(φ²) viscosity correction.
  • Two-dimensional models might therefore serve as efficient test beds for exploring emulsion behavior before moving to more expensive three-dimensional calculations.

Load-bearing premise

The assumptions of small droplet deformation and dilute suspension are invoked to obtain closed analytical expressions for viscosity and deformation after the flow fields are known.

What would settle it

A direct numerical simulation at low capillary number that shows the steady Taylor deformation parameter differing from the capillary number for any viscosity ratio in the range 0.01 to 100.

Figures

Figures reproduced from arXiv: 2508.13022 by Maziyar Jalaal, Thomas Appleford, Vatsal Sanjay.

Figure 1
Figure 1. Figure 1: Configuration of the simulations of the droplet under shear problem. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) The apparent viscosity, µ ∗, of the system as a function of the viscosity ratio, λ. The black line shows the result derived from the analytical solution. The red circular points show the results obtained from numerical simulations. Here, L = 20R such that the volume fraction is ϕ ≈ 0.008. Also shown are the asymptotic limits λ → 0 and λ → ∞. (b) The apparent viscosity, µ ∗, of the system as a function … view at source ↗
Figure 3
Figure 3. Figure 3: The norm, ∥E∥ of the strain rate tensor for non-deformable droplets (Ca = 0.01) with different values of λ. Panels (a), (b) and (c) show λ = 0.01, 1 and 100 respectively. The streamlines depict the vector field, u − R∞ · x, that is the velocity field with the solid-body rotational part of the undisturbed flow subtracted. (a) (b) (c) 0 0.5 1 1.5 2 kEk/ ˙γ [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The norm, ∥E∥, of the strain rate tensor for deformable droplets (Ca = 0.3) with different values of λ. Panels (a), (b) and (c) show λ = 0.01, 1 and 100, respectively. The streamlines depict the vector field, u − R∞ · x, that is the velocity field with the solid-body rotational part of the undisturbed flow subtracted. V. CONCLUSIONS We have derived an analytical solution for the two-dimensional droplet und… view at source ↗
Figure 5
Figure 5. Figure 5: The dependence of the Taylor deformation parameter, [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Validation of numerical simulations (circles) against the analytical solution Eqs. ( [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: (a) Convergence of the apparent viscosity [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
read the original abstract

The single droplet under shear is a foundational problem in fluid mechanics. In computational fluid dynamics, the two-dimensional (2D) formulation offers advantages in both computational efficiency and relevance, yet its theoretical treatment remains relatively underdeveloped. In this brief note, we present an analytical treatment of this problem, beginning with a derivation of the Lamb solution for 2D Stokes flows, which in turn is used to obtain the flow fields around a droplet in a purely extensional flow. Using these flow fields, expressions are obtained for the apparent viscosity, $\mu^*$, of a dilute emulsion as well as a small deformation theory. We show that $\mu^* = \mu( 1 + f(\lambda) \phi) + \mathcal{O}(\phi^2)$ with $f(\lambda) = (2\lambda + 1)/(\lambda + 1)$ where $\lambda$ is the ratio of the droplet viscosity to the matrix viscosity and $\phi$ is the area fraction covered by the suspended phase. Also the steady state value of the Taylor deformation parameter $D_T^\infty$, in the capillarity-dominated regime, obeys $D_T^\infty = g(\lambda)\,\text{Ca}$, where Ca is the capillary number and $g(\lambda) = 1$. This contrasts with the 3D case, where $g(\lambda)$ depends on $\lambda$. These results are then validated through direct numerical simulations across a wide range of viscosity ratios ($0.01 < \lambda < 100$). Our results provide a basic theoretical framework for interpreting 2D droplet simulations and provide clear benchmarks for computational fluid dynamics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript derives analytical expressions for the rheology of two-dimensional dilute emulsions by solving the 2D Stokes equations for a circular droplet in far-field extensional flow using the Lamb general solution. It obtains the apparent viscosity μ* = μ (1 + ((2λ + 1)/(λ + 1)) φ) + O(φ²) with f(λ) = (2λ + 1)/(λ + 1) and the steady Taylor deformation D_T^∞ = Ca (independent of λ) in the capillarity-dominated regime, validated by direct numerical simulations over 0.01 < λ < 100.

Significance. If the central results hold, this brief note supplies a parameter-free theoretical framework and clear benchmarks for 2D droplet simulations, which offer computational advantages. The direct derivation from the 2D Stokes solution without fitted constants, together with the DNS confirmation across a wide viscosity-ratio range, constitutes a useful contribution that highlights the λ-independence of deformation in 2D (in contrast to 3D).

minor comments (2)
  1. [deformation analysis] The manuscript would benefit from a short statement of the expected range of validity for the small-deformation closure used to obtain D_T^∞ = Ca (e.g., an order-of-magnitude estimate of Ca below which the O(Ca) result remains accurate).
  2. [validation section] DNS comparisons would be strengthened by reporting quantitative measures of agreement (e.g., relative error or L2 deviation) between the simulated viscosity and deformation values and the analytical predictions, rather than qualitative statements of match.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript, accurate summary of our results, and positive recommendation to accept. We appreciate their recognition of the utility of the parameter-free 2D analytical framework and the DNS benchmarks across a wide range of viscosity ratios.

Circularity Check

0 steps flagged

Derivation is self-contained from first-principles 2D Stokes solution

full rationale

The central results for apparent viscosity μ* and steady Taylor deformation D_T^∞ are obtained by solving the 2D Stokes equations inside and outside the droplet using the general Lamb expansion for 2D flow, then enforcing velocity continuity, tangential stress continuity, and normal stress jump due to curvature. The small-deformation limit and far-field stresslet directly yield f(λ) = (2λ + 1)/(λ + 1) and g(λ) = 1 with no fitted parameters or reduction to prior results. DNS validation is supplied as independent confirmation. No self-citations, ansatzes, or fitted inputs are load-bearing for the O(φ) and O(Ca) expressions.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard low-Reynolds-number incompressible flow assumptions and the small-deformation approximation; no new free parameters or postulated entities are introduced.

axioms (2)
  • standard math Stokes equations govern the flow at vanishing Reynolds number
    Basis for the Lamb solution derivation in the opening sections.
  • domain assumption Droplet deformation remains small so that shape can be treated perturbatively
    Required to obtain the steady Taylor deformation expression.

pith-pipeline@v0.9.0 · 5831 in / 1139 out tokens · 39172 ms · 2026-05-18T22:17:44.120619+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 53 canonical work pages

  1. [1]

    Rheology of emulsions

    S. R. Derkach, “Rheology of emulsions”, Adv. Colloid Interface Sci. 151, 1–23 (2009)

  2. [2]

    The deformation of small viscous drops and bubbles in shear flows

    J. M. Rallison, “The deformation of small viscous drops and bubbles in shear flows”, Annu. Rev. Fluid Mech. 16, 45–66 (1984)

  3. [3]

    Stretching and breakup of droplets in chaotic flows

    M. Tjahjadi and J. M. Ottino, “Stretching and breakup of droplets in chaotic flows”, J. Fluid Mech. 232, 191–219 (1991)

  4. [4]

    The viscosity of a fluid containing small drops of another fluid

    G. I. Taylor, “The viscosity of a fluid containing small drops of another fluid”, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 138, 41–48 (1932)

  5. [5]

    The formation of emulsions in definable fields of flow

    G. I. Taylor, “The formation of emulsions in definable fields of flow”, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 146, 501–523 (1934)

  6. [6]

    Lamb, Cambridge mathematical library: hydrodynamics , en, 6th ed

    H. Lamb, Cambridge mathematical library: hydrodynamics , en, 6th ed. (Cambridge University Press, Cambridge, England, Nov. 1993)

  7. [7]

    Eine neue bestimmung der molek¨ uldimensionen

    A. Einstein, “Eine neue bestimmung der molek¨ uldimensionen”, PhD thesis (ETH Zurich, 1905)

  8. [8]

    The deformation of a drop in a general time-dependent fluid flow

    R. G. Cox, “The deformation of a drop in a general time-dependent fluid flow”, Journal of Fluid Mechanics 37, 601–623 (1969)

  9. [9]

    A second-order theory for shear deformation of drops

    C. E. Chaffey and H. Brenner, “A second-order theory for shear deformation of drops”, Journal of Colloid and Interface Science 24, 258–269 (1967)

  10. [10]

    Equation of change for ellipsoidal drops in viscous flow

    P. Maffettone and M. Minale, “Equation of change for ellipsoidal drops in viscous flow”, Journal of Non-Newtonian Fluid Mechanics 78, 227–241 (1998)

  11. [11]

    Note on the time-dependent deformation of a viscous drop which is almost spherical

    J. M. Rallison, “Note on the time-dependent deformation of a viscous drop which is almost spherical”, Journal of Fluid Mechanics 98, 625–633 (1980)

  12. [12]

    Droplet dynamics in confinement

    N. Ioannou, H. Liu, and Y. Zhang, “Droplet dynamics in confinement”, Journal of Computational Science 17, 463–474 (2016)

  13. [13]

    The effects of surfactants on drop deformation and breakup

    H. A. Stone and L. G. Leal, “The effects of surfactants on drop deformation and breakup”, Journal of Fluid Mechanics 220, 161–186 (1990)

  14. [14]

    Small-deformation theory for a surfactant-covered drop in linear flows

    P. M. Vlahovska, J. B/suppress lawzdziewicz, and M. Loewenberg, “Small-deformation theory for a surfactant-covered drop in linear flows”, Journal of Fluid Mechanics 624, 293–337 (2009)

  15. [15]

    Deformation of clean and surfactant-laden droplets in shear flow

    G. Soligo, A. Roccon, and A. Soldati, “Deformation of clean and surfactant-laden droplets in shear flow”, Meccanica 55, 371–386 (2020)

  16. [16]

    Numerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method

    J. Li, Y. Y. Renardy, and M. Renardy, “Numerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method”, Physics of Fluids 12, 269–282 (2000)

  17. [17]

    An experimental investigation of drop deformation and breakup in steady, two-dimensional linear flows

    B. Bentley and L. G. Leal, “An experimental investigation of drop deformation and breakup in steady, two-dimensional linear flows”, Journal of Fluid Mechanics 167, 241–283 (1986)

  18. [18]

    Shear-induced droplet deformation: effects of confined geometry and viscoelasticity

    S. Guido, “Shear-induced droplet deformation: effects of confined geometry and viscoelasticity”, en, Curr. Opin. Colloid Interface Sci. 16, 61–70 (2011)

  19. [19]

    Shear-induced rupturing of a viscous drop in a bingham liquid

    J. Li and Y. Y. Renardy, “Shear-induced rupturing of a viscous drop in a bingham liquid”, Journal of Non-Newtonian Fluid Mechanics 95, 235–251 (2000)

  20. [20]

    Droplet deformation and breakup in shear-thinning viscoelastic fluid under simple shear flow

    D. Wang, N. Wang, and H. Liu, “Droplet deformation and breakup in shear-thinning viscoelastic fluid under simple shear flow”, Journal of Rheology 66, 585–603 (2022)

  21. [21]

    The breakup of small drops and bubbles in shear flows

    A. Acrivos, “The breakup of small drops and bubbles in shear flows”, Annals of the New York Academy of Sciences 404, 1–11 (1983)

  22. [22]

    Dynamics of drop deformation and breakup in viscous fluids

    H. A. Stone, “Dynamics of drop deformation and breakup in viscous fluids”, Annu. Rev. Fluid Mech. 26, 65–102 (1994)

  23. [23]

    Models for the deformation of a single ellipsoidal drop: a review

    M. Minale, “Models for the deformation of a single ellipsoidal drop: a review”, Rheologica acta 49, 789–806 (2010)

  24. [24]

    Modeling drop deformations and rheology of dilute to dense emulsions

    R. B. Reboucas, N. N. Nikolova, and V. Sharma, “Modeling drop deformations and rheology of dilute to dense emulsions”, Current Opinion in Colloid & Interface Science, 101904 (2025)

  25. [25]

    Flow of deformable droplets: discontinuous shear thinning and velocity oscillations

    M. Foglino, A. Morozov, O Henrich, and D Marenduzzo, “Flow of deformable droplets: discontinuous shear thinning and velocity oscillations”, Physical review letters 119, 208002 (2017)

  26. [26]

    Yield-stress transition in suspensions of deformable droplets

    G. Negro, L. N. Carenza, G. Gonnella, F. Mackay, A. Morozov, and D. Marenduzzo, “Yield-stress transition in suspensions of deformable droplets”, Science Advances 9, eadf8106 (2023)

  27. [27]

    Two-dimensional bubbles in slow viscous flows

    S Richardson, “Two-dimensional bubbles in slow viscous flows”, Journal of Fluid Mechanics 33, 475–493 (1968)

  28. [28]

    Two-dimensional bubbles in slow viscous flows. part 2

    S Richardson, “Two-dimensional bubbles in slow viscous flows. part 2”, Journal of Fluid Mechanics 58, 115–127 (1973)

  29. [29]

    The bursting of two-dimensional drops in slow viscous flow

    J. D. Buckmaster and J. E. Flaherty, “The bursting of two-dimensional drops in slow viscous flow”, J. Fluid Mech. 60, 625–639 (1973)

  30. [30]

    A boundary integral method for two-dimensional (non)-newtonian drops in slow viscous flow

    E. Toose, B. Geurts, and J. Kuerten, “A boundary integral method for two-dimensional (non)-newtonian drops in slow viscous flow”, Journal of non-newtonian fluid mechanics 60, 129–154 (1995)

  31. [31]

    Emulsion droplet deformation and breakup with lattice boltzmann model

    R. G. van der Sman and S Van der Graaf, “Emulsion droplet deformation and breakup with lattice boltzmann model”, Computer Physics Communications 178, 492–504 (2008)

  32. [32]

    Hydrodynamic lift of a two-dimensional liquid domain with odd viscosity

    Y. Hosaka, S. Komura, and D. Andelman, “Hydrodynamic lift of a two-dimensional liquid domain with odd viscosity”, Physical Review E 104, 064613 (2021)

  33. [33]

    Odd droplets: fluids with odd viscosity and highly deformable interfaces

    H. Fran¸ ca and M. Jalaal, “Odd droplets: fluids with odd viscosity and highly deformable interfaces”, arXiv preprint arXiv:2503.21649 (2025)

  34. [34]

    The translational and rotational motions of an n-dimensional hypersphere through a viscous fluid at small reynolds numbers

    H. Brenner, “The translational and rotational motions of an n-dimensional hypersphere through a viscous fluid at small reynolds numbers”, Journal of Fluid Mechanics 111, 197–215 (1981)

  35. [35]

    G. K. Batchelor, An introduction to fluid dynamics (Cambridge university press, 2000). 14

  36. [36]

    The stress system in a suspension of force-free particles

    G. Batchelor, “The stress system in a suspension of force-free particles”, Journal of fluid mechanics 41, 545–570 (1970)

  37. [37]

    Capillary-scale solid rebounds: experiments, modelling and simulations

    C. A. Galeano-Rios, R. Cimpeanu, I. A. Bauman, A. MacEwen, P. A. Milewski, and D. M. Harris, “Capillary-scale solid rebounds: experiments, modelling and simulations”, J. Fluid Mech. 912, A17 (2021)

  38. [38]

    The einstein viscosity correction in n dimensions

    J. F. Brady, “The einstein viscosity correction in n dimensions”, International Journal of Multiphase Flow 10, 113–114 (1983)

  39. [39]

    Gerris: a tree-based adaptive solver for the incompressible euler equations in complex geometries

    S. Popinet, “Gerris: a tree-based adaptive solver for the incompressible euler equations in complex geometries”, Journal of Computational Physics 190, 572–600 (2003)

  40. [40]

    An accurate adaptive solver for surface-tension-driven interfacial flows

    S. Popinet, “An accurate adaptive solver for surface-tension-driven interfacial flows”, Journal of Computational Physics 228, 5838–5866 (2009)

  41. [41]

    A quadtree-adaptive multigrid solver for the serre–green–naghdi equations

    S. Popinet, “A quadtree-adaptive multigrid solver for the serre–green–naghdi equations”, Journal of Computational Physics 302, 336–358 (2015)

  42. [42]

    Viscous free-surface flows

    V. Sanjay, “Viscous free-surface flows”, PhD thesis (Univ. Twente, 2022)

  43. [43]

    Fluid mechanics and rheology of dense suspensions

    J. J. Stickel and R. L. Powell, “Fluid mechanics and rheology of dense suspensions”, Annual Review of Fluid Mechanics 37, 129–149 (2005)

  44. [44]

    Deformation and burst of a liquid droplet freely suspended in a linear shear field

    D. Barth` es-Biesel and A. Acrivos, “Deformation and burst of a liquid droplet freely suspended in a linear shear field”, Journal of Fluid Mechanics 61, 1–22 (1973)

  45. [45]

    Deformation and breakup of a viscoelastic drop in a newtonian matrix under steady shear

    N. AGGAR W AL and K. SARKAR, “Deformation and breakup of a viscoelastic drop in a newtonian matrix under steady shear”, Journal of Fluid Mechanics 584, 1–21 (2007)

  46. [46]

    Effects of matrix viscoelasticity on viscous and viscoelastic drop deformation in a shear flow

    N. AGGAR W AL and K. SARKAR, “Effects of matrix viscoelasticity on viscous and viscoelastic drop deformation in a shear flow”, Journal of Fluid Mechanics 601, 63–84 (2008)

  47. [47]

    Rheology of an emulsion of viscoelastic drops in steady shear

    N. Aggarwal and K. Sarkar, “Rheology of an emulsion of viscoelastic drops in steady shear”, Journal of Non-Newtonian Fluid Mechanics 150, 19–31 (2008)

  48. [48]

    Dynamics of an elastoviscoplastic droplet in a newtonian medium under shear flow

    D. Izbassarov and O. Tammisola, “Dynamics of an elastoviscoplastic droplet in a newtonian medium under shear flow”, Phys. Rev. Fluids 5, 113301 (2020)

  49. [49]

    Elasto-viscoplastic spreading: from plastocapillarity to elastocapillarity

    H. L. Fran¸ ca, M. Jalaal, and C. M. Oishi, “Elasto-viscoplastic spreading: from plastocapillarity to elastocapillarity”, Physical Review Research 6, 013226 (2024)

  50. [50]

    The effect of surfactants on drop deformation and on the rheology of dilute emulsions in stokes flow

    X. Li and C Pozrikidis, “The effect of surfactants on drop deformation and on the rheology of dilute emulsions in stokes flow”, Journal of fluid mechanics 341, 165–194 (1997)

  51. [51]

    Shape and rheology of droplets with viscous surface moduli

    V. Narsimhan, “Shape and rheology of droplets with viscous surface moduli”, Journal of Fluid Mechanics 862, 385–420 (2019)

  52. [52]

    Odd viscosity

    J. E. A vron, “Odd viscosity”, Journal of Statistical Physics 92, 543–557 (1998)

  53. [53]

    Odd viscosity and odd elasticity

    M. Fruchart, C. Scheibner, and V. Vitelli, “Odd viscosity and odd elasticity”, Annual Review of Condensed Matter Physics 14, 471–510 (2023). 15