pith. sign in

arxiv: 2604.11034 · v1 · submitted 2026-04-13 · ⚛️ physics.flu-dyn

Patchy Polymeric Scalar Turbulence

Rahul K. Singh , Marco E. Rosti This is my paper

Pith reviewed 2026-05-10 16:19 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords polymeric turbulencescalar mixingpassive scalarsSchmidt numberturbulent mixingnon-Newtonian fluidspatchy fluctuations
0
0 comments X

The pith

Polymeric turbulence forms small patchy scalar fluctuations instead of large islands, reducing mixing efficiency at moderate diffusivities.

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

This paper examines mixing of passive scalar fields in turbulent polymeric flows versus Newtonian turbulence at varying Schmidt numbers. It shows that polymeric scalar turbulence develops small interspersed patches of strong fluctuations rather than the large contiguous islands seen in Newtonian scalar turbulence. These patches occupy a larger volume fraction but have weaker gradients and lower average flux across boundaries. The patches also exhibit smoother, more space-filling boundaries along with stronger yet less intermittent spatial changes. The overall result indicates that polymeric turbulence mixes scalars less efficiently than Newtonian turbulence at small to moderate Schmidt numbers.

Core claim

In polymeric scalar turbulence the scalar field organizes into small interspersed patches of strong but less intermittent fluctuations that occupy a larger volume fraction than the large islands in Newtonian scalar turbulence. This structure produces smaller scalar gradients and therefore smaller average flux across boundaries, even as spatial changes of the scalar remain stronger overall with slower self-similar growth and reduced intermittency revealed by kurtosis of scalar differences.

What carries the argument

The patchy organization of scalar fluctuations in polymeric turbulence versus contiguous island fronts in Newtonian turbulence, measured through volume fractions of strong fluctuations, box-counting dimensions of boundaries, and kurtosis of scalar increments.

If this is right

  • Mixing times for dissolved substances in polymer-laden turbulent flows would increase at moderate diffusivities.
  • Transport models for non-Newtonian turbulence must include the effect of patch formation on net scalar flux.
  • Engineering designs relying on polymer additives for flow control should anticipate slower scalar homogenization.
  • The distinction between patches and islands provides a diagnostic for identifying polymeric effects in scalar fields.

Where Pith is reading between the lines

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

  • The same patch mechanism may operate in other viscoelastic fluids and could explain broader reductions in mixing for non-Newtonian turbulence.
  • At sufficiently high Schmidt numbers the patchy regime might cross over to Newtonian-like island behavior.
  • Direct visualization of scalar fields in laboratory polymer solutions could confirm whether the reduced flux persists outside simulations.

Load-bearing premise

The numerical simulations with the selected constitutive model and parameters for polymeric flows faithfully represent physical behavior without artifacts that would artificially produce patchiness instead of islands.

What would settle it

Laboratory measurements of scalar gradient statistics or average flux in a real polymeric turbulent flow at moderate Schmidt number that show mixing rates equal to or higher than in Newtonian turbulence.

Figures

Figures reproduced from arXiv: 2604.11034 by Marco E. Rosti, Rahul K. Singh.

