pith. sign in

arxiv: 2605.12883 · v1 · pith:C6JH6NCBnew · submitted 2026-05-13 · 🧮 math.AP

Mixing and Small-Scale Formation in a Passive Divergence-Free Vector Field

Anuj Kumar , Franziska Weber This is my paper

Pith reviewed 2026-05-14 19:01 UTC · model grok-4.3

classification 🧮 math.AP
keywords mixingdivergence-free vector fieldspassive transportSobolev normsH^{-α} decaynumerical optimizationsmall-scale formationincompressible flows
0
0 comments X

The pith

Numerical simulations indicate that passive divergence-free vector fields mix at least exponentially when the advecting field is chosen at each instant to maximize instantaneous decay of the negative Sobolev norm.

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

The paper studies the mixing of a divergence-free passive vector field u transported by a divergence-free field U on the torus, with U constrained to have bounded W^{1,q} norm. Mixing is measured by the decay of the homogeneous H^{-α} norm of u, and the pressure term in the evolution equation prevents a simple Lagrangian approach. The authors first give conditions ensuring existence and uniqueness of solutions, then obtain lower bounds on mixing rates for ranges of q and α. Numerical simulations that select U at each time to maximize the instantaneous decay rate supply evidence that the optimal mixing is at least exponential in time. This construction supplies a tractable model for examining how small scales form in divergence-free fields and for testing simplified versions of questions that arise in the incompressible Euler and Navier-Stokes equations.

Core claim

When the advecting field U is chosen at every instant to be divergence-free, to obey the W^{1,q} bound, and to maximize the instantaneous rate of decay of the H^{-α} norm of the passive field u, numerical evolution on the torus produces at least exponential decay of that norm.

What carries the argument

The time-dependent instantaneous maximizer of the decay rate of ||u||_{H^{-α}} subject to the divergence-free constraint and the W^{1,q} bound on U, which is used to generate the mixing flow.

If this is right

  • Existence and uniqueness of solutions hold when U satisfies the stated regularity conditions.
  • Lower bounds on the mixing rate are derived for different ranges of the parameters q and α.
  • The construction supplies a simplified setting in which to study small-scale formation mechanisms inside divergence-free vector fields.
  • The same framework can be used to formulate reduced versions of open questions about the incompressible Euler and Navier-Stokes equations.

Where Pith is reading between the lines

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

  • If the exponential rate persists, the model would imply that rapid small-scale generation occurs even without diffusion or forcing.
  • The same instantaneous-optimization idea could be tested on the diffusive version of the equation mentioned in the paper.
  • The approach may connect to questions of enhanced dissipation or optimal control in incompressible transport.
  • It would be interesting to see whether the exponential rate survives when the W^{1,q} bound is replaced by a weaker constraint.

Load-bearing premise

That a divergence-free U with bounded W^{1,q} norm exists at each instant which instantaneously maximizes the decay of the H^{-α} norm of u, and that the resulting evolution remains well-defined.

What would settle it

A longer numerical run of the same optimization procedure in which the H^{-α} norm of u stops decaying exponentially or the solution ceases to exist while the W^{1,q} bound is respected.

Figures

Figures reproduced from arXiv: 2605.12883 by Anuj Kumar, Franziska Weber.

