pith. sign in

arxiv: 2604.09507 · v1 · submitted 2026-04-10 · ⚛️ physics.flu-dyn

Early emergence of ultimate-like transport in two-dimensional turbulent thermomagnetic convection

Paolo Capobianchi This is my paper

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

classification ⚛️ physics.flu-dyn
keywords thermomagnetic convectionultimate scalingturbulent transportthermal plumesmagnetic forcedirect numerical simulationhigh Prandtl numbersquare cavity
0
0 comments X

The pith

Magnetic force enables early ultimate-like scalings in two-dimensional thermomagnetic convection.

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

Direct numerical simulations of a high-Prandtl fluid in a square cavity under thermomagnetic convection show that after the flow turns turbulent, both the Nusselt number for heat transport and the Reynolds number for flow speed scale as the square root of the magnetic Rayleigh number. This regime begins right after the laminar-to-turbulent transition and continues across more than an order of magnitude in the driving parameter. The authors link the behavior to the magnetic force detaching thermal plumes from the walls and carrying them through the fluid interior. A reader would care because ordinary thermal convection reaches such efficient transport only at much higher driving strengths, while here an additional force makes it possible sooner.

Core claim

In two-dimensional direct numerical simulations of thermomagnetic convection, the magnetic force promotes the ejection of thermal plumes from the boundary layers and their advection across the fluid bulk, producing an ultimate-like regime in which Nu ∼ Ra_m^{1/2} and Re ∼ Ra_m^{1/2} that appears immediately after the laminar-to-turbulent transition and persists for more than an order of magnitude.

What carries the argument

The magnetic body force that detaches thermal plumes from the walls and advects them through the fluid bulk.

If this is right

  • Heat and momentum transport follow ultimate scalings once turbulence is reached.
  • The magnetic force directly controls plume ejection and bulk crossing, sustaining the scalings over a wide parameter range.
  • The regime appears at lower driving strengths than in standard Rayleigh-Bénard convection without magnetic effects.
  • Theoretical arguments based on plume dynamics confirm the observed power-law exponents.

Where Pith is reading between the lines

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

  • Three-dimensional simulations could test whether additional flow instabilities disrupt the plume-advection mechanism.
  • The same magnetic assistance might allow efficient transport in other buoyancy-driven systems at moderate driving levels.
  • Laboratory experiments in thin-layer or quasi-two-dimensional setups could verify the scalings before full three-dimensional checks.

Load-bearing premise

The scaling relations found in two-dimensional simulations reflect the true physical mechanism and will continue to hold in three-dimensional geometries and laboratory experiments.

What would settle it

A three-dimensional simulation or controlled experiment at high magnetic Rayleigh numbers that measures a clear deviation from Nu scaling as the square root of Ra_m would falsify the claim that the regime emerges and persists as described.

Figures

Figures reproduced from arXiv: 2604.09507 by Paolo Capobianchi.

Figure 1
Figure 1. Figure 1: Comparison of viscous (δu) and thermal (δθ) bound￾ary layers along z at the mid-plane (“Mid”, solid) and cor￾ner (“Cor”, dotted) of the large-scale circulation for Ram = 2.08 × 107 . Thicknesses are defined as δu = 1/|∂z(u/umax)|w and δθ = (θw −θb)/|∂zθ|w, where subscript w denotes wall val￾ues. Inset table: average thicknesses for all Ram, calculated across all computational cells of both heated boundarie… view at source ↗
Figure 2
Figure 2. Figure 2: (a) Normalized time-averaged velocity vector field [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) Nu and (b) Re as functions of the magnetic Rayleigh number Ram. Dashed and dash-dotted lines indicate the theoretical Ra1/2 m scaling, with prefactors obtained via best-fit. The inset in (a) shows the compensated plot Nu/Ra0.4 m versus Ram. such as the presence of distinct corner rolls (see Sup￾plemental Material [33]), the structure of the LSCs does not always show evident columnar structures as in th… view at source ↗
read the original abstract