Figure 1
Figure 1. Figure 1: Passive scalar fluctuations about the mean for Sc = 1 in polymeric turbulence, with elasticity increasing from left [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Main panels The probability distribution functions (pdfs) of the passive scalar fluctuations about the mean δϕ. A PT background results in stronger fluctuations about the mean, especially at an optimal De = 1. Insets The standard deviation ||δϕ|| vs De. The first glimpses into the nature of polymeric scalar turbulence (PST) is shown in [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Main panels The fraction of points S constituting the boundary of regions within which δϕ ≥ ϕ0. Insets Volume fraction V of all the points in the bulk of the same regions. Larger S and V in PT means regions with large scalar concentration occupy a maximal volume. Thus, mixing is less efficient in PT, and more so at De = 1. trations in PST. The width of the distributions, given by their variance ∥δϕ∥ = p ⟨δ… view at source ↗
Figure 4
Figure 4. Figure 4: The pdfs of normalised scalar gradients ∇^ϕ(x) = L|∇ϕ(x)| measured at the boundaries enclosing the regions described by δϕ ⩾ 6 for (a) Sc = 0.3, (b) 1.0, (c) 3.0. boxes N(r) of size r that cover an iso-level set where N(r) ∼ r −D [40, 53–55]. For instance, a space filling and well mixed scalar in three dimensions has D = 3 for the isolevels δϕ = 0, i.e N(r) ∼ r −3 . A D < 3 im￾plies that fluctuations are c… view at source ↗
Figure 5
Figure 5. Figure 5: Number N(r) of cubes with side r required to cover the boundaries of regions where δϕ ⩾ 6 are maximal at De = 1 for all Sc. A steeper fall-off of N(r) shows patch boundaries are more space filling in PST. The three sets of curves for each Sc have been moved vertically for visual clarity with only Sc = 1 showing the correct number of boxes. themselves, suggesting that there exists an upper bound to mixing i… view at source ↗
Figure 6
Figure 6. Figure 6: Log-log plots of second-order scalar structure function [PITH_FULL_IMAGE:figures/full_fig_p005_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The scaling exponents ζp for scalar structure func￾tions plotted in [PITH_FULL_IMAGE:figures/full_fig_p005_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Flatness F(r) as a function of r for different De and Sc as the non-dimensionalised ratio of S ϕ 6 (r) and S ϕ 2 (r). Intermittency remains subduedx in PST. patches and the fluctuations scale with a spectrum of ex￾ponents in space. A quantification on how important are these deviations is given by the (kurtosis) flatness F(r) of the δrϕ distributions F(r) = ⟨(δrϕ) 6 ⟩ ⟨(δrϕ) 2 ⟩ 3 = S ϕ 6 (r) [S ϕ 2 (r)] 3… view at source ↗
read the original abstract

Turbulent polymeric flows show strong deviations from Kolomogorov-like behaviour resulting from more complex dynamics compared to Newtonian turbulence. We now study the nature of mixing in polymeric turbulence via Eulerian passive scalar fields of varying molecular diffusivities, given by the Schmidt number Sc. We show that polymeric turbulence is a less efficient mixer than the Newtonian one at small to moderate Sc numbers. Newtonian scalar turbulence (NST) forms large islands of fluctuations with extended, contiguous fronts. In contrast, polymeric scalar turbulence (PST) is marked by small, interspersed patches of strong but less intermittent fluctuations. These patches collectively comprise a larger volume fraction of strong fluctuations, indicating a less efficient mixing, alongwith smaller scalar gradients and therefore smaller average flux across their boundaries. Box counting dimensions reveal a smoother and more space filling nature of patch boundaries in PST compared to NST fronts. Finally, spatial changes of the scalar are stronger in PST, but with a slower self-similar growth and less intermittency as revealed by the kurtosis of scalar differences. Overall, these observations hint at reduced mixing in PST where fluctuations are typically stronger while the average scalar flux is smaller in a stationary state.

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

3 major / 2 minor

Summary. The manuscript examines passive scalar mixing in polymeric turbulence (PST) versus Newtonian turbulence (NST) via direct numerical simulations of Eulerian scalar fields at varying Schmidt numbers Sc. It reports that PST produces small, interspersed patches of strong fluctuations (larger volume fraction but weaker gradients and boundary fluxes) rather than the large contiguous islands seen in NST, concluding that polymeric flows are less efficient mixers at small-to-moderate Sc; supporting diagnostics include box-counting dimensions of patch boundaries and kurtosis of scalar increments.

Significance. If the reported morphological and statistical distinctions prove robust, the work would establish a concrete mechanism by which polymer elasticity suppresses scalar mixing efficiency, with potential relevance to drag-reduced flows and industrial mixing processes. The absence of any resolution study, parameter sweep, or cross-model validation, however, leaves the central claim on a weak empirical footing.

major comments (3)
  1. [Abstract] Abstract and simulation description: no grid resolution, time-stepping scheme for the conformation tensor, or statistical convergence diagnostics (error bars, run times, or ensemble sizes) are provided for either the NST or PST cases. Without these, it is impossible to rule out that the reported patch/island morphology and reduced flux arise from under-resolved viscoelastic stresses or differing numerical diffusion between the two runs.
  2. [Abstract] The claim that PST exhibits 'smaller scalar gradients and therefore smaller average flux across their boundaries' is central to the reduced-mixing conclusion, yet the manuscript supplies no explicit computation of the scalar flux (e.g., via the surface integral of |∇θ| or the volume-averaged |u·∇θ|), nor any normalization that guarantees identical mean dissipation rates between NST and PST. This leaves open the possibility that the observed difference is an artifact of unmatched velocity statistics rather than polymer physics.
  3. [Abstract] The FENE-P (or equivalent) constitutive model parameters (Weissenberg number, polymer viscosity ratio) are not varied; the entire PST versus NST comparison rests on a single parameter point. A load-bearing robustness test would require at least a modest sweep in Wi or β to confirm that the patchiness and flux reduction persist outside the chosen regime.
minor comments (2)
  1. [Abstract] The abstract refers to 'Box counting dimensions' and 'kurtosis of scalar differences' without defining the precise quantities or the scale range over which they are computed.
  2. Notation for the scalar field (θ or c) and the precise definition of 'strong fluctuations' (e.g., a threshold on |θ'|) should be stated explicitly in the main text.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each major comment point by point below, indicating the revisions we will make.

read point-by-point responses

  1. Referee: [Abstract] Abstract and simulation description: no grid resolution, time-stepping scheme for the conformation tensor, or statistical convergence diagnostics (error bars, run times, or ensemble sizes) are provided for either the NST or PST cases. Without these, it is impossible to rule out that the reported patch/island morphology and reduced flux arise from under-resolved viscoelastic stresses or differing numerical diffusion between the two runs.

    Authors: We acknowledge that these numerical details are not summarized in the abstract. We will revise the manuscript to include explicit statements on grid resolution, the time-stepping scheme employed for the conformation tensor, run times, ensemble sizes, and statistical convergence diagnostics (including error bars) in both the abstract and a dedicated methods subsection. We will also add a short resolution study confirming that the reported morphological and flux differences persist under increased resolution and are not attributable to under-resolved stresses or mismatched numerical diffusion. revision: yes

  2. Referee: [Abstract] The claim that PST exhibits 'smaller scalar gradients and therefore smaller average flux across their boundaries' is central to the reduced-mixing conclusion, yet the manuscript supplies no explicit computation of the scalar flux (e.g., via the surface integral of |∇θ| or the volume-averaged |u·∇θ|), nor any normalization that guarantees identical mean dissipation rates between NST and PST. This leaves open the possibility that the observed difference is an artifact of unmatched velocity statistics rather than polymer physics.

    Authors: We agree that direct computation of the scalar flux would strengthen the central claim. The current manuscript infers the reduced flux from the observed smaller gradients and patch statistics, but we will add explicit calculations of the surface integral of |∇θ| and the volume-averaged |u·∇θ| in the revised version. We will also report the normalization procedure used to match mean dissipation rates between the NST and PST cases, thereby confirming that the flux reduction is due to polymer elasticity rather than unmatched velocity statistics. revision: yes

  3. Referee: [Abstract] The FENE-P (or equivalent) constitutive model parameters (Weissenberg number, polymer viscosity ratio) are not varied; the entire PST versus NST comparison rests on a single parameter point. A load-bearing robustness test would require at least a modest sweep in Wi or β to confirm that the patchiness and flux reduction persist outside the chosen regime.

    Authors: The comparison is performed at a single, representative FENE-P parameter set corresponding to moderate elasticity. While a full parameter sweep would be desirable for robustness, the associated computational expense for additional DNS runs is high. In the revised manuscript we will expand the discussion to justify the chosen parameters, explain why the qualitative features of patchiness and reduced mixing are expected to hold for Wi > 1 on the basis of existing literature, and explicitly note the limitation of the single-point study. We view this as a partial revision that addresses the concern without new simulations. revision: partial

Circularity Check

0 steps flagged

No circularity: results are direct DNS outputs with no derivation chain

full rationale

The paper presents empirical observations from direct numerical simulations of scalar fields in Newtonian versus polymeric turbulence. Claims about patchiness, volume fractions of strong fluctuations, scalar gradients, average flux, box-counting dimensions, and kurtosis of differences are computed statistics extracted from the Eulerian fields, not predictions or derivations that reduce to the simulation inputs by construction. No equations, fitted parameters renamed as predictions, self-citations invoked as uniqueness theorems, or ansatzes smuggled via prior work appear in the load-bearing steps. The central distinction between NST islands and PST patches follows directly from the simulated data under the stated constitutive model and parameters, making the analysis self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

Abstract-only review provides no explicit equations or parameter lists; inferred standard assumptions for polymeric turbulence simulations are noted below.

free parameters (2)
  • Schmidt number Sc
    Varied parameter controlling molecular diffusivity of the scalar; central to the comparison at small to moderate values.
  • Polymer relaxation time or Weissenberg number
    Key control parameter for polymeric effects, implied but not quantified in abstract.
axioms (2)
  • domain assumption Polymeric fluid obeys a viscoelastic constitutive model (e.g., Oldroyd-B or FENE-P) coupled to incompressible Navier-Stokes equations
    Standard modeling choice for direct numerical simulations of polymeric turbulence.
  • standard math The scalar field is passive and advected by the velocity field with molecular diffusion
    Standard passive scalar transport equation in turbulence studies.

pith-pipeline@v0.9.0 · 5497 in / 1351 out tokens · 58881 ms · 2026-05-10T16:19:58.534819+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

27 extracted references · 27 canonical work pages

  1. [1]

    M.andMungal M

    [1]White C. M.andMungal M. G.,Annu. Rev. Fluid Mech.,40(2008)

    work page 2008

  2. [2]

    [2]Benzi R.andChing E. S.,Annu. Rev. Condens. Matter Phys.,9(2018)

    work page 2018

  3. [3]

    E.andHof B.,Annual Re- view of Fluid Mechanics,55(2023)

    [3]Dubief Y., Terrapon V. E.andHof B.,Annual Re- view of Fluid Mechanics,55(2023)

    work page 2023

  4. [4]

    p-6 Inefficient Mixing in Polymeric Turbulence at small Sc [4]Toms B. A.,Proc. 1st Int. Cong. Rheology, Amsterdam, 1948,135(1948) . [5]Lumley J. L.,Journal of Polymer Science: Macromolec- ular Reviews,7(1973)

    work page 1948

  5. [5]

    [6]Perlekar P., Mitra D.andPandit R.,Phys. Rev. Lett.,97(2006) 264501. [7]Bos W. J., Shao X., Wu T.andFang L.,Journal of Fluid Mechanics,1019(2025) A32. [8]Serafini F., Battista F., Gualtieri P.andCasciola C.,Journal of Fluid Mechanics,1007(2025) A10. [9]Samanta D., Dubief Y., Holzner M., Sch ¨afer C., Morozov A. N., Wagner C.andHof B.,Proceedings of the ...

    work page 2006

  6. [6]

    [13]Berti S., Bistagnino A., Boffetta G., Celani A.and Musacchio S.,Phys. Rev. E,77(2008) 055306. [14]Steinberg V.,Annual Review of Fluid Mechanics,53 (2021)

    work page 2008

  7. [7]

    E.,International Journal of Multiphase Flow, (2023) 104546

    [15]Soligo G.andRosti M. E.,International Journal of Multiphase Flow, (2023) 104546. [16]Singh R. K., Perlekar P., Mitra D.andRosti M. E., Nat. Comm.,15(2024)

    work page 2023

  8. [8]

    R., Picardo J

    [17]Yerasi S. R., Picardo J. R., Gupta A.andVincenzi D.,Journal of Fluid Mechanics,1000(2024) A37. [18]Foggi Rota G., Amor C., Le Clainche S.andRosti M. E.,Phys. Rev. Fluids,9(2024) L122602. [19]Foggi Rota G., Singh R. K., Chiarini A., Amor C., Soligo G., Mitra D.andRosti M. E.,International Journal of Multiphase Flow,197(2026) 105630. [20]Zhang Y.-B., Bo...

    work page 2024

  9. [9]

    E.,Annual Review of Fluid Mechanics,37 (2005)

    [29]Dimotakis P. E.,Annual Review of Fluid Mechanics,37 (2005)

    work page 2005

  10. [10]

    I.andSreenivasan K

    [30]Abarzhi S. I.andSreenivasan K. R.,Philosophical Transactions of the Royal Society A: Mathematical, Phys- ical and Engineering Sciences,368(2010)

    work page 2010

  11. [11]

    [31]Gregg M., D’Asaro E., Riley J.andKunze E.,An- nual Review of Marine Science,10(2018)

    work page 2018

  12. [12]

    [32]Villermaux E.,Annual Review of Fluid Mechanics,51 (2019)

    work page 2019

  13. [13]

    [33]Caulfield C.,Annual Review of Fluid Mechanics,53 (2021)

    work page 2021

  14. [14]

    [34]Warhaft Z.,Annual Review of Fluid Mechanics,32 (2000)

    work page 2000

  15. [15]

    [35]Obukhov A. M.,Izv. Akad. Nauk SSSR, Ser. Geogr. Ge- ofiz.,13(1949)

    work page 1949

  16. [16]

    [36]Corrsin S.,Journal of Applied Physics,22(1951)

    work page 1951

  17. [17]

    [37]Falkovich G., Gawe ¸dzki K.andVergassola M.,Rev. Mod. Phys.,73(2001)

    work page 2001

  18. [18]

    [38]Chertkov M.andLebedev V.,Phys. Rev. Lett.,90 (2003) 034501. [39]Sreenivasan K. R.andSchumacher J.,Philosophical Transactions of the Royal Society A: Mathematical, Phys- ical and Engineering Sciences,368(2010)

    work page 2003

  19. [19]

    P., Schumacher J., Sreenivasan K

    [40]Iyer K. P., Schumacher J., Sreenivasan K. R.and Yeung P. K.,Phys. Rev. Lett.,121(2018) 264501. [41]Sreenivasan K. R.,Proceedings of the National Academy of Sciences,116(2019) 18175. [42]Buaria D., Clay M. P., Sreenivasan K. R.andYe- ung P. K.,Phys. Rev. Lett.,126(2021) 074501. [43]Buaria D., Clay M. P., Sreenivasan K. R.andYe- ung P. K.,Phys. Rev. Let...

    work page 2018

  20. [20]

    [45]Burghelea T., Segre E., Bar-Joseph I., Groisman A.andSteinberg V.,Phys. Rev. E,69(2004) 066305. [46]Thomases B., Shelley M.andThiffeault J.-L., Physica D: Nonlinear Phenomena,240(2011) 1602 spe- cial Issue: Fluid Dynamics: From Theory to Experiment. [47]Li D.-Y., Li X.-B., Zhang H.-N., Li F.-C., Qian S.and Joo S. W.,Microfluidics and Nanofluidics,21(2017)

    work page 2004

  21. [21]

    C., Pinho F

    [48]Guimar ˜aes M. C., Pinho F. T.andda Silva C. B., Phys. Rev. Fluids,10(2025) 023303. [49]Dzanic V., From C. S.andSauret E.,PNAS Nexus, 4(2025) pgaf231. [50]Schumacher J.andSreenivasan K. R.,Physics of Flu- ids,17(2005) 125107. [51]Pope S.,International Journal of Engineering Science, 26(1988)

    work page 2025

  22. [22]

    R., Emerson B., Lieuwen T

    [52]Roy A., Picardo J. R., Emerson B., Lieuwen T. C. andSujith R.,Journal of Fluid Mechanics,957(2023) A21. [53]Constantin P., Procaccia I.andSreenivasan K. R., Phys. Rev. Lett.,67(1991)

    work page 2023

  23. [23]

    [54]Grossmann S.andLohse D.,Europhysics Letters,27 (1994)

    work page 1994

  24. [24]

    P., Schumacher J., Sreenivasan K

    [55]Iyer K. P., Schumacher J., Sreenivasan K. R.and Yeung P. K.,Phys. Rev. Fluids,5(2020) 044501. [56]Kolmogorov A. N.,Doklady Akademii Nauk SSSR,30 (1941)

    work page 2020

  25. [25]

    R., Hunt J

    [57]Sreenivasan K. R., Hunt J. C. R., Phillips O. M. andWilliams D.,Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences, 434(1991)

    work page 1991

  26. [26]

    D.,Physics of Fluids,6(1994)

    [58]Holzer M.andSiggia E. D.,Physics of Fluids,6(1994)

    work page 1994

  27. [27]

    [59]Celani A., Lanotte A., Mazzino A.andVergassola M.,Physics of Fluids,13(2001)

    work page 2001