pith. sign in

arxiv: 2410.04629 · v1 · submitted 2024-10-06 · ❄️ cond-mat.soft · physics.flu-dyn

Solute dispersion in pre-turbulent confined active nematics

Tom\'as Alvim , Margarida M. Telo da Gama , Rodrigo C. V. Coelho This is my paper

Pith reviewed 2026-05-23 19:52 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.flu-dyn
keywords active nematicssolute dispersionpre-turbulent regimesnematohydrodynamicsvelocity momentsconfined channelsoscillatory flowdancing flow
0
0 comments X

The pith

Longitudinal solute dispersion in pre-turbulent confined active nematics is set by velocity second moments and scales with activity.

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

Simulations of nematohydrodynamics examine solute spreading inside narrow channels filled with active nematic fluid. Two pre-turbulent regimes are compared: oscillatory flow that carries net mass flux and dancing flow that does not. In both cases the longitudinal dispersion coefficient is fixed by the second moments of the longitudinal and transverse velocity components and grows linearly with activity. This common dependence indicates the same underlying mechanism operates despite the different large-scale flow patterns. The resulting dispersion reaches up to ten times the molecular value in dancing flow and is also boosted in oscillatory flow.

Core claim

The longitudinal dispersion of solutes, in both the oscillatory and dancing pre-turbulent regimes, is determined by the second moments of the longitudinal and transverse components of the velocity field and proportional to the activity, suggesting a common mechanism for dispersion in these distinct flows.

What carries the argument

Second moments of the longitudinal and transverse velocity components, which fix the dispersion coefficient and enforce its linear scaling with activity.

If this is right

  • Dispersion coefficient increases up to one order of magnitude relative to molecular diffusion in the dancing flow regime.
  • Dispersion is also enhanced in the oscillatory flow regime.
  • A single velocity-moment mechanism governs dispersion across both distinct pre-turbulent flows.
  • The scaling opens routes for controlled mixing in biological environments and microfluidic devices.

Where Pith is reading between the lines

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

  • Direct imaging of velocity fluctuations could test the moment scaling without needing to track individual solute particles.
  • The same moment-based relation may apply to other confined active-matter systems whose velocity statistics resemble these regimes.
  • Varying channel aspect ratio or boundary conditions could be used to tune the dispersion coefficient predictably.

Load-bearing premise

The nematohydrodynamics model and chosen simulation parameters faithfully reproduce the pre-turbulent regimes without unaccounted physical effects that would alter the reported velocity-moment scaling.

What would settle it

A measurement showing that solute dispersion in either regime fails to scale linearly with activity or deviates from the value predicted by the velocity second moments would falsify the claimed relation.

Figures

Figures reproduced from arXiv: 2410.04629 by Margarida M. Telo da Gama, Rodrigo C. V. Coelho, Tom\'as Alvim.