Scaling laws for turbulent thermomagnetic convection of a high-Pr fluid in a square cavity are obtained through direct numerical simulations and formulated via theoretical arguments informed by the numerical data. A regime consistent with an ultimate-like scaling $Nu \sim Ra_m^{1/2} and Re \sim Ra_m^{1/2} emerges after the laminar-to-turbulent transition and persists for more than an order of magnitude. Evidence is provided that this heretofore unseen behavior stems from the ability of the magnetic force to facilitate the ejection and advection of thermal plumes across the fluid bulk.

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 reports results from direct numerical simulations of two-dimensional thermomagnetic convection of a high-Prandtl-number fluid in a square cavity. It identifies a post-transition regime in which both the Nusselt number and Reynolds number follow an ultimate-like scaling Nu ∼ Ra_m^{1/2} and Re ∼ Ra_m^{1/2} that persists over more than an order of magnitude in the magnetic Rayleigh number. The authors attribute this behavior to the magnetic body force enabling efficient ejection and advection of thermal plumes across the bulk and formulate theoretical scaling arguments informed by the simulation data.

Significance. If the reported scalings and mechanism are robust, the work would demonstrate an unusually early onset of ultimate transport in thermomagnetic convection, providing a concrete example of how an external body force can alter the transition to high-transport regimes. The data-driven formulation of theoretical arguments is a methodological strength that could be extended to other MHD or buoyancy-driven flows.

major comments (2)
  1. [Numerical methods and results sections] The central claim that an ultimate-like regime with exponent 1/2 emerges and persists for more than an order of magnitude rests on 2D DNS, yet the manuscript provides no information on grid resolution, time-averaging windows, or convergence checks with respect to spatial or temporal discretization. Without these, it is impossible to rule out that the observed exponent is influenced by under-resolved small-scale dissipation or insufficient statistical sampling.
  2. [Discussion of the mechanism] In two-dimensional turbulence the inverse energy cascade produces large-scale coherent structures that can dominate bulk transport even in the absence of the proposed magnetic-plume mechanism. The paper does not present comparisons with purely thermal (non-magnetic) 2D runs at the same Ra or with 3D simulations to isolate the magnetic force contribution from 2D condensation effects; such tests are load-bearing for the attribution of the scaling to magnetic plume advection.
minor comments (2)
  1. [Abstract and introduction] The abstract and introduction should explicitly state the Prandtl number(s) simulated and the precise definition of the magnetic Rayleigh number Ra_m upon first use.
  2. [Figures] Figure captions and axis labels should indicate the precise range of Ra_m over which the 1/2 scaling is claimed and whether the exponent is obtained from a fit or from visual inspection.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive and insightful comments, which have helped us improve the clarity and robustness of the manuscript. We address each major comment below and have revised the manuscript to incorporate additional details and clarifications where possible.

read point-by-point responses

  1. Referee: [Numerical methods and results sections] The central claim that an ultimate-like regime with exponent 1/2 emerges and persists for more than an order of magnitude rests on 2D DNS, yet the manuscript provides no information on grid resolution, time-averaging windows, or convergence checks with respect to spatial or temporal discretization. Without these, it is impossible to rule out that the observed exponent is influenced by under-resolved small-scale dissipation or insufficient statistical sampling.

    Authors: We agree that explicit documentation of the numerical parameters is essential for validating the reported scalings. In the revised manuscript, we have expanded the Numerical Methods section to include a table summarizing the grid resolutions (in terms of grid points per thermal boundary layer and Kolmogorov scale) for each Ra_m, the duration of the time-averaging windows used for Nu and Re statistics (typically exceeding 100 free-fall times after discarding transients), and the results of convergence tests performed by doubling the grid resolution and halving the time step at selected Ra_m values. These tests confirm that the observed Nu ~ Ra_m^{1/2} and Re ~ Ra_m^{1/2} scalings remain unchanged within statistical uncertainty, indicating that the regime is not an artifact of under-resolution or inadequate sampling. revision: yes

  2. Referee: [Discussion of the mechanism] In two-dimensional turbulence the inverse energy cascade produces large-scale coherent structures that can dominate bulk transport even in the absence of the proposed magnetic-plume mechanism. The paper does not present comparisons with purely thermal (non-magnetic) 2D runs at the same Ra or with 3D simulations to isolate the magnetic force contribution from 2D condensation effects; such tests are load-bearing for the attribution of the scaling to magnetic plume advection.

    Authors: We acknowledge that 2D inverse cascades can generate large-scale structures, and we have revised the Discussion section to explicitly address this. We now include a direct comparison with non-magnetic 2D Rayleigh-Bénard convection at equivalent Ra (drawn from our own auxiliary simulations and literature data), demonstrating that the ultimate-like scaling does not appear in the purely thermal case at these parameters; the magnetic body force is required to sustain the efficient plume ejection and cross-bulk advection observed in our thermomagnetic runs. Additional visualizations and conditional statistics in the revised manuscript further isolate the magnetic force contribution by showing enhanced plume detachment correlated with local Lorentz force peaks. However, systematic 3D simulations at the highest Ra_m remain computationally prohibitive and are noted as a limitation of the present study; we argue that the high-Pr regime and the specific form of the magnetic force make the 2D results representative of the dominant transport mechanism. revision: partial

