REVIEW 4 major objections 5 minor 38 references
How interfacial tension enhances drag in turbulent Taylor-Couette flow with neutrally buoyant and equally viscous droplets
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Droplet drag in Taylor–Couette flow traced to interfacial tension on stretched droplet fore heads.
desk verdict Solid simulation study with an exact flux decomposition and a plausible, clearly-stated fore-head/rear-end mechanism; the high-volume-fraction causal chain rests on unverified droplet statistics, so review but with careful attention to the coalescence caveats. 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 tool is the angular velocity flux decomposition $J^{\omega} = J_T^{\omega}(r) + J_V^{\omega}(r) + J_\sigma^{\omega}(r) = \mathrm{const.}$, which splits the conserved flux into turbulent, viscous, and interfacial-tension transport. The interfacial term is $J_\sigma^{\omega}(r) = -\int_{r_i}^r \langle r^2 f_\sigma^{\theta} \rangle\,dr$, so it directly converts the azimuthal component of interfacial tension into a drag contribution. The causal mechanism is the fore-head/rear-end asymmetry: near the inner cylinder, shear stretches droplets, and interfacial tension on the fore head acts against the flow while that on the rear end acts with it; the fore-head effect dominates, producing a net hindering force that slows the continuous phase and raises wall viscous stress.
What would settle it
A direct test would be to run the same geometry and parameters with interfacial tension switched off ($\sigma = 0$) while keeping droplets deformable by another mechanism; if drag enhancement persists, the claim that interfacial tension is the source fails. Alternatively, tracking individual droplets near the inner cylinder and measuring the azimuthal interfacial tension summed over fore heads should give a net against-flow value that grows with volume fraction; observing zero or positive fore-head forces would falsify the mechanism.
Extended reading notes
Core claim
The central claim is that the droplet-induced drag enhancement in this system originates from the contribution of interfacial tension, specifically the hindering of the continuous phase by the fore head of stretched droplets. In a Taylor-Couette cell at Re=5200 with neutrally buoyant, equally viscous droplets at volume fractions 0-40%, the torque required to drive the inner cylinder increases with the volume fraction, in line with experiment. Decomposing the conserved angular velocity flux into turbulent, viscous, and interfacial-tension parts shows that the interfacial part grows monotonically with volume fraction, while the radial averages of the other two stay essentially unchanged. Within the viscous and buffer layers near the inner cylinder, shear deforms droplets in the streamwise direction so that the rear end lags behind the fore head; interfacial tension then acts against the flow on the fore head and with the flow on the rear end, with the fore-head effect dominant. This braking effect lowers the angular velocity near the wall, disrupts high-speed streaks, increases the occurrence of low-speed streaks, and thereby raises viscous stress and drag.
Load-bearing premise
The load-bearing assumption is that the modified multi-marker VOF method with marker-group volume fraction 5% produces droplet sizes, deformation, and coalescence behavior representative of the real experiments at all volume fractions up to 40%; the paper itself notes that numerical coalescence remains and that no resolution test was run at 40%.
Editorial extensions
If this is right
- If the central claim is correct, models of turbulent two-phase drag must include interfacial tension in the near-wall region, and omitting it will underpredict drag at moderate volume fractions.
- The linear wall law $u^+ = y^+$ breaks down in the viscous sublayer of two-phase flow, because interfacial tension acts as an additional near-wall stress rather than a simple roughness effect.
- Drag enhancement should grow with droplet volume fraction at fixed Reynolds and Weber numbers, consistent with experiments up to about 30% volume fraction, with the deviation at 40% attributed to residual numerical coalescence.
- The mechanism explains why deformable droplets enhance drag less than rigid particles: deformability lets droplets adjust shape and exert a weaker hindering effect on the surrounding flow.
- In industrial Euler-Lagrange type predictions, the boundary-layer interfacial tension contribution must be added to correctly capture the torque increase reported here.
Reading between the lines
- A testable extension is to vary the Weber number at fixed Reynolds number: the fore-head braking effect should weaken as interfacial tension becomes weaker relative to inertia, giving a shear-thinning trend in effective viscosity.
- If the mechanism is general, the same fore-head braking should appear in other wall-bounded shear flows with deformable droplets, including turbulent channel and pipe flows, not only Taylor-Couette geometry.
- Comparing droplets with progressively higher surface tension against rigid-particle suspensions could sharpen the claim that deformability is what caps the drag enhancement.
- Direct experimental measurement of droplet interface curvature near the wall could reveal the sign of the azimuthal interfacial tension on fore heads and rear ends, providing a non-simulation check on the proposed causal chain.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies turbulent Taylor–Couette flow laden with neutrally buoyant, equally viscous droplets at volume fractions up to 40%, using a modified multi-marker volume-of-fluid method with the dispersed phase split into marker groups of 5% volume fraction each. The main claim is that the observed drag enhancement originates from interfacial tension: an exact decomposition of the angular velocity flux shows that the interfacial-tension contribution J_sigma grows with phi while the turbulent and viscous contributions do not, and near-wall diagnostics indicate that stretched droplets near the inner cylinder exert a fore-head interfacial-tension force opposing the flow, which slows the continuous phase, disrupts high-speed streaks, and raises wall viscous stress. The paper reports agreement with previous experiments on global torque and presents supporting statistics: radial profiles of J_sigma, joint PDFs of viscous stress and angular velocity, phase-resolved angular velocity PDFs, and per-droplet fore-head/rear-end interfacial tension PDFs.
Significance. If the proposed mechanism is correct, it offers a concrete physical picture for droplet-induced drag enhancement in turbulent emulsions, connecting interfacial tension in the viscous sublayer to streak modification and wall stress. The work combines an exact flux budget with multiple independent diagnostics and includes open-source code, a comparison against experiments, and a spurious-current assessment, which are strengths. The significance is moderate: it is a mechanism-focused numerical study at a single Reynolds and Weber number, and its broader reach depends on the fidelity of the simulated droplet population at high volume fractions.
major comments (4)
- [Appendix B, Fig. 10] The manuscript states that the modified multi-marker VOF method 'does not completely resolve the issue of numerical coalescence' and that at phi=40% the simulated droplet-size PDF deviates from experimental results. Because the causal mechanism in §3 is explicitly deformation- and size-dependent (fore-head versus rear-end interfacial tension on stretched droplets), the mechanism is least verified exactly in the regime where the drag enhancement is largest. The standard-VOF comparison in Fig. 11 shows quantitatively different droplet and continuous phase velocity PDFs at phi=40%, confirming that the numerical method affects the droplet statistics on which the mechanism rests. The authors should either provide evidence that the remaining numerical coalescence does not change the sign or magnitude of the fore-head PDF in Fig. 5(b), or explicitly restrict the mechanism claim to volume fractions where the droplet size distribution is validated.
- [Appendix A] The resolution test is conducted only for phi=0 and phi=10%, and the paper acknowledges that a resolution test at phi=40% is impractical. The proposed mechanism depends on interface curvature, droplet deformation, and interfacial tension in the viscous sublayer; the Kolmogorov-scale grid estimate does not address interface-resolution requirements such as curvature accuracy and parasitic currents at high phi. A resolution study at an intermediate high volume fraction, or at least a quantitative estimate of how the fore-head interfacial tension statistics change with grid refinement, is needed to support the deformation-based mechanism at the volume fractions where it is invoked.
- [§2, choice of phi_mvf] The marker-group volume fraction phi_mvf=5% is a free parameter that controls the maximum droplet size and the degree of numerical coalescence, and no sensitivity study with respect to this parameter is reported. Since the fore-head/rear-end mechanism depends on droplet size and deformation, and since the choice is justified only by reproduction of the global drag from experiments, the robustness of the mechanism to different phi_mvf values should be demonstrated to rule out a method-induced artifact.
- [§3, discussion near Fig. 5] The fore-head/rear-end division by the droplet center of mass is acknowledged by the authors as 'somewhat idealized.' The proposed causal chain relies on this division being physically meaningful, but the paper does not test whether the result is sensitive to the definition of the dividing surface, for example by using a different split criterion or by weighting droplets by their deformation. Such a sensitivity check would strengthen the causal interpretation beyond the exact flux budget.
minor comments (5)
- [§1 and Appendix B] The phrase 'fore head' is written inconsistently with the standard 'forehead' and appears without hyphen in most places; please unify the spelling.
- [§2, notation near Eq. (3.1)] The operator 'Í' used for radial averaging is nonstandard and not defined in the main text; please define it explicitly at first use.
- [Fig. 5(b) caption] The normalization of the PDFs by the absolute value of the total azimuthal interfacial tension experienced by all droplets in the range y+<24 should be stated more clearly in the caption, since the PDFs are not conventional probability densities.
- [§2, text after Fig. 2] The sentence 'The datasets agree well with each other at phi ⩽ 30%, with a minor deviation observed at phi=40%' should be supported by an explicit error metric, as the deviation at 40% is later used to qualify the method's fidelity.
- [Appendix C] The statement that |u|_max/u_i < 0.02 is acceptable would be more informative with a comparison to the typical turbulent velocity fluctuations in the simulations, rather than only to the inner-cylinder velocity.
Circularity Check
No significant circularity: the interfacial-tension mechanism is diagnosed from the simulated surface-tension force, not imposed by the calibration of phi_mvf.
full rationale
The only potentially self-referential element is the choice of marker-group volume fraction phi_mvf=5%, which the paper states 'faithfully reproduces the global drag of the system reported in our previous experimental study (Yi et al. 2021)' in Section 2. This is a calibration of one numerical parameter to a global benchmark, not a fit of the mechanism. The central claim that drag enhancement originates from interfacial tension is obtained from the exact angular-velocity-flux decomposition in Eq. (3.1), where J_sigma(r) is computed directly from the simulated azimuthal surface-tension force via J_sigma(r) = -∫_{r_i}^{r} <r^2 f_sigma^theta> dr. Nothing in that budget is set equal to the experimental drag increment; the turbulent and viscous contributions are independently measured and shown not to grow with phi. The causal story is supported by separate diagnostics: the fore-head/rear-end PDFs of summed interfacial tension (Fig. 5b), the negative joint PDF between angular velocity and viscous stress (Fig. 3c), the droplet/continuous phase velocity PDFs (Fig. 4), and the standard-VOF comparison in Appendix B. The acknowledged limitations—numerical coalescence at 30-40% (Appendix B) and the absence of a resolution test at phi=40% (Appendix A)—are validity and verification concerns about the droplet population, not circular reductions of the conclusion to its inputs. The paper therefore has no significant circularity.
Assumptions & free parameters
free parameters (2)
- Marker-group volume fraction phi_mvf =
5%
- Domain size parameters =
azimuthal angle pi/3, axial aspect ratio 2pi/3
assumptions (5)
- standard math The Navier-Stokes equations with a continuum surface force model govern the two-phase flow (Eqs. 2.1-2.4).
- domain assumption A 6-fold rotational symmetry and axial aspect ratio 2pi/3 adequately represent the full Taylor-Couette annulus.
- domain assumption The droplet population reaches a statistically stationary state with breakup-coalescence balance.
- ad hoc to paper The division of a droplet into fore head and rear end at its center of mass is physically meaningful for the mechanism.
- standard math The angular velocity flux decomposition (Eq. 3.1) is exact and the radial average of the three contributions is the appropriate diagnostic.
Cite this review
Pith. "Pith review of How interfacial tension enhances drag in turbulent Taylor-Couette flow with neutrally buoyant and equally viscous droplets." pith.science (2026). https://pith.science/paper/SCVDO63Q
@misc{pith2026241113115,
author = {Pith},
title = {Pith review of: How interfacial tension enhances drag in turbulent Taylor-Couette flow with neutrally buoyant and equally viscous droplets},
year = {2026},
howpublished = {\url{https://pith.science/paper/SCVDO63Q}},
note = {Machine review of arXiv:2411.13115}
}
abstract
The presence of dispersed-phase droplets can result in a notable increase in the system's drag. However, our understanding of the mechanism underlying this phenomenon remains limited. In this study, we use three-dimensional direct numerical simulations with a modified multi-marker volume-of-fluid method to investigate liquid-liquid two-phase turbulence in a Taylor-Couette geometry. The dispersed phase has the same density and viscosity as the continuous phase. The Reynolds number $Re\equiv r_i\omega_i d/\nu$ is fixed at 5200, the volume fraction of the dispersed phase is up to $40\%$, and the Weber number $We\equiv \rho u^2_\tau d/\sigma$ is around 8. It is found that the increase in the system's drag originates from the contribution of interfacial tension. Specifically, droplets experience significant deformation and stretching in the streamwise direction due to shear near the inner cylinder. Consequently, the rear end of the droplets lags behind the fore head. This causes opposing interfacial tension effects on the fore head and rear end of the droplets. For the fore head of the droplets, the effect of interfacial tension appears to act against the flow direction. For the rear end, the effect appears to act in the flow direction. The increase in the system's drag is primarily attributed to the effect of interfacial tension on the fore head of the droplets which leads to the hindering effect of the droplets on the surrounding continuous phase. This hindering effect disrupts the formation of high-speed streaks, favoring the formation of low-speed ones, which are generally associated with higher viscous stress and drag of the system. This study provides new insights into the mechanism of drag enhancement reported in our previous experiments.
Figures
Figures from the paper (9 more)
Reference graph
Works this paper leans on
-
[1]
Assen, Martin PA, Ng, Chong Shen, Will, Jelle B, Stevens, Richard JAM, Lohse, Detlef & Verzicco, Roberto 2022 Strong alignment of prolate ellipsoids in Taylor–Couette flow. J. Fluid Mech. 935, A7
work page 2022
-
[2]
Huisman, Sander G 2021 Catastrophic phase inversion in high-Reynolds-number turbulent Taylor– Couette flow. Phys. Rev. Lett. 126 (6), 064501. Balc´azar, N ´estor, Lehmkuhl, Oriol, Rigola, Joaquim & Oliva, Assensi 2015 A multiple marker level-set method for simulation of deformable fluid particles. Int. J. Multiphase Flow 74, 125–142
work page 2021
-
[3]
Barcilon, AJMFR, Brindley, J, Lessen, M & Mobbs, FR 1979 Marginal instability in Taylor–Couette flows at a very high Taylor number.J. Fluid Mech. 94 (3), 453–463
work page 1979
-
[4]
Brauckmann, Hannes J & Eckhardt, Bruno 2013 Direct numerical simulations of local and global torque in Taylor–Couette flow up to Re= 30 000.J. Fluid Mech. 718, 398–427
work page 2013
-
[5]
Chen, Songtao, Zhao, Weiwen & Wan, Decheng 2022 Turbulent structures and characteristics of flows past a vertical surface-piercing finite circular cylinder. Phys. Fluids 34 (1), 015115
work page 2022
-
[6]
Coyajee, Emil & Boersma, Bendiks Jan 2009 Numerical simulation of drop impact on a liquid–liquid interface with a multiple marker front-capturing method. J. Comput. Phys. 228 (12), 4444–4467
work page 2009
-
[7]
Crialesi-Esposito, Marco, Rosti, Marco Edoardo, Chibbaro, Sergio & Brandt, Luca 2022 Modulation of homogeneous and isotropic turbulence in emulsions. J. Fluid Mech. 940, A19
work page 2022
-
[8]
Dong, S 2007 Direct numerical simulation of turbulent Taylor–Couette flow.J. Fluid Mech. 587, 373–393
work page 2007
Show all 38 references
-
[9]
Eckhardt, Bruno, Grossmann, Siegfried & Lohse, Detlef 2007 Torque scaling in turbulent Taylor– Couette flow between independently rotating cylinders. J. Fluid Mech. 581, 221–250
2007
-
[10]
Elghobashi, Said 2019 Direct numerical simulation of turbulent flows laden with droplets or bubbles. Annu. Rev. Fluid Mech. 51, 217–244
2019
-
[11]
Loeben, Germany: Holzmann CFD
Holzmann, Tobias 2016 Mathematics, numerics, derivations and openfoam ®. Loeben, Germany: Holzmann CFD
2016
-
[12]
Sun, Chao 2013 Logarithmic boundary layers in strong Taylor–Couette turbulence.Phys. Rev. Lett. 110 (26), 264501
2013
-
[13]
Innocenti, Alessio, Jaccod, Alice, Popinet, St ´ephane & Chibbaro, Sergio 2021 Direct numerical simulation of bubble-induced turbulence. J. Fluid Mech. 918, A23
2021
-
[14]
Energy Fuels 26 (7), 4017–4026
Kilpatrick, Peter K 2012 Water-in-crude oil emulsion stabilization: review and unanswered questions. Energy Fuels 26 (7), 4017–4026
2012
-
[15]
Maxey, Martin R 2017 Droplets in turbulence: a new perspective. J. Fluid Mech. 816, 1–4
2017
-
[16]
Mcclements, David Julian 2007 Critical review of techniques and methodologies for characterization of emulsion stability. Crit. Rev. Food Sci. Nutr.47 (7), 611–649
2007
-
[17]
Mukherjee, Siddhartha, Safdari, Arman, Shardt, Orest, Kenjere ˇs, Sa ˇsa & Van den Akker, Harry EA 2019 Droplet–turbulence interactions and quasi-equilibrium dynamics in turbulent emulsions. J. Fluid Mech. 878, 221–276
2019
-
[18]
Nemati, Hassan, Breugem, Wim-Paul, Kwakkel, Marcel & Boersma, Bendiks Jan 2021 Direct numerical simulation of turbulent bubbly down flow using an efficient CLSVOF method. Int. J. Multiphase Flow 135, 103500
2021
-
[19]
Ni, Rui 2024 Deformation and breakup of bubbles and drops in turbulence. Annu. Rev. Fluid Mech. 56, 319–347
2024
-
[20]
Fluid Mech.719, 14–46
Ostilla, Rodolfo, Stevens, Richard JAM, Grossmann, Siegfried, Verzicco, Roberto & Lohse, Detlef 2013 Optimal Taylor–Couette flow: direct numerical simulations.J. Fluid Mech.719, 14–46
2013
-
[21]
Picano, Francesco, Breugem, Wim-Paul & Brandt, Luca 2015 Turbulent channel flow of dense suspensions of neutrally buoyant spheres. J. Fluid Mech. 764, 463–487
2015
-
[22]
Popinet, St´ephane 2018 Numerical models of surface tension. Annu. Rev. Fluid Mech. 50 (1), 49–75
2018
-
[23]
Rosti, Marco E, Brandt, Luca & Mitra, Dhrubaditya 2018 Rheology of suspensions of viscoelastic spheres: Deformability as an effective volume fraction. Phys. Rev. Fluids 3 (1), 012301
2018
-
[24]
Rosti, Marco E, Ge, Zhouyang, Jain, Suhas S, Dodd, Michael S & Brandt, Luca 2019 Droplets in homogeneous shear turbulence. J. Fluid Mech. 876, 962–984
2019
-
[25]
PhD thesis, Imperial College London (University of London)
Rusche, Henrik 2003 Computational fluid dynamics of dispersed two-phase flows at high phase fractions. PhD thesis, Imperial College London (University of London)
2003
-
[26]
Scarbolo, Luca & Soldati, Alfredo 2015 Wall drag modification by large deformable droplets in turbulent channel flow. Comput. Fluids 113, 87–92
2015
-
[27]
Scheufler, Henning & Roenby, Johan 2019 Accurate and efficient surface reconstruction from volume fraction data on general meshes. J. Comput. Phys. 383, 1–23
2019
-
[28]
Spandan, Vamsi, Verzicco, Roberto & Lohse, Detlef 2018 Physical mechanisms governing drag reduction in turbulent Taylor–Couette flow with finite-size deformable bubbles. J. Fluid Mech. 849, R3. 20 J Su, Y Zhang, C Wang, L Yi, F Xu, Y Fan, J Wang, and C Sun
2018
-
[29]
Spernath, Aviram & Aserin, Abraham 2006 Microemulsions as carriers for drugs and nutraceuticals. Adv. Colloid Interface Sci. 128, 47–64
2006
-
[30]
Su, Jinghong, Wang, Cheng, Zhang, Yi-bao, Xu, Fan, Wang, Junwu & Sun, Chao 2024a Turbulence modulation in liquid-liquid two-phase Taylor-Couette turbulence, arXiv: 2407.00992
-
[31]
Thomas, S, Esmaeeli, A & Tryggvason, G 2010 Multiscale computations of thin films in multiphase flows. Int. J. Multiphase Flow 36 (1), 71–77. de Vita, Francesco, Rosti, Marco Edoardo, Caserta, Sergio & Brandt, Luca2019 On the effect of coalescence on the rheology of emulsions....
2010
-
[32]
Fluid Mech.952, A39
Wang, Cheng, Yi, Lei, Jiang, Linfeng & Sun, Chao 2022 Turbulence drag modulation by dispersed droplets in Taylor–Couette flow: the effects of the dispersed phase viscosity.J. Fluid Mech.952, A39
2022
-
[33]
Wang, Guiquan, Abbas, Micheline & Climent, ´Eric 2017 Modulation of large-scale structures by neutrally buoyant and inertial finite-size particles in turbulent Couette flow. Phys. Rev. Fluids 2 (8), 084302
2017
-
[34]
Wang, Lijuan, Li, Xuefeng, Zhang, Gaoyong, Dong, Jinfeng & Eastoe, Julian 2007 Oil-in-water nanoemulsions for pesticide formulations. J. Colloid Interface Sci. 314 (1), 230–235
2007
-
[35]
Xu, Fan, Su, Jinghong, Lan, Bin, Zhao, Peng, He, Yurong, Sun, Chao & Wang, Junwu 2023 Direct numerical simulation of Taylor–Couette flow with vertical asymmetric rough walls. J. Fluid Mech. 975, A30
2023
-
[36]
Xu, Fan, Zhao, Peng, Sun, Chao, He, Yurong & Wang, Junwu 2022 Direct numerical simulation of Taylor-Couette flow: Regime-dependent role of axial walls.Chem. Eng. Sci. 263, 118075
2022
-
[37]
Yi, Lei, Toschi, Federico & Sun, Chao 2021 Global and local statistics in turbulent emulsions. J. Fluid Mech. 912, A13
2021
-
[38]
Yi, Lei, Wang, Cheng, van Vuren, Thomas, Lohse, Detlef, Risso, Fr´ed´eric, Toschi, Federico & Sun, Chao 2022 Physical mechanisms for droplet size and effective viscosity asymmetries in turbulent emulsions. J. Fluid Mech. 951, A39
2022
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.