Phase separation dynamics and active turbulence in a binary fluid mixture
Pith reviewed 2026-05-18 01:54 UTC · model grok-4.3
The pith
Activity in a binary active-passive fluid mixture arrests phase domain coarsening, producing a finite length scale that shrinks as activity rises.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the two-fluid model, active stress enhances velocity and vorticity differences between phases and increased active concentration strengthens inter-fluid coupling through drag. Activity amplifies turbulent fluctuations yet arrests domain coarsening, producing a finite characteristic length scale that decreases with rising activity. The active parameter, tumbling parameter, and Frank elastic constant each influence the resulting flow scale.
What carries the argument
Two distinct momentum equations for the active and passive phases, coupled solely by viscous drag, together with Cahn-Hilliard dynamics and Beris-Edwards nematohydrodynamics.
If this is right
- Domain coarsening halts at a finite scale whose size falls as activity is increased.
- Velocity and vorticity contrasts between the active and passive phases grow with active stress.
- Inter-phase coupling strengthens at higher active concentrations.
- The characteristic scale varies with the active parameter, tumbling parameter, and Frank elastic constant.
Where Pith is reading between the lines
- The same drag-mediated arrest could appear in bacterial suspensions inside polymeric fluids where phase separation must remain controlled.
- Varying the viscosity contrast between phases in experiments would test whether relative motion alone sets the observed length scale.
- Similar finite-scale stabilization might occur in other active-passive emulsions if the coupling mechanism is replaced by different interfacial forces.
Load-bearing premise
All relevant momentum transfer and relative motion between the two phases is captured by viscous drag alone, without additional interfacial forces or non-viscous couplings.
What would settle it
A simulation or experiment in which separated domains continue to coarsen without bound even at high activity, or in which the measured length scale shows no systematic decrease with increasing activity strength.
read the original abstract
Active matter, encompassing natural systems, converts surrounding energy to sustain autonomous motion, exhibiting unique non-equilibrium behaviors such as active turbulence and phase separation. In this study, we develop a continuum two-fluid model for a binary mixture of an active nematic and a passive Newtonian fluid, coupling Cahn-Hilliard dynamics for phase separation with Beris-Edwards nematohydrodynamics and two distinct momentum equations connected by viscous drag. A phase field-based lattice Boltzmann method is used to investigate the existence of active turbulence and phase separation in the binary mixture. We find that active stress enhances velocity and vorticity differences between phases, and that increased active concentration promotes stronger inter-fluid coupling. Activity not only amplifies turbulent fluctuations but also arrests domain coarsening, leading to a finite characteristic length scale that decreases with increasing activity. Key parameters, like active parameter, tumbling parameter and Frank elastic constant, affect the characteristic scale of flow. These results highlight the role of relative motion and drag-mediated momentum transfer in active binary mixtures, providing a framework for studying systems such as bacterial suspensions in polymeric fluids or active emulsions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a two-fluid continuum model for a binary mixture consisting of an active nematic and a passive Newtonian fluid. It couples Cahn-Hilliard dynamics for phase separation to Beris-Edwards nematohydrodynamics, employing two distinct momentum equations linked only by a viscous drag term. Phase-field lattice Boltzmann simulations are used to show that activity increases velocity and vorticity differences between phases, strengthens inter-fluid coupling at higher active concentrations, amplifies turbulent fluctuations, and arrests domain coarsening, producing a finite characteristic length scale that decreases with increasing activity. The dependence of this scale on the active parameter, tumbling parameter, and Frank elastic constant is also examined.
Significance. If the central claim holds, the work supplies a useful modeling framework for active binary mixtures (e.g., bacterial suspensions in polymeric fluids or active emulsions) by isolating the effects of relative motion and drag-mediated momentum transfer on both turbulence and the saturation of phase separation.
major comments (1)
- [model description / equations section] Model description (paragraph following the abstract and the equations section): the two momentum equations are stated to be connected solely by viscous drag. It is not clear whether the chemical-potential gradient force from the Cahn-Hilliard free energy or the nematic elastic stresses at the interface are included as body forces in each phase’s momentum balance. If these standard interfacial contributions are omitted, the reported arrest of coarsening to a finite scale may be an artifact of the reduced force set rather than a generic consequence of activity-induced relative motion; this point is load-bearing for the central claim.
minor comments (2)
- [abstract] Abstract: the statement that “increased active concentration promotes stronger inter-fluid coupling” would benefit from a quantitative metric (e.g., a correlation function or drag-induced velocity alignment measure) rather than a qualitative description.
- [results] Results section: the characteristic length scale is reported to decrease with activity, but no error bars, ensemble averages, or sensitivity to grid resolution / time-step are provided; these should be added to establish robustness.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for identifying this important point of clarification in the model description. We address the concern directly below and will revise the manuscript accordingly to improve clarity without altering the underlying model.
read point-by-point responses
-
Referee: [model description / equations section] Model description (paragraph following the abstract and the equations section): the two momentum equations are stated to be connected solely by viscous drag. It is not clear whether the chemical-potential gradient force from the Cahn-Hilliard free energy or the nematic elastic stresses at the interface are included as body forces in each phase’s momentum balance. If these standard interfacial contributions are omitted, the reported arrest of coarsening to a finite scale may be an artifact of the reduced force set rather than a generic consequence of activity-induced relative motion; this point is load-bearing for the central claim.
Authors: We appreciate the referee drawing attention to the precise formulation of the interfacial forces. In the two-fluid model, the momentum equations for the active nematic and passive phases include the chemical-potential gradient force derived from the Cahn-Hilliard free energy as a body force acting on each phase (weighted by the local phase concentration). The nematic elastic stresses from the Beris-Edwards description enter the divergence of the stress tensor in the active-phase momentum equation, with their interfacial contributions transmitted through the phase-field coupling and the viscous drag term. The phrase 'connected solely by viscous drag' in the current text was intended to emphasize the inter-phase momentum exchange mechanism but inadvertently omitted explicit mention of these standard contributions. The arrest of coarsening arises from activity-driven velocity differences amplified by the drag coupling, in the presence of the full set of interfacial forces. We will revise the model description and equations section to state the body forces explicitly and to remove any ambiguity about the force balance. revision: yes
Circularity Check
No circularity: results obtained from direct numerical solution of explicitly stated model equations
full rationale
The paper constructs a two-fluid continuum model by coupling Cahn-Hilliard dynamics, Beris-Edwards nematohydrodynamics, and two momentum equations linked by a viscous drag term, then integrates the system numerically via a phase-field lattice Boltzmann method. All reported behaviors, including the arrest of domain coarsening and the activity-dependent characteristic length scale, are presented as direct outcomes of these simulations rather than as fitted quantities or quantities derived by re-expressing the input equations. No self-definitional steps, fitted-input predictions, or load-bearing self-citations appear in the derivation chain; the central claims remain independent of the numerical results themselves.
Axiom & Free-Parameter Ledger
free parameters (3)
- active parameter
- tumbling parameter
- Frank elastic constant
axioms (2)
- domain assumption Continuum two-fluid description with viscous drag coupling is adequate for the length and time scales of interest.
- domain assumption Beris-Edwards nematohydrodynamics plus Cahn-Hilliard phase separation capture the essential physics of the active nematic component.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
two distinct momentum equations connected by viscous drag... Fdrag = γ φ (1-φ)(u2 - u1)
-
IndisputableMonolith/Foundation/BranchSelection.leanbranch_selection unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
activity arrests domain coarsening, leading to a finite characteristic length scale that decreases with increasing activity
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
F. Yang, S. Liu, H. J. Lee, R. Phillips and M. Thomson, Nature Materials, 2025, 24, 615--625
work page 2025
- [2]
- [3]
-
[4]
M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao and R. A. Simha, Rev. Mod. Phys., 2013, 85, 1143--1189
work page 2013
-
[5]
Ramaswamy, Annual Review of Condensed Matter Physics, 2010, 1, 323--345
S. Ramaswamy, Annual Review of Condensed Matter Physics, 2010, 1, 323--345
work page 2010
-
[6]
H. H. Wensink, J. Dunkel, S. Heidenreich, K. Drescher, R. E. Goldstein, H. Löwen and J. M. Yeomans, Proceedings of the National Academy of Sciences, 2012, 109, 14308--14313
work page 2012
- [7]
-
[8]
C. Dombrowski, L. Cisneros, S. Chatkaew, R. E. Goldstein and J. O. Kessler, Phys. Rev. Lett., 2004, 93, 098103
work page 2004
- [9]
-
[10]
S. Mukherjee, R. K. Singh, M. James and S. S. Ray, Nature Physics, 2023, 19, 891--897
work page 2023
- [11]
-
[12]
R. Alert, J.-F. Joanny and J. Casademunt, Nature Physics, 2020, 16, 682--688
work page 2020
-
[13]
A. Doostmohammadi, J. Ignes-Mullol, J. Yeomans and F. Sagués, Nature Communications, 2018, 9,
work page 2018
- [14]
- [15]
-
[16]
H. Li, X. qing Shi, M. Huang, X. Chen, M. Xiao, C. Liu, H. Chaté and H. P. Zhang, Proceedings of the National Academy of Sciences, 2019, 116, 777--785
work page 2019
-
[17]
A. Doostmohammadi and J. Yeomans, The European Physical Journal Special Topics, 2019, 227, 2401--2411
work page 2019
-
[18]
S. P. Thampi, Current Opinion in Colloid and Interface Science, 2022, 61, 101613
work page 2022
- [19]
-
[20]
R. Assante, D. Corbett, D. Marenduzzo and A. Morozov, Soft Matter, 2023, 19, 189--198
work page 2023
-
[21]
S. Bhattacharyya and J. M. Yeomans, Phys. Rev. Lett., 2023, 130, 238201
work page 2023
-
[22]
C. P. Brangwynne, C. R. Eckmann, D. S. Courson, A. Rybarska, C. Hoege, J. Gharakhani, F. Jülicher and A. A. Hyman, Science, 2009, 324, 1729--1732
work page 2009
-
[23]
C. P. Brangwynne, T. J. Mitchison and A. A. Hyman, Proceedings of the National Academy of Sciences, 2011, 108, 4334--4339
work page 2011
-
[24]
M. E. Cates and J. Tailleur, Annual Review of Condensed Matter Physics, 2015, 6, 219--244
work page 2015
- [25]
-
[26]
R. Wittkowski, A. Tiribocchi, J. Stenhammar, R. J. Allen, D. Marenduzzo and M. E. Cates, Nature Communications, 2014, 5, 4351
work page 2014
- [27]
-
[28]
J. Stenhammar, A. Tiribocchi, R. J. Allen, D. Marenduzzo and M. E. Cates, Phys. Rev. Lett., 2013, 111, 145702
work page 2013
-
[29]
H.-S. Kuan, W. P\"onisch, F. J\"ulicher and V. Zaburdaev, Phys. Rev. Lett., 2021, 126, 018102
work page 2021
-
[30]
A. Tiribocchi, R. Wittkowski, D. Marenduzzo and M. E. Cates, Phys. Rev. Lett., 2015, 115, 188302
work page 2015
- [31]
- [32]
-
[33]
S. Saha, J. Agudo-Canalejo and R. Golestanian, Phys. Rev. X, 2020, 10, 041009
work page 2020
- [34]
-
[35]
M. L. Blow, S. P. Thampi and J. M. Yeomans, Phys. Rev. Lett., 2014, 113, 248303
work page 2014
-
[36]
P. Yue, J. FENG, C. Liu and J. Shen, Journal of Fluid Mechanics, 2004, 515, 293 -- 317
work page 2004
-
[37]
D. Marenduzzo, E. Orlandini, M. Cates and J. Yeomans, Physical review. E, Statistical, nonlinear, and soft matter physics, 2007, 76, 031921
work page 2007
-
[38]
C. Denniston, E. Orlandini and J. M. Yeomans, Phys. Rev. E, 2001, 63, 056702
work page 2001
-
[39]
S. P. Thampi, R. Golestanian and J. M. Yeomans, Europhysics Letters, 2014, 105, 18001
work page 2014
-
[40]
A. Doostmohammadi, M. F. Adamer, S. P. Thampi and J. M. Yeomans, Nature Communications, 2016, 7, 10566
work page 2016
-
[41]
A. Dadvand, M. Bagheri, N. Samkhaniani, H. Marschall and M. Wörner, Physics of Fluids, 2021, 33, 053311
work page 2021
-
[42]
S. P. Thampi, A. Doostmohammadi, R. Golestanian and J. M. Yeomans, Europhysics Letters, 2015, 112, 28004
work page 2015
-
[43]
R. Saghatchi, D. C. Kolukisa and M. Yildiz, International Journal for Numerical Methods in Engineering, 2023, 124, 159--182
work page 2023
-
[44]
M. E. Cates and E. Tjhung, Journal of Fluid Mechanics, 2018, 836, P1
work page 2018
-
[45]
R. Mueller and A. Doostmohammadi, Phase field models of active matter, 2021, https://arxiv.org/abs/2102.05557
-
[46]
S. P. Thampi, R. Golestanian and J. M. Yeomans, Philosophical transactions. Series A, Mathematical, physical, and engineering sciences, 2014, 372, 20130366
work page 2014
-
[47]
A. N. Kolmogorov, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences, 1991, 434, 9--13 mcitethebibliography rsc-articletemplate-softmatter.tex0000664000000000000000000011557415101340733016444 0ustar rootroot [twoside,twocolumn,9pt] article extsizes [super,sort&compress,comma] natbib [version=3] mhchem [left=1.5cm, rig...
-
[48]
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION format.date year empty "" 'year if FU...
work page 2059
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.