standing simulated objections not resolved
  • Full 3D DNS across the entire Ra_m range to definitively separate 2D condensation effects from the magnetic-plume mechanism is not feasible with current resources.

Circularity Check

0 steps flagged

No circularity: scaling obtained from DNS data then interpreted theoretically without reduction to fitted inputs or self-citations by construction

full rationale

The paper states that scaling laws are obtained through direct numerical simulations and formulated via theoretical arguments informed by the numerical data. The ultimate-like regime Nu ∼ Ra_m^{1/2} and Re ∼ Ra_m^{1/2} is reported as emerging from the simulations after transition, with the magnetic-plume mechanism offered as an interpretation of the observed behavior. No load-bearing equation reduces the reported exponents to a parameter fitted from the same data by construction, nor does any step invoke a self-citation chain or uniqueness theorem that collapses the central claim. The derivation remains self-contained against external benchmarks because the primary evidence is the DNS output itself rather than a tautological re-expression of inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on direct numerical simulations of a high-Prandtl-number fluid and on the interpretation that magnetic forces control plume ejection. No explicit free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption The flow is two-dimensional and the fluid has high Prandtl number
    Stated directly in the abstract as the setup for the simulations.

pith-pipeline@v0.9.0 · 5382 in / 1288 out tokens · 54805 ms · 2026-05-10T16:37:13.784831+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