Figure 1
Figure 1. Figure 1: Comparison of the nematic director field, flow velocity, and vorticity in oscillatory (a, c, e) and dancing flow (b, d, f), in a subsection [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of a solute drop in oscillatory a) and dancing flow b) regimes, in the section of the channel where the solute drop is [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Time evolution of the mean squared displacement (MSD) [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: Effective diffusivity in oscillatory and dancing flow, as a [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: Average velocity of the flow a) and effective diffusivity [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: Dynamics of tracers and their trajectories. In a), the log-log plot shows the MSD of tracers, calculated from Eq. 16. The initial [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
read the original abstract

We investigate the dispersion of solutes in active nematic fluids confined in narrow channels based on simulations of nematohydrodynamics. The study focuses on two pre-turbulent regimes: oscillatory flow, with net mass flux, and dancing flow, without net flux. Our findings reveal that the longitudinal dispersion of solutes, in both regimes, is determined by the second moments of the longitudinal and transverse components of the velocity field and proportional to the activity, suggesting a common mechanism for dispersion in these distinct flows. Finally, we found that the dispersion coefficient may increase up to one order of magnitude with respect to molecular diffusion in the dancing flow regime and that it is also enhanced in oscillatory flows, with potential for applications in natural biological environments and engineered devices.

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

2 major / 2 minor

Summary. The manuscript uses direct numerical simulations of nematohydrodynamics equations to study solute dispersion in narrow-channel confined active nematics, focusing on two pre-turbulent regimes: oscillatory flow (with net mass flux) and dancing flow (without net flux). It reports that longitudinal dispersion in both regimes is set by the second moments of the longitudinal and transverse velocity components, scales proportionally with activity, and shares a common mechanism; dispersion is enhanced by up to an order of magnitude relative to molecular diffusion in the dancing regime.

Significance. If the reported scaling holds after verification, the work provides a simulation-grounded link between velocity statistics and effective transport in active fluids, with potential relevance to biological mixing and microfluidic applications. The direct solution of the governing equations and comparison across distinct flow regimes are strengths, though the proportionality claim requires explicit checks on activity dependence of correlation times.

major comments (2)
  1. [Abstract and §4] Abstract and §4 (results on scaling): the claim that dispersion 'is determined by the second moments ... and proportional to the activity' with a 'common mechanism' requires showing that the velocity autocorrelation time remains independent of activity (or identical across regimes). Without this, D_eff ~ <u_x'^2> * tau does not reduce to proportionality with second moments alone when activity varies.
  2. [§3.2] §3.2 (simulation parameters): the chosen activity strengths, channel aspect ratios, and boundary conditions are stated to produce the pre-turbulent regimes, but no systematic check is reported on whether unaccounted effects (e.g., inertia or boundary-induced secondary flows) alter the velocity-moment scaling at the highest activities examined.
minor comments (2)
  1. [Figure 5] Figure 5 caption: the error bars on the dispersion coefficient vs. activity are not defined (standard deviation across runs?); clarify the number of independent realizations used.
  2. [§2] Notation: the definition of the effective dispersion coefficient D_eff should be stated explicitly in the methods (e.g., via mean-squared displacement or flux moments) rather than only referenced to prior literature.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful review and constructive suggestions. We address each major comment below and will revise the manuscript accordingly to strengthen the claims.

read point-by-point responses

  1. Referee: [Abstract and §4] Abstract and §4 (results on scaling): the claim that dispersion 'is determined by the second moments ... and proportional to the activity' with a 'common mechanism' requires showing that the velocity autocorrelation time remains independent of activity (or identical across regimes). Without this, D_eff ~ <u_x'^2> * tau does not reduce to proportionality with second moments alone when activity varies.

    Authors: We agree that explicit verification of activity-independent correlation times is needed to rigorously support the proportionality claim. In our simulations, the velocity autocorrelation times remain approximately constant across the studied activity range in both regimes. We will add this verification, including a plot of correlation time versus activity, to §4 and the supplementary material in the revised manuscript. This addition will also reinforce the evidence for a common mechanism. revision: yes

  2. Referee: [§3.2] §3.2 (simulation parameters): the chosen activity strengths, channel aspect ratios, and boundary conditions are stated to produce the pre-turbulent regimes, but no systematic check is reported on whether unaccounted effects (e.g., inertia or boundary-induced secondary flows) alter the velocity-moment scaling at the highest activities examined.

    Authors: The parameters were chosen following established literature to remain in the overdamped, low-Reynolds-number regime characteristic of active nematics, with no observed secondary flows. To address the concern directly, we will include in the revised §3.2 a brief analysis confirming that the Reynolds number stays ≪1 across the activity range and that velocity-moment scaling is unaffected by inertial effects. revision: yes

Circularity Check

0 steps flagged

No circularity; claims from direct numerical simulation of nematohydrodynamics

full rationale

The paper obtains its results via direct numerical solution of the nematohydrodynamics equations for two pre-turbulent regimes. The reported relation between longitudinal dispersion, velocity second moments, and activity is an observed outcome of the simulations rather than a quantity defined by construction, fitted from the target data, or reduced via self-citation. No load-bearing self-citations, uniqueness theorems, or ansatzes appear in the abstract or context. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete. The work relies on the standard nematohydrodynamics framework whose parameters (activity, elastic constants, viscosity ratios) are not enumerated here.

free parameters (1)
  • activity strength
    Activity is the control parameter that sets the flow regime and is expected to be chosen or scanned to reach the pre-turbulent states described.
axioms (1)
  • domain assumption nematohydrodynamics equations govern the dynamics of the confined active nematic
    The study is explicitly based on simulations of nematohydrodynamics.

pith-pipeline@v0.9.0 · 5663 in / 1284 out tokens · 32551 ms · 2026-05-23T19:52:33.964454+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Active nematic pumps

    cond-mat.soft 2024-07 unverdicted novelty 7.0

    Triangular inclusions in active nematic gels stabilize self-pumping flows for wall-free microfluidic cargo transport and mixing.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    Ahrens, B

    J. Ahrens, B. Geveci, and C. Law, ParaView, Kitware, Inc. (2024), available at https://www.paraview.org

    work page 2024

  2. [2]

    I. C. Engelhardt, D. Patko, Y . Liu, M. Mimault, G. De Las Heras Martinez, T. S. George, M. MacDonald, M. Ptashnyk, T. Sukhodub, N. R. Stanley-Wall, N. Holden, T. J. Daniell, and L. X. Dupuy, The ISME Journal 16, 2337 (2022)

    work page 2022

  3. [3]

    Kim and K

    M. Kim and K. Breuer, Small 4, 111 (2008)

    work page 2008

  4. [4]

    Marenduzzo, E

    D. Marenduzzo, E. Orlandini, M. E. Cates, and J. M. Yeomans, Physical Review E 76, 031921 (2007)

    work page 2007

  5. [5]

    LIN, J.-L

    Z. LIN, J.-L. THIFFEAULT, and S. CHILDRESS, Journal of Fluid Mechanics 669, 167–177 (2011)

    work page 2011

  6. [6]

    Wu and A

    X.-L. Wu and A. Libchaber, Physical Review Letters 84, 3017 (2000)

    work page 2000

  7. [7]

    K. C. Leptos, J. S. Guasto, J. P. Gollub, A. I. Pesci, and R. E. Goldstein, Physical Review Letters 103, 198103 (2009)

    work page 2009

  8. [8]

    Aditi Simha and S

    R. Aditi Simha and S. Ramaswamy, Physical Review Letters 89, 058101 (2002)

    work page 2002

  9. [9]

    Wolke, E

    U. Wolke, E. A. Jezuit, and J. R. Priess, Development134, 2227 (2007)

    work page 2007

  10. [10]

    Alert, J

    R. Alert, J. Casademunt, and J. F. Joanny, Annual Review of Condensed Matter Physics 13, 10.1146/annurev-conmatphys- 082321-035957 (2022)

    work page doi:10.1146/annurev-conmatphys- 2022

  11. [11]

    Chandragiri, A

    S. Chandragiri, A. Doostmohammadi, J. M. Yeomans, and S. P. Thampi, Soft Matter 15, 1597 (2019)

    work page 2019

  12. [12]

    Hardo ¨uin, R

    J. Hardo ¨uin, R. Hughes, A. Doostmohammadi, J. Laurent, T. Lopez-Leon, J. M. Yeomans, J. Ign´es-Mullol, and F. Sagu´es, Communications Physics 2, 121 (2019)

    work page 2019

  13. [13]

    Wioland, E

    H. Wioland, E. Lushi, and R. E. Goldstein, New Journal of Physics 18, 075002 (2016)

    work page 2016

  14. [14]

    R. J. Henshaw, O. G. Martin, and J. S. Guasto, Physical Review Fluids 8, 023101 (2023)

    work page 2023

  15. [15]

    S. P. Thampi, Current Opinion in Colloid & Interface Science 61, 101613 (2022)

    work page 2022

  16. [16]

    T. N. Shendruk, A. Doostmohammadi, K. Thijssen, and J. M. Yeomans, Soft Matter 13, 3853 (2017)

    work page 2017

  17. [17]

    A. J. Tan, E. Roberts, S. A. Smith, U. A. Olvera, J. Arteaga, S. Fortini, K. A. Mitchell, and L. S. Hirst, Nature Physics 15, 1033 (2019)

    work page 2019

  18. [18]

    F. L. Memarian, D. Hammar, M. M. H. Sabbir, M. Elias, K. A. Mitchell, and L. S. Hirst, Physical Review Letters 132, 228301 (2024)

    work page 2024

  19. [19]

    J. P. Hernandez-Ortiz, C. G. Stoltz, and M. D. Graham, Physical Review Letters 95, 204501 (2005)

    work page 2005

  20. [20]

    P. T. Underhill, J. P. Hernandez-Ortiz, and M. D. Graham, Phys- ical Review Letters 100, 248101 (2008)

    work page 2008

  21. [21]

    A. S. Mikhailov, Y . Koyano, and H. Kitahata, Journal of the Physical Society of Japan 86, 101013 (2017)

    work page 2017

  22. [22]

    Dehkharghani, N

    A. Dehkharghani, N. Waisbord, J. Dunkel, and J. S. Guasto, Proceedings of the National Academy of Sciences 116, 11119 (2019)

    work page 2019

  23. [23]

    T. E. Bate, M. E. Varney, E. H. Taylor, J. H. Dickie, C.-C. Chueh, M. M. Norton, and K.-T. Wu, Nature Communications 13, 6573 (2022)

    work page 2022

  24. [24]

    K. J. Modica, A. K. Omar, and S. C. Takatori, Soft Matter 19, 1890 (2023)

    work page 2023

  25. [25]

    Ge and G

    Z. Ge and G. J. Elfring, Hydrodynamic diffusion in apolar ac- tive suspensions (2024), arXiv:2311.09448 [cond-mat.soft]

    work page arXiv 2024

  26. [26]

    Basu, J.-F

    A. Basu, J.-F. Joanny, F. J ¨ulicher, and J. Prost, New Journal of Physics 14, 115001 (2012)

    work page 2012

  27. [27]

    Das, Langmuir 40, 18773 (2024)

    S. Das, Langmuir 40, 18773 (2024)

    work page 2024

  28. [28]

    Aris, Proceedings of the Royal Society of London

    R. Aris, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 235, 67 (1956). 8

    work page 1956

  29. [29]

    W. Xi, T. B. Saw, D. Delacour, C. T. Lim, and B. Ladoux, Na- ture Reviews Materials 4, 23 (2018)

    work page 2018

  30. [30]

    Gupta, B

    M. Gupta, B. R. Sarangi, J. Deschamps, Y . Nematbakhsh, A. Callan-Jones, F. Margadant, R.-M. M`ege, C. T. Lim, R. V oi- turiez, and B. Ladoux, Nature Communications6, 7525 (2015)

    work page 2015

  31. [31]

    S. P. Thampi, A. Doostmohammadi, R. Golestanian, and J. M. Yeomans, EPL (Europhysics Letters) 112, 28004 (2015)

    work page 2015

  32. [32]

    Santhosh, M

    S. Santhosh, M. R. Nejad, A. Doostmohammadi, J. M. Yeo- mans, and S. P. Thampi, Journal of Statistical Physics 180, 699 (2020)

    work page 2020

  33. [33]

    Kr ¨uger, H

    T. Kr ¨uger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, and E. M. Viggen, The Lattice Boltzmann method: principles and practice, Graduate texts in physics (Springer International Publishing, Cham, 2017)

    work page 2017

  34. [34]

    Samui, J

    A. Samui, J. M. Yeomans, and S. P. Thampi, Soft Matter 17, 10640 (2021)

    work page 2021

  35. [35]

    Gulati, S

    P. Gulati, S. Shankar, and M. C. Marchetti, Frontiers in Physics 10, 948415 (2022)

    work page 2022

  36. [36]

    Li, D.-Q

    Z.-Y . Li, D.-Q. Zhang, and B. Li, Physical Review Research 3, 023253 (2021)

    work page 2021

  37. [37]

    C. G. Wagner, R. H. Pallock, J. S. Park, M. M. Norton, and P. Grover, Physical Review Fluids8, 124401 (2023)

    work page 2023

  38. [38]

    Doostmohammadi, T

    A. Doostmohammadi, T. N. Shendruk, K. Thijssen, and J. M. Yeomans, Nature Communications 8, 15326 (2017)

    work page 2017

  39. [39]

    R. C. V . Coelho, N. A. M. Ara ´ujo, and M. M. Telo da Gama, Soft Matter 15, 6819 (2019)

    work page 2019

  40. [40]

    Chandrakar, M

    P. Chandrakar, M. Varghese, S. Aghvami, A. Baskaran, Z. Dogic, and G. Duclos, Physical Review Letters 125, 257801 (2020)

    work page 2020

  41. [41]

    E. J. Hemingway, P. Mishra, M. C. Marchetti, and S. M. Field- ing, Soft Matter 12, 7943 (2016)

    work page 2016

  42. [42]

    V oituriez, J

    R. V oituriez, J. F. Joanny, and J. Prost, Europhysics Letters (EPL) 70, 404 (2005)

    work page 2005

  43. [43]

    M. L. Parks and L. A. Romero, Mathematical and Computer Modelling 46, 699 (2007)

    work page 2007

  44. [44]

    Active nematic pumps

    I. V ´elez-Ceron, R. C. V . Coelho, P. Guillamat, M. T. da Gama, F. Sagu´es, and J. Ign ´es-Mullol, Active nematic pumps (2024), arXiv:2407.09960 [cond-mat.soft]

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  45. [45]

    R. C. Coelho, N. A. Ara ´ujo, and M. M. Telo da Gama, Soft Matter 18, 7642 (2022)

    work page 2022