Figure 1
Figure 1. Figure 1: Snapshots of U (column 1), u (column 2) and p (column 3) at various times in the case α = 1. In columns 1 and 2, the colorbars represent the magnitudes of the velocity fields U and u, respectively, while in column 3 the colorbar represents the pressure p. The plots of U and u are overlaid with streamlines and arrows indicating the direction of the flow. 19 [PITH_FULL_IMAGE:figures/full_fig_p019_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Same as Figure [PITH_FULL_IMAGE:figures/full_fig_p020_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: b). This is consistent with the fact that H−1 places comparatively less emphasis on smaller scales than H−1/2 . (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p021_3.png] view at source ↗
read the original abstract

We study mixing for a divergence-free passive vector field $u$ transported by another divergence-free vector field $U$, where $u$ evolves according to $ \partial_t u + (U \cdot \nabla) u + \nabla p = 0.$ In recent years, a lot of attention has been given to the question of optimal mixing in the scalar case, where there is a Sobolev constraint on the advecting velocity. In the vector setting considered here, however, the pressure term introduces substantial difficulties, since the simple Lagrangian perspective available in the scalar case is no longer applicable. In this paper, we investigate mixing on a torus $\mathbb{T}^d$ under the assumption that the field $U$ satisfies $ \|U(t)\|_{W^{1,q}} \leq C $ and we quantify mixing through the decay of the homogeneous $ H^{-\alpha}$ norm of $u$. We start with establishing conditions on $U$ that guarantee existence and uniqueness of solutions. We then derive lower bounds on the mixing rate for various ranges of $q$ and $\alpha$. In addition, we carry out numerical simulations of mixing by choosing, at each time instant, a field $U$ that maximizes the instantaneous decay of the $ H^{-\alpha}$ norm. These simulations provide evidence that the optimal mixing rate is at least exponential in time. More broadly, we view the present model and its diffusive analogue, as a useful framework for probing mechanisms of small-scale formation in divergence-free vector fields and for formulating simplified versions of open questions related to the incompressible Euler and Navier--Stokes equations.

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

1 major / 3 minor

Summary. The paper studies mixing for a divergence-free passive vector field u transported by another divergence-free vector field U with pressure term, evolving as ∂_t u + (U · ∇) u + ∇p = 0 on the torus T^d. Under the assumption ||U(t)||_{W^{1,q}} ≤ C, it establishes conditions guaranteeing existence and uniqueness of solutions, derives lower bounds on the mixing rate via decay of the homogeneous H^{-α} norm of u for various ranges of q and α, and conducts numerical simulations in which U is chosen at each instant to maximize the instantaneous decay rate of ||u||_{H^{-α}}. These simulations are presented as evidence that the optimal mixing rate is at least exponential in time. The work frames the model as a framework for probing small-scale formation in divergence-free vector fields and simplified versions of questions in incompressible Euler and Navier-Stokes equations.

Significance. The lower bounds on mixing rates under Sobolev constraints extend scalar mixing results to the vector setting with pressure, a non-trivial step given the loss of the Lagrangian perspective. If the numerical evidence for exponential decay can be rigorously shown to remain within the paper's analytic regime, it would offer valuable insight into optimal rates and mechanisms of small-scale formation. The positioning as a simplified model for open problems in fluid equations adds broader relevance, though the strength of the exponential claim depends on validating the closed-loop numerics.

major comments (1)
  1. [Numerical Simulations section] Numerical Simulations section: The claim that simulations provide evidence for at least exponential optimal mixing relies on selecting, at each time, a divergence-free U maximizing the instantaneous decay of ||u||_{H^{-α}} subject to ||U||_{W^{1,q}} ≤ C. This closed-loop choice makes U a (possibly nonlocal) functional of the current u. It is not shown that the resulting time-dependent U automatically satisfies the regularity hypotheses used to prove existence and uniqueness for the transport equation with pressure, nor that the W^{1,q} bound remains uniformly controlled. Without such justification or a posteriori verification that the discrete maximizer stays admissible, the observed exponential decay may occur outside the regime where the analytic lower bounds apply, undermining the central numerical conclusion.
minor comments (3)
  1. The description of the numerical scheme lacks details on discretization (e.g., finite elements or spectral methods), time-stepping, spatial resolution, and the concrete optimization procedure used to compute the maximizing U at each step; these are needed to assess reliability of the exponential decay observation.
  2. [Existence and Uniqueness section] Clarify the precise ranges of q and α for which existence/uniqueness hold and how they align with the parameter regimes used for the lower bounds and the numerical experiments.
  3. Figure captions and the numerical results section should explicitly report the observed decay rates (e.g., fitted exponents) and any checks performed to confirm that ||U(t)||_{W^{1,q}} remained bounded by C throughout the runs.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive comments. We address the major concern regarding the numerical simulations below.

read point-by-point responses

  1. Referee: [Numerical Simulations section] Numerical Simulations section: The claim that simulations provide evidence for at least exponential optimal mixing relies on selecting, at each time, a divergence-free U maximizing the instantaneous decay of ||u||_{H^{-α}} subject to ||U||_{W^{1,q}} ≤ C. This closed-loop choice makes U a (possibly nonlocal) functional of the current u. It is not shown that the resulting time-dependent U automatically satisfies the regularity hypotheses used to prove existence and uniqueness for the transport equation with pressure, nor that the W^{1,q} bound remains uniformly controlled. Without such justification or a posteriori verification that the discrete maximizer stays admissible, the observed exponential decay may occur outside the regime where the analytic lower bounds apply, undermining the central numerical conclusion.

    Authors: We agree that a more explicit justification is needed to connect the numerics to the analytic regime. By construction, the optimization at each step is performed over divergence-free fields satisfying the W^{1,q} bound, so the constraint is enforced directly. The existence/uniqueness hypotheses in the paper require precisely that U be divergence-free and bounded in W^{1,q}; these are preserved in the discrete finite-dimensional maximization space we employ. In the revised version we will add a detailed description of the discretization and optimization procedure, together with a posteriori verification (where computationally feasible) that the computed maximizers remain admissible for the entire simulation interval. This will strengthen the link between the observed decay and the analytic lower bounds. revision: yes

Circularity Check

0 steps flagged

No significant circularity; analytic bounds and numerical evidence remain independent

full rationale

The paper first proves existence/uniqueness and derives explicit lower bounds on the decay of ||u||_{H^{-α}} under the standing assumption ||U(t)||_{W^{1,q}} ≤ C (quoted in abstract and section on well-posedness). These bounds are obtained directly from the transport equation and Sobolev embedding estimates without fitting or redefinition. The numerical component separately constructs, at each instant, a divergence-free U that maximizes the instantaneous decay rate of the H^{-α} norm while respecting the same W^{1,q} bound; the observed exponential decay is an output of the time-stepping simulation, not a parameter fitted to the target rate and then relabeled as a prediction. No self-citation chain is invoked to justify the central claims, and the derivation chain does not collapse any quantity to its own input by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard existence theory for linear transport equations on the torus together with the imposed uniform bound on the first derivatives of U; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption Both u and U are divergence-free
    Required to close the system with the pressure term and to preserve the incompressibility constraint.
  • domain assumption U belongs to a bounded set in W^{1,q} for some q
    This regularity assumption is used to guarantee existence/uniqueness and to control the mixing rate.

pith-pipeline@v0.9.0 · 5588 in / 1468 out tokens · 50253 ms · 2026-05-14T19:01:50.766280+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

34 extracted references · 34 canonical work pages

  1. [1]

    Alberti, G

    G. Alberti, G. Crippa, and A. L. Mazzucato. Exponential self-similar mixing by incompressible flows.J. Amer. Math. Soc., 32(2):445–490, 2019

    work page 2019

  2. [2]

    A. Bressan. A lemma and a conjecture on the cost of rearrangements.Rendiconti del Seminario Matematico della Universita di Padova, 110:97–102, 2003

    work page 2003

  3. [3]

    Brezis and P

    H. Brezis and P. Mironescu. Gagliardo-Nirenberg inequalities and non-inequalities: the full story. Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire, 35(5):1355–1376, 2018. 24 Mixing and Small-Scale Formation in a Passive VectorA. Kumar and F. Weber

    work page 2018

  4. [4]

    Brezis and P

    H. Brezis and P. Mironescu. Where Sobolev interacts with Gagliardo-Nirenberg.J. Funct. Anal., 277(8):2839–2864, 2019

    work page 2019

  5. [5]

    Bru´ e, M

    E. Bru` e, M. Colombo, and A. Kumar. Sharp Nonuniqueness in the Transport Equation with Sobolev Velocity Field.arXiv preprint arXiv:2405.01670, 2024

    work page arXiv 2024

  6. [6]

    D. R. Caldwell and J. N. Mourn. Turbulence and mixing in the ocean.Rev. Geophys., 33(S2):1385–1394, 1995

    work page 1995

  7. [7]

    Crippa and C

    G. Crippa and C. De Lellis. Regularity and compactness for the diperna–lions flow.Hyperbolic problems: theory, numerics, applications, pages 423–430, 2008

    work page 2008

  8. [8]

    N. Depauw. Non unicit´ e des solutions born´ ees pour un champ de vecteurs BV en dehors d’un hyperplan.Comptes rendus. Math´ ematique, 337(4):249–252, 2003

    work page 2003

  9. [9]

    Di Nezza, G

    E. Di Nezza, G. Palatucci, and E. Valdinoci. Hitchhiker’s guide to the fractional Sobolev spaces. Bull. Sci. Math., 136(5):521–573, 2012

    work page 2012

  10. [10]

    C. R. Doering and I. Tobasco. On the optimal design of wall-to-wall heat transport.Comm. Pure Appl. Math., 72(11):2385–2448, 2019

    work page 2019

  11. [11]

    T. D. Drivas. Self-regularization in turbulence from the Kolmogorov 4/5-law and alignment. Philosophical Transactions of the Royal Society A, 380(2226):20210033, 2022

    work page 2022

  12. [12]

    T. D. Drivas, T. M. Elgindi, G. Iyer, and I.-J. Jeong. Anomalous dissipation in passive scalar transport.Arch. Ration. Mech. Anal., 243(3):1151–1180, 2022

    work page 2022

  13. [13]

    T. M. Elgindi and K. Liss. Norm growth, non-uniqueness, and anomalous dissipation in passive scalars.Archive for Rational Mechanics and Analysis, 248(6):120, 2024

    work page 2024

  14. [14]

    T. M. Elgindi and A. Zlatoˇ s. Universal mixers in all dimensions.Adv. Math., 356:106807, 33, 2019

    work page 2019

  15. [15]

    Escauriaza, G

    L. Escauriaza, G. Seregin, and V. Sverak. On l 3, infinity-solutions to the navier-stokes equations and backward uniqueness. 2002

    work page 2002

  16. [16]

    Frisch.Turbulence: the legacy of AN Kolmogorov

    U. Frisch.Turbulence: the legacy of AN Kolmogorov. Cambridge university press, 1995

    work page 1995

  17. [17]

    G. Iyer, A. Kiselev, and X. Xu. Lower bounds on the mix norm of passive scalars advected by incompressible enstrophy-constrained flows.Nonlinearity, 27(5):973, 2014

    work page 2014

  18. [18]

    Kaneda, T

    Y. Kaneda, T. Ishihara, M. Yokokawa, K. Itakura, and A. Uno. Energy dissipation rate and energy spectrum in high resolution direct numerical simulations of turbulence in a periodic box. Physics of Fluids, 15(2):L21–L24, 2003. 25 Mixing and Small-Scale Formation in a Passive VectorA. Kumar and F. Weber

    work page 2003

  19. [19]

    A. Kumar. Almost anomalous dissipation in advection-diffusion of a divergence-free passive vector.arXiv preprint arXiv:2409.15000, 2024

    work page arXiv 2024

  20. [20]

    A. Kumar. Nonuniqueness of trajectories on a set of full measure for sobolev vector fields.Archive for Rational Mechanics and Analysis, 248(6):114, 2024

    work page 2024

  21. [21]

    A. Kumar. Three dimensional branching pipe flows for optimal scalar transport between walls. Nonlinearity, 37(11):115011, 2024

    work page 2024

  22. [22]

    Lin, J.-L

    Z. Lin, J.-L. Thiffeault, and C. R. Doering. Optimal stirring strategies for passive scalar mixing. J. Fluid Mech., 675:465–476, 2011

    work page 2011

  23. [23]

    Munk and C

    W. Munk and C. Wunsch. Abyssal recipes II: Energetics of tidal and wind mixing.Deep Sea Research Part I: Oceanographic Research Papers, 45(12):1977–2010, 1998

    work page 1977

  24. [24]

    J. M. Ottino. Mixing and chemical reactions a tutorial.Chemical Engineering Science, 49(24):4005–4027, 1994

    work page 1994

  25. [25]

    E. L. Paul, V. A. Atiemo-Obeng, S. M. Kresta, et al.Handbook of industrial mixing. Wiley Online Library, 2004

    work page 2004

  26. [26]

    B. R. Pearson, P.- ˚A. Krogstad, and W. van de Water. Measurements of the turbulent energy dissipation rate.Physics of fluids, 14(3):1288–1290, 2002

    work page 2002

  27. [27]

    C. Seis. Maximal mixing by incompressible fluid flows.Nonlinearity, 26(12):3279, 2013

    work page 2013

  28. [28]

    J. Simon. Compact sets in the spaceL p(0, T;B).Ann. Mat. Pura Appl. (4), 146:65–96, 1987

    work page 1987

  29. [29]

    K. R. Sreenivasan. On the scaling of the turbulence energy dissipation rate.Phys. Fluids, 27(5):1048–1051, 1984

    work page 1984

  30. [30]

    K. R. Sreenivasan. An update on the energy dissipation rate in isotropic turbulence.Physics of Fluids, 10(2):528–529, 1998

    work page 1998

  31. [31]

    Thiffeault

    J.-L. Thiffeault. Using multiscale norms to quantify mixing and transport.Nonlinearity, 25(2):R1–R44, 2012

    work page 2012

  32. [32]

    I. Tobasco. Optimal cooling of an internally heated disc.Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 380(2225), 2022

    work page 2022

  33. [33]

    Wunsch and R

    C. Wunsch and R. Ferrari. Vertical mixing, energy, and the general circulation of the oceans. Annu. Rev. Fluid Mech., 36:281–314, 2004

    work page 2004

  34. [34]

    Yao and A

    Y. Yao and A. Zlatoˇ s. Mixing and un-mixing by incompressible flows.J. Eur. Math. Soc. (JEMS), 19(7):1911–1948, 2017. 26

    work page 1911