37 extracted references · 37 canonical work pages

  1. [1]

    M. D. Cowley and R. E. Rosensweig, The interfacial sta- bilityof a ferromagnetic fluid, Journal ofFluidMechanics 30, 671–688 (1967)

    work page 1967

  2. [2]

    Bénard, Les tourbillons cellulaires dans une nappe liquide.-méthodes optiques d’observation et d’enregistrement, Journal de Physique Théorique et Ap- pliquée10, 254 (1901)

    H. Bénard, Les tourbillons cellulaires dans une nappe liquide.-méthodes optiques d’observation et d’enregistrement, Journal de Physique Théorique et Ap- pliquée10, 254 (1901)

    work page 1901

  3. [3]

    L. Rayleigh, On convection currents in a horizontal layer of fluid, when the higher temperature is on the under side, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science32, 529 (1916)

    work page 1916

  4. [4]

    B. A. Finlayson, Convective instability of ferromagnetic fluids, Journal of Fluid Mechanics40, 753 (1970)

    work page 1970

  5. [5]

    Yamaguchi, I

    H. Yamaguchi, I. Kobori, and Y. Uehata, Heat transfer in natural convection of magnetic fluids, Journal of Ther- mophysics and Heat Transfer13, 501 (1999)

    work page 1999

  6. [6]

    Krakov and I

    M. Krakov and I. Nikiforov, To the influence of uniform magnetic field on thermomagnetic convection in square cavity, Journal of Magnetism and Magnetic Materials 252, 209 (2002), proceedings of the 9th International Conference on Magnetic Fluids

    work page 2002

  7. [7]

    Yamaguchi, Z

    H. Yamaguchi, Z. Zhang, S. Shuchi, and K. Shimada, Heat transfer characteristics of magnetic fluid in a parti- tioned rectangular box, Journal of Magnetism and Mag- netic Materials252, 203 (2002), proceedings of the 9th International Conference on Magnetic Fluids

    work page 2002

  8. [8]

    Ganguly, S

    R. Ganguly, S. Sen, and I. K. Puri, Heat transfer aug- mentation using a magnetic fluid under the influence of a line dipole, Journal of Magnetism and Magnetic Mate- rials271, 63 (2004)

    work page 2004

  9. [9]

    H. L. G. Couto, N. B. Marcelino, and F. R. Cunha, A study on magnetic convection in a narrow rectangular cavity, Magnetohydrodynamics43, 421 (2007)

    work page 2007

  10. [10]

    Engler and S

    H. Engler and S. Odenbach, Parametric modulation of thermomagnetic convection in magnetic fluids, Journal of Physics: Condensed Matter20, 204135 (2008)

    work page 2008

  11. [11]

    P. S. Szabo and W.-G. Früh, The transition from natural convection to thermomagnetic convection of a magnetic fluid in a non-uniform magnetic field, Journal of Mag- netism and Magnetic Materials447, 116 (2018)

    work page 2018

  12. [12]

    Alegretti and R

    C. Alegretti and R. Gontijo, An experimental investiga- tionofthermomagneticconvectioninatallenclosuresub- jected to progressive field gradients, International Com- munications in Heat and Mass Transfer158, 107846 (2024)

    work page 2024

  13. [13]

    Oberbeck, Ueber die wärmeleitung der flüssigkeiten bei berücksichtigung der strömungen infolge von temper- aturdifferenzen, Annalen der Physik243, 271 (1879)

    A. Oberbeck, Ueber die wärmeleitung der flüssigkeiten bei berücksichtigung der strömungen infolge von temper- aturdifferenzen, Annalen der Physik243, 271 (1879)

    work page

  14. [14]

    J. Boussinesq,Théorie analytique de la chaleur mise en harmonie avec la thermodynamique et avec la théorie mé- canique de la lumière, Cours de physique mathématique de la Faculté des sciences (Gauthier-Villars, Paris, 1901)

    work page 1901

  15. [15]

    W. V. R. Malkus, The Heat Transport and Spectrum of Thermal Turbulence, Proceedings of the Royal Society of London Series A225, 196 (1954)

    work page 1954

  16. [16]

    C. H. B. Priestley, Convection from a Large Horizontal Surface, Australian Journal of Physics7, 176 (1954)

    work page 1954

  17. [17]

    Chandrasekhar,Hydrodynamic and hydromagnetic stability(Courier Corporation, 2013)

    S. Chandrasekhar,Hydrodynamic and hydromagnetic stability(Courier Corporation, 2013)

    work page 2013

  18. [18]

    R. H. Kraichnan, Turbulent thermal convection at arbi- trary Prandtl number, Physics of Fluids5, 1374 (1962)

    work page 1962

  19. [19]

    B. I. Shraiman and E. D. Siggia, Heat transport in high-rayleigh-number convection, Phys. Rev. A42, 3650 (1990)

    work page 1990

  20. [20]

    Grossmann and D

    S. Grossmann and D. Lohse, Scaling in thermal convec- tion: a unifying theory, Journal of Fluid Mechanics407, 27–56 (2000)

    work page 2000

  21. [21]

    Castaing, G

    B. Castaing, G. Gunaratne, F. Heslot, L. Kadanoff, A. Libchaber, S. Thomae, X.-Z. Wu, S. Zaleski, and G. Zanetti, Scaling of hard thermal turbulence in rayleigh-bénard convection, Journal of Fluid Mechanics 6 204, 1–30 (1989)

    work page 1989

  22. [22]

    Chavanne, F

    X. Chavanne, F. Chillà, B. Castaing, B. Hébral, B. Chabaud, and J. Chaussy, Observation of the ulti- mate regime in rayleigh-bénard convection, Phys. Rev. Lett.79, 3648 (1997)

    work page 1997

  23. [23]

    Johnston and C

    H. Johnston and C. R. Doering, Comparison of turbulent thermal convection between conditions of constant tem- perature and constant flux, Phys. Rev. Lett.102, 064501 (2009)

    work page 2009

  24. [24]

    Shishkina and C

    O. Shishkina and C. Wagner, Analysis of thermal dis- sipation rates in turbulent rayleigh–bénard convection, Journal of Fluid Mechanics546, 51–60 (2006)

    work page 2006

  25. [25]

    E. P. van der Poel, R. J. A. M. Stevens, and D. Lohse, Comparison between two- and three- dimensional Rayleigh–Bénard convection, Journal of Fluid Mechanics736, 177 (2013)

    work page 2013

  26. [26]

    X. Zhu, V. Mathai, R. J. A. M. Stevens, R. Verzicco, and D. Lohse, Transition to the ultimate regime in two- dimensional rayleigh-bénard convection, Phys. Rev. Lett. 120, 144502 (2018)

    work page 2018

  27. [27]

    Ahlers, S

    G. Ahlers, S. Grossmann, and D. Lohse, Heat transfer and large scale dynamics in turbulent rayleigh-bénard convection, Rev. Mod. Phys.81, 503 (2009)

    work page 2009

  28. [28]

    D.LohseandO.Shishkina,Ultimaterayleigh-bénardtur- bulence, Rev. Mod. Phys.96, 035001 (2024)

    work page 2024

  29. [29]

    Gotoh and M

    K. Gotoh and M. Yamada, Thermal convection in a hor- izontal layer of magnetic fluids, Journal of Physical Soci- ety of Japan51, 3042 (1982)

    work page 1982

  30. [30]

    Capobianchi, M

    P. Capobianchi, M. Lappa, M. S. N. Oliveira, and F. T. Pinho, Shear rheology of a dilute emulsion of ferrofluid droplets dispersed in a nonmagnetizable carrier fluid un- der the influence of a uniform magnetic field, Journal of Rheology65, 925 (2021)

    work page 2021

  31. [31]

    J. L. Neuringer and R. E. Rosensweig, Ferrohydrodynam- ics, The Physics of Fluids7, 1927 (1964)

    work page 1927

  32. [32]

    S. A. Suslov, Thermomagnetic convection in a vertical layer of ferromagnetic fluid, Physics of Fluids20, 084101 (2008)

    work page 2008

  33. [33]

    Capobianchi, See Supplemental Material at [URL] for technical details

    P. Capobianchi, See Supplemental Material at [URL] for technical details

    work page

  34. [34]

    Heckert, L

    M. Heckert, L. Sprenger, A. Lange, and S. Odenbach, Ex- perimental determination of the critical rayleigh number for thermomagnetic convection with focus on fluid com- position, Journal of Magnetism and Magnetic Materials 381, 337 (2015)

    work page 2015

  35. [35]

    T. P. Bednarz, W. Lin, J. C. Patterson, C. Lei, and S. W. Armfield, Scaling for unsteady thermo-magnetic convec- tion boundary layer of paramagnetic fluids of pr>1 in micro-gravity conditions, International Journal of Heat and Fluid Flow30, 1157 (2009)

    work page 2009

  36. [36]

    Chavanne, F

    X. Chavanne, F. Chillà, B. Chabaud, B. Castaing, and B. Hébral, Turbulent rayleigh–bénard convection in gaseous and liquid he, Physics of Fluids13, 1300 (2001)

    work page 2001

  37. [37]

    X. He, D. Funfschilling, H. Nobach, E. Bodenschatz, and G. Ahlers, Transition to the ultimate state of turbulent rayleigh-bénardconvection,Phys.Rev.Lett.108,024502 (2012)

    work page 2012