REVIEW 3 major objections 5 minor 35 references
On the effect of coalescence on the rheology of emulsions
T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Coalescence, not just droplet deformation, is what makes emulsion viscosity rise to a peak near 20% volume fraction and then fall.
desk verdict A solid VoF study showing coalescence flips the curvature of the emulsion viscosity curve, but the artificial repulsive force that defines the 'no-coalescence' branch is not a fully clean control. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing device is an Eulerian repulsive force of lubrication form, Fc = μ0 r U (a/ψ + b/ψ²) n, written in terms of the signed distance ψ from the interface; it is applied only where two different droplet indices are detected within a stencil, and it models the effect of surfactants by preventing film drainage and merging. The force is hand-tuned (a = 55, b = 3.5) so that coalescence is fully suppressed with minimal side effects. The second piece is the stress-budget identity that equates the wall shear stress with the bulk average of viscous, interfacial-tension, and collision-force stresses; combined with the time history of the total interface area, this identity is what lets the authors attribute the viscosity change to a reduction of interfacial area. A third element is the flow-topology parameter Q = (D² − Ω²)/(D² + Ω²), used to show that the matrix flow is nearly pure shear when droplets merge and gains extensional regions when merging is suppressed.
What would settle it
A decisive check is to repeat the volume-fraction series with a second, force-free way of suppressing coalescence—for example an interface model with insoluble surfactant transport and strong Marangoni resistance—and compare the effective-viscosity curve and stress budget. If the positive curvature and the roughly 10% collision-force stress contribution do not reproduce, the reported curvature flip is an artifact of the repulsive model. A more local test: in the two-droplet collision of figure 4, the forced case differs by about 6% in final vertical displacement from the double-resolution unforced case; refining the grid should reduce this difference if the force is inert, or expose that the force itself alters the collision outcome.
Extended reading notes
Core claim
The central claim is that droplet coalescence causes the negative curvature of the effective-viscosity versus volume-fraction curve in emulsions. In the simulations, free merging produces a maximum of the effective viscosity near φ ≈ 0.2; for viscosity ratio λ = 0.01 the normalized viscosity even falls below 1 at high volume fraction. Applying the repulsive collision force (3.1) with coefficients a = 55 and b = 3.5 prevents merging altogether, and the same system then shows a monotonic, positively curved viscosity curve that is well described by the Eilers formula once droplet deformation is folded into an effective volume fraction. The paper traces this to the stress budget: with coalescence the total interface area drops by up to 80% of its initial value, lowering the interfacial-tension contribution to the shear stress; without coalescence the interface area stays nearly constant, and the interface-tension term accounts for about half of the effective viscosity, with the collision force itself contributing about 10%.
Load-bearing premise
The load-bearing premise is that the hand-tuned Eulerian repulsive force with coefficients a = 55 and b = 3.5 prevents coalescence completely without otherwise changing the droplet dynamics, interface deformation, or the stress budget; the paper itself reports that this force contributes about 10% of the total stress, so a perturbing effect large enough to alter the comparison cannot be excluded from the reported evidence.
Editorial extensions
If this is right
- Emulsion constitutive curves must be reported with coalescence efficiency controlled; two emulsions with identical composition, drop size, and shear rate can have qualitatively different viscosity curves depending on whether merging is allowed.
- Analytical viscosity formulas that assume positive curvature, such as Pal's equation and the Eilers formula, should only be applied to non-coalescing emulsions; the paper shows they fail when coalescence is active.
- Suppressing coalescence makes an emulsion behave like a suspension of deformable particles: the viscosity data collapse onto the Eilers formula using an effective volume fraction based on droplet deformation, with about 6% error at the highest concentration.
- The ratio of first to second normal stress difference stays roughly constant with volume fraction, and the first normal stress difference changes sign at the lowest viscosity ratio—an experimentally checkable rheological fingerprint.
- Flow topology changes with coalescence: merging droplets leave the matrix in nearly pure shear, whereas non-coalescing emulsions develop extensional flow in the gaps between droplets.
Reading between the lines
- If the mechanism is generic, coalescence efficiency should be treated as a rheological state variable as important as volume fraction and capillary number; any processing change that alters drop-size distribution—surfactant dose, pre-shear history, compatibilizer—will slide the system along the family of curves between the two limits computed here.
- The reported 10% stress contribution of the artificial repulsive force is a directly testable artifact check: a second, force-free way to inhibit coalescence (for instance an interface with strong Marangoni resistance) should reproduce the positive curvature and the same stress-budget decomposition.
- The area-reduction explanation predicts a quantitative scaling: the drop in effective viscosity should track the loss of total interfacial area, so simultaneous measurement of drop size and shear stress during coalescence would give a direct test of the mechanism.
- For viscosity ratios below one, coalescence can make an emulsion less viscous than its continuous phase; engineering controlled coalescence could therefore be a purely microstructural route to viscosity reduction without changing chemistry.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. De Vita et al. present volume-of-fluid simulations of shear-driven emulsions in a Couette cell at vanishing Reynolds number, scanning volume fraction φ=0.00164–0.3, viscosity ratio λ=0.01–1, and capillary number Ca=0.05–0.2. Two branches are compared: one in which droplets coalesce naturally, and one in which an Eulerian repulsive force (Eq. 3.1) suppresses coalescence, intended to mimic surfactant-stabilized emulsions. The main finding is that the effective viscosity versus φ curve has negative curvature when coalescence is allowed, with a maximum near φ≈0.2 and values below unity for λ<1, and positive curvature when coalescence is suppressed, resembling suspensions of deformable particles. The authors attribute the difference to the reduction of interfacial area upon merging, which lowers the interfacial-tension contribution to the shear stress. Normal stress differences, flow-topology statistics, and droplet size distributions are also reported.
Significance. If the central comparison were clean, the paper would be a valuable contribution: it gives a mechanistic explanation of the experimentally observed negative curvature (Caserta & Guido 2012) and demonstrates that neglecting coalescence in emulsion simulations changes the qualitative rheology. Strengths include the careful stress-budget derivation (Eqs. 4.2–4.8), grid-resolution checks, and benchmarks against Einstein, Taylor, and the Pal model in the dilute regime. The coalescing branch is compared to experiments with reasonable qualitative agreement. The principal caveat is that the non-coalescing control is not a passive variation: the collision force contributes directly to the stress balance and alters the wall-normal droplet distribution, so the strength of the causal claim (that coalescence causes the curvature reversal) depends on the inertness of that force model.
major comments (3)
- [§4.1, Eqs. (4.1), (4.8), Fig. 8] The effective viscosity is computed from the wall shear stress via Eq. (4.1), while the bulk stress balance in Eq. (4.8) includes the collision-force integral C. Because the collision force is a numerical body force whose streamwise component is about 10% of the total shear stress at φ=0.3, the reported µ_e for the non-coalescing branch is not strictly a material property of the emulsion. The authors should provide a convergence study in the force strength (e.g., increasing a and b beyond 55 and 3.5) to demonstrate that the suspension viscosity is independent of the force amplitude as long as coalescence is prevented. Without this, the positive curvature of the non-coalescing curve could be partly an artifact of the model force.
- [§4.1, Fig. 9(right)] The wall-normal profile of the average volume fraction is substantially modified by the collision force, with a marked increase of <φ> near the walls. Since µ_e is evaluated from wall quantities, this force-induced segregation directly contributes to the measured viscosity increase, independently of the interfacial-area mechanism emphasized in the paper. The statement in §4 that the dispersed phase is 'approximately homogeneous' is difficult to reconcile with the changes in near-wall <φ> shown in Fig. 9(right). The authors should quantify the wall-region contribution to µ_e (e.g., by computing µ_e from the bulk stress instead and comparing) or demonstrate that the results are insensitive to the redistribution.
- [§3, Fig. 4] The calibration of the collision force coefficients a and b (Eq. 3.1) is done by trial and error on a single two-drop collision, and the only reported sensitivity test compares two force amplitudes on the post-collision trajectory. This does not establish that the force is dynamically inert in the dense many-drop suspensions: the force acts in a shell of thickness 3Δ around every interface, can alter the local flow in the gaps, and changes the deformation of droplets. A sensitivity test of the suspension rheology and of the stress budget to a and b in a dense case (e.g., φ=0.3) is needed to support the claim that the difference between the two branches is due solely to the presence or absence of coalescence.
minor comments (5)
- [Eq. (2.3)] The advection equation for the VoF function appears to contain a typo: ∂(u_j H)/∂x_j should likely be ∂(u_j T)/∂x_j, or the right-hand side should involve H ∂u_j/∂x_j, since H is the color function and T is the cell average. Please check the notation.
- [Fig. 3 caption] The rows in the caption are mislabeled: the third row is described twice as '(d)-(e)-(f)', and the fourth row should be '(g)-(h)-(i)' or '(j)-(k)-(l)' accordingly.
- [Author affiliations] The affiliation contains a typo: 'Stockhom' should be 'Stockholm'.
- [Fig. 14 caption] The black diamonds representing the experiments of Caserta & Guido are not identified in the caption; please state explicitly what they denote.
- [§4.2, Fig. 13] The statement that the ratio of N1 and N2 is 'almost constant' with volume fraction is supported by the figure, but for λ=1 the scatter is non-negligible; consider providing a quantitative measure such as a best-fit slope with confidence intervals.
Circularity Check
No significant circularity: the rheology is simulated, not fitted; the non-coalescing control uses a tuned collision force, but the viscosity result is emergent.
full rationale
The paper's central claim—that coalescence produces negative curvature in the effective-viscosity-versus-volume-fraction curve while suppressing coalescence produces positive curvature—is obtained from direct numerical simulation, not from a fitted or self-referential construction. The coalescing branch is validated against external benchmarks: the Einstein/Taylor dilute limits, Pal's formula in the dilute regime, and the Caserta & Guido experiments. The non-coalescing branch is generated by adding the Eulerian repulsive force of Eq. (3.1), with coefficients a=55 and b=3.5 chosen by trial and error only to prevent merging; the viscosity curve itself is not a fit target, so the curvature sign is an emergent simulation outcome. The stress budget in Eq. (4.8) explicitly includes the collision-force contribution C, and the paper reports it is about 10% of the total stress; this is a possible physical confound in comparing the two branches, but it is not circularity because the measured viscosity is not defined in terms of the force amplitude and the force was not tuned to reproduce any rheological curve. The self-citations (Rosti et al. 2019 for the VoF code, Rosti & Brandt 2018 and Rosti et al. 2018 for the deformable-particle Eilers comparison) are used as methodological references and external comparisons, not as premises that entail the coalescence result. No equation is defined in terms of the claimed finding, no fitted parameter is renamed as a prediction, and no uniqueness theorem is imported from the authors' prior work. The paper is therefore self-contained against external benchmarks and its central derivation chain is not circular, though the non-coalescing control deserves scrutiny as a modeling artifact rather than a circularity issue.
Assumptions & free parameters
free parameters (1)
- collision force coefficients a, b =
a=55, b=3.5
assumptions (5)
- standard math Incompressible Navier-Stokes equations accurately describe the two-fluid flow.
- domain assumption Reynolds number Re=0.1 is sufficiently small for inertia to be negligible.
- domain assumption The VoF/MTHINC method with the given grid accurately tracks the interface and resolves the stresses.
- ad hoc to paper The Eulerian repulsive force (equation 3.1) models the effect of surfactants and prevents coalescence without otherwise altering the flow physics.
- domain assumption The simulation domain (16x16x10 r) is large enough to represent the constitutive behavior of the emulsion, i.e. finite-size and wall effects are negligible for the reported viscosities.
invented entities (1)
-
Eulerian collision force Fc
Cite this review
Pith. "Pith review of On the effect of coalescence on the rheology of emulsions." pith.science (2026). https://pith.science/paper/5GCVSP5M
@misc{pith2026190808383,
author = {Pith},
title = {Pith review of: On the effect of coalescence on the rheology of emulsions},
year = {2026},
howpublished = {\url{https://pith.science/paper/5GCVSP5M}},
note = {Machine review of arXiv:1908.08383}
}
read the original abstract
We present a numerical study of the rheology of a two-fluid emulsion in dilute and semidilute conditions. The analysis is performed for different capillary numbers, volume fraction and viscosity ratio under the assumption of negligible inertia and zero buoyancy force. The effective viscosity of the system increases for low values of the volume fraction and decreases for higher values, with a maximum for about 20 % concentration of the disperse phase. When the dispersed fluid has lower viscosity, the normalised effective viscosity becomes smaller than 1 for high enough volume fractions. To single out the effect of droplet coalescence on the rheology of the emulsion we introduce an Eulerian force which prevents merging, effectively modelling the presence of surfactants in the system. When the coalescence is inhibited the effective viscosity is always greater than 1 and the curvature of the function representing the emulsion effective viscosity vs. the volume fraction becomes positive, resembling the behaviour of suspensions of deformable particles. The reduction of the effective viscosity in the presence of coalescence is associated to the reduction of the total surface of the disperse phase when the droplets merge, which leads to a reduction of the interface tension contribution to the total shear stress. The probability density function of the flow topology parameter shows that the flow is mostly a shear flow in the matrix phase, with regions of extensional flow when the coalescence is prohibited. The flow in the disperse phase, instead, always shows rotational components. The first normal stress difference is positive whereas the second normal difference is negative, with their ratio being constant with the volume fraction. Our results clearly show that the coalescence efficiency strongly affects the system rheology and neglecting droplet merging can lead to erroneous predictions.
Figures
Figures from the paper (12 more)
Reference graph
Works this paper leans on
-
[1]
Albadawi, A., Donoghue, D. B., Robinson, A. J., Murray, D. B. & Delaur ´e, Y. M C 2013 Influence of surface tension implementation in Volume o f Fluid and coupled Volume of Fluid with Level Set methods for bubble growth and detachm ent. International Journal of Multiphase Flow 53, 11–28
work page 2013
-
[2]
Journal of Fluid Mechanics 41, 545—-&
Batchelor, G K 1970 Stress System In A Suspension of Force-free Particles. Journal of Fluid Mechanics 41, 545—-&
work page 1970
-
[3]
Journal of Fluid Mechanics 56 (NOV28), 375—-&
Batchelor, G K & Green, J T 1972 Hydrodynamic Interaction of 2 Small Freely-moving Spheres In A Linear Flow Field. Journal of Fluid Mechanics 56 (NOV28), 375—-&
work page 1972
-
[4]
Bolotnov, Igor A., Jansen, Kenneth E., Drew, Donald A., Ober ai, Assad A., Lahey, Richard T. & Podowski, Michael Z. 2011 Detached direct numerical simulations of turbulent two-phase bubbly channel flow. International Journal of Multiphase Flow 37 (6), 647–659
work page 2011
-
[5]
Brackbill, J. U., Kothe, D. B. & Zemach, C. 1992 A continuum method for modeling surface tension. Journal of Computational Physics 100 (2), 335–354, arXiv: fld.1
work page 1992
-
[6]
Caserta, Sergio & Guido, Stefano 2012 Vorticity banding in biphasic polymer blends. Langmuir 28 (47), 16254–16262
work page 2012
-
[7]
Caserta, S., Sabetta, L., Simeone, M. & Guido, S. 2005 Shear-induced coalescence in aqueous biopolymer mixtures. Chemical Engineering Science 60 (4), 1019–1027. Chesters, a. K 1991 The modelling of coalescence process in fluid-liquid di spersions: A Review of Current Understanding. Chemical Engineering Research and Design 69 (4), 259–270
work page 2005
-
[8]
Choi, Seung Jung & Schow alter, W. R. 1975 Rheological properties of nondilute suspensions of deformable particles. Physics of Fluids 18 (4), 420–427. De Vita, F., Rosti, M E, Izbassarov, D, Du, L, Tammisola, O, Ho rmozi, S & Brandt, L 2018 Elastoviscoplastic flows in porous media. Journal of Non-Newtonian Fluid Mechanics 258 (April), 10–21
work page 1975
Show all 35 references
-
[9]
1941 Die Viskositat von Emulsionen hochviskoser Stoffe als F unktion der Konzentration
Eilers, Von H. 1941 Die Viskositat von Emulsionen hochviskoser Stoffe als F unktion der Konzentration. Kolloidal-Z 97, 313–321
1941
-
[10]
1906 Eine neue Bestimmung der Molekuldimensionen
Einstein, A. 1906 Eine neue Bestimmung der Molekuldimensionen. Annalen der Physik 19, 289–306
1906
-
[11]
1911 Berichtingung zu meiner Arbeit:”Eine neue Bestimmung der Molekuldimensionen”
Einstein, A. 1911 Berichtingung zu meiner Arbeit:”Eine neue Bestimmung der Molekuldimensionen”. Annalen der Physik 34, 591–592
1911
-
[12]
Advances in Colloid and Interface Science 220, 78–91
Foudazi, Reza, Qavi, Sahar, Masalova, Irina & Malkin, Alexa nder Ya 2015 Physical chemistry of highly concentrated emulsions. Advances in Colloid and Interface Science 220, 78–91
2015
-
[13]
Journal of Fluid Mechanics 852 (P1), 1–73, arXiv: 1801.03805
Guazzelli, Elisabeth & Pouliquen, Olivier 2018 Rheology of dense granular suspensions under extensional flow. Journal of Fluid Mechanics 852 (P1), 1–73, arXiv: 1801.03805
2018 arXiv
-
[14]
Journal of Fluid Mechanics 357, 1–20
Guido, S & Simeone, M 1998 Binary collision of drops in simple shear flow. Journal of Fluid Mechanics 357, 1–20
1998
-
[15]
Hudson, S. D. 2003 Wall migration and shear-induced diffusion of fluid drop lets in emulsions. Physics of Fluids 15 (5), 1106–1113
2003
-
[16]
Journal of Computational Physics 231 (5), 2328–2358
Ii, Satoshi, Sugiyama, Kazuyasu, Takeuchi, Shintaro, Taka gi, Shu & Matsumoto, Yoichiro 2012 An interface capturing method with a continuous functi on : The THINC method with multi-dimensional reconstruction. Journal of Computational Physics 231 (5), 2328–2358
2012
-
[17]
Computers and Fluids 49 (1), 150–165
Khatri, Shilpa & Tornberg, Anna-karin 2011 Computers & Fluids A numerical method for two phase flows with insoluble surfactants. Computers and Fluids 49 (1), 150–165
2011
-
[18]
& Morris, Jeffrey F
Kulkarni, P andurang M. & Morris, Jeffrey F. 2008 Suspension properties at finite Reynolds number from simulated shear flow. Physics of Fluids 20 (4)
2008
-
[19]
Loewenberg, M & Hinch, E J 1996 Numerical simulations of a concentrated emulsion in shear flow. J. Fluid Mech. 321, 395–419
1996
-
[20]
& Denn, Morto n M
Mari, Romain, Seto, Ryohei, Morris, Jeffrey F. & Denn, Morto n M. 2014 Shear thickening, frictionless and frictional rheologies in non -Brownian suspensions 1693 (2014), arXiv: 1403.6793
2014 arXiv
-
[21]
& Ishikaw a, T
Matsunaga, D., Imai, Y., Yamaguchi, T. & Ishikaw a, T. 2015 Rheology of a dense suspension of spherical capsules under simple shear flow. Journal of Fluid Mechanics 786, 110–127. P al, Rajinder 2002 Novel shear modulus equations for concentrated emulsi ons of two 22 F. De Vita,...
2015
-
[22]
Physical Review Letters 111 (9), 2–6
Picano, Francesco, Breugem, Wim P aul, Mitra, Dhrubaditya & Brandt, Luca 2013 Shear thickening in non-Brownian suspensions: An excluded volume effect. Physical Review Letters 111 (9), 2–6
2013
-
[23]
& Brandt, Luca 2018 Suspensions of deformable particles in a Couette flow
Rosti, Marco E. & Brandt, Luca 2018 Suspensions of deformable particles in a Couette flow. Journal of Non-Newtonian Fluid Mechanics 262 (C), 3–11
2018
-
[24]
Physical Review Fluids 3 (1), 012301
Rosti, Marco E., Brandt, Luca & Mitra, Dhrubaditya 2018 Rheology of suspensions of viscoelastic spheres: Deformability as an effective volume fraction. Physical Review Fluids 3 (1), 012301
2018
-
[25]
Acta Mechanica 230 (2), 667–682
Rosti, Marco Edoardo, De Vita, Francesco & Brandt, Luca 2019 Numerical simulations of emulsions in shear flows. Acta Mechanica 230 (2), 667–682
2019
-
[26]
Journal of Computational Physics 163 (1), 51–67
Russo, Giovanni & Smereka, Peter 2000 A Remark on Computing Distance Functions. Journal of Computational Physics 163 (1), 51–67. Schow alter, W. R., Chaffey, C. E. & Brenner, How ard 1968 Rheological behavior of a dilute emulsion. Journal of Colloid And Interface Science 26 (2...
2000
-
[27]
& Ivanov, Alexey O
Sega, Marcello, Sbragaglia, Mauro, Kantorovich, Sofia S. & Ivanov, Alexey O. 2013 Mesoscale structures at complex fluid-fluid interfaces : A novel lattice Boltzmann/molecular dynamics coupling. Soft Matter 9 (42), 10092–10107
2013
-
[28]
Shardt, Orest, Derksen, J. J. & Mitra, Sushanta K. 2013 Simulations of droplet coalescence in simple shear flow. Langmuir 29 (21), 6201–6212
2013
-
[29]
Physical Review Letters 97 (5), 1–4
Guido, Stefano 2006 Drop deformation in microconfined shear flow. Physical Review Letters 97 (5), 1–4
2006
-
[30]
Journal of Computational Physics 376, 1292–1311
Soligo, Giovanni, Roccon, Alessio & Soldati, Alfredo 2019 Coalescence of surfactant- laden drops by Phase Field Method. Journal of Computational Physics 376, 1292–1311
2019
-
[31]
Journal of Fluid Mechanics 805, 494–522
Srivastava, Priyesh, Malipeddi, Abhilash Reddy & Sarkar, K ausik 2016 Steady shear rheology of a viscous emulsion in the presence of finite inert ia at moderate volume fractions: Sign reversal of normal stress differences. Journal of Fluid Mechanics 805, 494–522
2016
-
[32]
Taylor, G. I. 1932 The Viscosity of a Fluid Containig Small Drops of Anothe r Fluid. Proc. R. Soc. Lond. A 138, 41–48. V ananroye, Anja, V an Puyvelde, Peter & Moldenaers, P aula 2006 Structure development in confined polymer blends: Steady-state shear flow and relaxation. Langm...
1932
-
[33]
Xia, By Younan, Gates, Byron, Yin, Yadong & Lu, Yu 2000 Spillover enhanced hydrogen storage in Pt doped Supporting information of MOF graphene o xide composite produced via an impregnation method of supporting pp. 693–713
2000
-
[34]
2016 Numerical simulations of the rheology of suspensions of rigid spheres at low volume fr action in a viscoelastic fluid under shear
Yang, Mengfei, Krishnan, Sreenath & Shaqfeh, Eric S.G. 2016 Numerical simulations of the rheology of suspensions of rigid spheres at low volume fr action in a viscoelastic fluid under shear. Journal of Non-Newtonian Fluid Mechanics 233, 181–197
2016
-
[35]
1993 The flow of ordered and random suspensions of two- dimensional drops in a channel
Zhou, Hua & Pozrikidis, C. 1993 The flow of ordered and random suspensions of two- dimensional drops in a channel. Journal of Fluid Mechanics 255, 103–127
1993
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.