pith. sign in

arxiv: 2512.01308 · v4 · submitted 2025-12-01 · ⚛️ physics.flu-dyn

Analysis and reformulation of the k--ω turbulence model for buoyancy-driven thermal convection

Da-Sol Joo This is my paper

Pith reviewed 2026-05-17 03:28 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords k-omega turbulence modelbuoyancy-driven convectionRayleigh-Benard convectionRANS modelingNusselt number scalingPrandtl number effectsnatural convectionheat transfer modeling
0
0 comments X

The pith

Two algebraic functions added to the k-ω model restore the observed heat transfer scaling with Prandtl number in buoyancy-driven convection.

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

The paper derives an analytical solution of the standard k-ω turbulence model for Rayleigh-Bénard convection in an infinite layer where buoyancy alone generates turbulence. This yields explicit scaling relations for the Nusselt number that differ from experimental trends across low and high Prandtl numbers. Informed by the mismatches, the authors introduce two dimensionless algebraic functions into the buoyancy source terms. The functions are constructed to vanish when buoyancy is absent so the model stays identical to the standard closure in non-buoyant flows. The reformulated model is tested on multiple buoyancy-driven configurations and produces closer agreement with measured mean temperatures and turbulent heat fluxes.

Core claim

Analysis of the standard k-ω equations for pure-buoyancy Rayleigh-Bénard convection in an infinite layer produces the scalings Nu ∼ Ra^{1/3} Pr^{1/3} for Pr ≪ 1 and Nu ∼ Ra^{1/3} Pr^{-0.415} for Pr ≫ 1. Reformulation of the buoyancy modeling terms with two algebraic functions recovers the measured trends Nu ∼ Pr^{1/8} for Pr ≪ 1 and Nu ∼ Pr^0 for Pr ≫ 1 at moderate Ra while preserving full compatibility with the standard closure.

What carries the argument

Two dimensionless algebraic functions inserted into the buoyancy production terms of the k and ω equations.

If this is right

  • Mean temperature fields improve in two-dimensional Rayleigh-Bénard convection.
  • Turbulent heat flux distributions improve in internally heated convection and unstably stratified Couette flow.
  • Predictions remain consistent with the original model in flows without buoyancy.
  • Results hold across different aspect ratios in vertically heated natural convection.

Where Pith is reading between the lines

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

  • The same functions could be tested in other two-equation RANS models to handle unstable stratification.
  • Adoption might reduce case-by-case tuning for engineering calculations of natural convection in enclosures or channels.
  • Extension to three-dimensional or higher-Rayleigh-number regimes would test whether the recovered scalings continue to hold.

Load-bearing premise

The two algebraic functions chosen to match infinite-layer trends remain accurate without adjustment when applied to the full range of tested buoyancy-driven flows.

What would settle it

An experimental measurement of Nusselt number versus Prandtl number at moderate Rayleigh number in Rayleigh-Bénard convection that deviates from the target scalings of Pr to the power 1/8 at low Pr and constant at high Pr.

Figures

Figures reproduced from arXiv: 2512.01308 by Da-Sol Joo.

Figure 1
Figure 1. Figure 1: FIG. 1: (a) Rayleigh–B´enard convection setup, (b) mean temperature distribution, and (c) [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: One-dimensional Rayleigh–B´enard convection predictions of the [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: One-dimensional Rayleigh–B´enard convection predictions of the [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: One-dimensional Rayleigh–B´enard convection predictions of the [PITH_FULL_IMAGE:figures/full_fig_p020_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: (a) Dependence of [PITH_FULL_IMAGE:figures/full_fig_p024_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Analytical distributions at [PITH_FULL_IMAGE:figures/full_fig_p026_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Simulation results for one-dimensional Rayleigh–B´enard convection with the [PITH_FULL_IMAGE:figures/full_fig_p027_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Stably stratified channel flow at [PITH_FULL_IMAGE:figures/full_fig_p030_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Internally heated convection setup for (a) top cooling and (b) top-and-bottom [PITH_FULL_IMAGE:figures/full_fig_p032_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Internally heated convection with top cooling at [PITH_FULL_IMAGE:figures/full_fig_p033_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: Internally heated convection with top-and-bottom cooling. Circles denote the [PITH_FULL_IMAGE:figures/full_fig_p034_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Simulation results for internally heated convection with top-and-bottom cooling [PITH_FULL_IMAGE:figures/full_fig_p036_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: Unstably stratified Couette flow simulation results of the corrected [PITH_FULL_IMAGE:figures/full_fig_p037_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: Rayleigh–B´enard convection simulation results of the corrected [PITH_FULL_IMAGE:figures/full_fig_p039_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15: Simulation results from the corrected [PITH_FULL_IMAGE:figures/full_fig_p040_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16: Natural convection in a rectangular cavity with differentially heated vertical [PITH_FULL_IMAGE:figures/full_fig_p041_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17: Simulation results of vertically heated natural convection for [PITH_FULL_IMAGE:figures/full_fig_p042_17.png] view at source ↗
read the original abstract

The representation of buoyancy-driven turbulence in Reynolds-averaged Navier--Stokes (RANS) models remains unresolved, with no widely accepted standard formulation. A key difficulty is the lack of analytical guidance for incorporating buoyant effects, particularly under unstable stratification. This study derives an analytical solution of the standard $k$--$\omega$ model for Rayleigh--B\'enard convection in an infinite layer, where turbulent kinetic energy is generated solely by buoyancy. The solution provides explicit scaling relations among the Rayleigh ($Ra$), Prandtl ($Pr$), and Nusselt ($Nu$) numbers that capture the simulation trends: $Nu \sim Ra^{1/3} Pr^{1/3}$ for $Pr \ll 1$ and $Nu \sim Ra^{1/3} Pr^{-0.415}$ for $Pr \gg 1$. This framework quantifies the discrepancies in the conventional buoyancy treatment and clarifies their origin. Informed by this analysis, the buoyancy-related modeling terms are reformulated to recover the measured trends: namely $Nu \sim Pr^{1/8}$ for $Pr \ll 1$ and $Nu \sim Pr^{0}$ for $Pr \gg 1$ at moderate $Ra$. Only two dimensionless algebraic functions are introduced, which vanish in the absence of buoyancy, ensuring full compatibility with the standard closure. The corrected model is validated across a range of buoyancy-driven flows, including two-dimensional Rayleigh--B\'enard convection, internally heated convection in two configurations, unstably stratified Couette flow, and vertically heated natural convection with varying aspect ratios. Across all cases, the corrected model provides significantly improved predictions of mean temperature fields and turbulent heat flux distributions.

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 derives an analytical solution for the standard k–ω model applied to Rayleigh–Bénard convection in an infinite horizontal layer with buoyancy-only turbulent kinetic energy production. This yields explicit scaling relations (Nu ∼ Ra^{1/3} Pr^{1/3} for Pr ≪ 1 and Nu ∼ Ra^{1/3} Pr^{-0.415} for Pr ≫ 1) that differ from measured trends. Informed by the discrepancies, the buoyancy-related terms are reformulated by introducing two dimensionless algebraic functions chosen to recover the observed scalings Nu ∼ Pr^{1/8} (Pr ≪ 1) and Nu ∼ Pr^0 (Pr ≫ 1) at moderate Ra; these functions are constructed to vanish without buoyancy. The modified closure is then validated on 2D Rayleigh–Bénard convection, two internally heated configurations, unstably stratified Couette flow, and vertically heated natural convection at varying aspect ratios, with reported improvements in mean temperature and turbulent heat flux predictions.

Significance. The analytical solution for the infinite layer constitutes a clear strength, supplying explicit, internally consistent scaling relations that can serve as a diagnostic benchmark for buoyancy closures. If the two algebraic functions prove transferable, the reformulation would offer a practical, backward-compatible improvement for RANS modeling of buoyancy-driven flows, an area that lacks a widely accepted standard. The multi-configuration validation adds weight, though the empirical calibration of the functions to target trends reduces the a priori character of the advance.

major comments (2)
  1. [§4 (reformulation)] §4 (reformulation): The two dimensionless algebraic functions are introduced explicitly to enforce the target Nu–Pr scalings (Pr^{1/8} at low Pr, Pr^0 at high Pr) identified from the infinite-layer analysis. While they vanish in the absence of buoyancy, their functional forms are selected to recover the measured trends by construction rather than from additional physical constraints or a parameter-free derivation; this makes the central modeling change partly empirical and requires clearer justification that the correction is not simply a data-guided adjustment.
  2. [§5 (validation)] §5 (validation): The same fixed algebraic functions, calibrated solely on the infinite horizontal layer with pure buoyancy production, are applied without retuning or sensitivity tests to configurations that include mean shear (unstably stratified Couette flow) and altered buoyancy–shear interactions (internally heated and vertical enclosures). No alternative functional forms or parameter variations are examined to assess whether the dependence remains accurate when the production mechanisms and boundary conditions differ from the calibration case.
minor comments (2)
  1. [Abstract] Abstract: The claim of 'significantly improved predictions' is stated without quantitative metrics (e.g., relative error reductions in Nu or temperature profiles); adding brief numerical indicators would make the improvement concrete.
  2. [Notation] Notation: The exact mathematical expressions for the two algebraic functions would benefit from being isolated in a single, numbered equation with all arguments and limiting behaviors shown explicitly.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and for highlighting the value of the analytical solution as a diagnostic benchmark. We respond to each major comment below, providing additional context from the derivation while noting where further clarification will be added to the manuscript.

read point-by-point responses

  1. Referee: §4 (reformulation): The two dimensionless algebraic functions are introduced explicitly to enforce the target Nu–Pr scalings (Pr^{1/8} at low Pr, Pr^0 at high Pr) identified from the infinite-layer analysis. While they vanish in the absence of buoyancy, their functional forms are selected to recover the measured trends by construction rather than from additional physical constraints or a parameter-free derivation; this makes the central modeling change partly empirical and requires clearer justification that the correction is not simply a data-guided adjustment.

    Authors: The functional forms were chosen as the minimal algebraic expressions that restore the experimentally observed Nu–Pr exponents once the discrepancies identified by the infinite-layer solution are quantified. Because the analysis supplies explicit relations for the standard model, the corrections are directly tied to those relations rather than to a broad empirical fit. The requirement that both functions vanish identically when buoyancy production is zero further constrains their structure and guarantees compatibility with the original closure. We agree that the motivation can be stated more explicitly and will revise §4 to include a step-by-step mapping from the derived scaling discrepancies to the arguments and exponents of the two functions. revision: partial

  2. Referee: §5 (validation): The same fixed algebraic functions, calibrated solely on the infinite horizontal layer with pure buoyancy production, are applied without retuning or sensitivity tests to configurations that include mean shear (unstably stratified Couette flow) and altered buoyancy–shear interactions (internally heated and vertical enclosures). No alternative functional forms or parameter variations are examined to assess whether the dependence remains accurate when the production mechanisms and boundary conditions differ from the calibration case.

    Authors: The decision to keep the functions fixed across all test cases was intentional: it constitutes a direct test of transferability to flows in which shear and buoyancy production coexist or in which the buoyancy–shear balance differs from the calibration configuration. The reported improvements in mean temperature and turbulent heat flux without any retuning already provide evidence that the corrections are not narrowly tuned to the infinite-layer case. We acknowledge that a systematic exploration of alternative functional forms or parameter variations would strengthen the claim of robustness; the present manuscript does not contain such a study. In the revised version we will add a short discussion of the observed robustness and the limitations of the current functional dependence. revision: partial

Circularity Check

1 steps flagged

Reformulation introduces two algebraic functions calibrated to enforce target Nu-Pr scalings

specific steps
  1. fitted input called prediction [Abstract (reformulation paragraph)]
    "Informed by this analysis, the buoyancy-related modeling terms are reformulated to recover the measured trends: namely Nu ∼ Pr^{1/8} for Pr ≪ 1 and Nu ∼ Pr^0 for Pr ≫ 1 at moderate Ra. Only two dimensionless algebraic functions are introduced, which vanish in the absence of buoyancy, ensuring full compatibility with the standard closure."

    The two algebraic functions are introduced specifically to enforce the target Nu-Pr scalings that match measured data. Consequently, the claim that the reformulated model 'recovers the measured trends' is achieved by the choice of the functions rather than emerging as an independent prediction from the equations.

full rationale

The paper first derives explicit Ra-Pr-Nu scalings from the unmodified k-ω equations under buoyancy-only production in an infinite layer. It then introduces two dimensionless algebraic functions whose explicit purpose is to alter the buoyancy terms so that the model reproduces the experimentally observed Nu ∼ Pr^{1/8} (low Pr) and Nu ∼ Pr^0 (high Pr) behaviors. Because these functions are constructed and tuned precisely to recover the measured trends (while vanishing in the non-buoyant limit), the central modeling improvement reduces in part to an empirical adjustment whose success on the calibration geometry is guaranteed by construction. Subsequent validation on other geometries therefore tests transferability rather than an independent first-principles derivation. This produces moderate circularity (score 6) without rendering the entire analysis tautological.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard k-ω closure equations plus the choice of two algebraic functions tuned to observed trends; no new physical entities are postulated.

free parameters (1)
  • two dimensionless algebraic functions for buoyancy terms
    Introduced explicitly to recover measured Nusselt-Prandtl trends in the reformulation step.
axioms (1)
  • domain assumption Standard k-ω model equations remain valid when buoyancy is the sole turbulence production mechanism
    The analytical solution begins from the unmodified model equations applied to the infinite-layer Rayleigh-Bénard setup.

pith-pipeline@v0.9.0 · 5612 in / 1397 out tokens · 50999 ms · 2026-05-17T03:28:38.460981+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

3 extracted references · 3 canonical work pages

  1. [1]

    Computational modeling of turbulent flows,

    1J. L. Lumley, “Computational modeling of turbulent flows,” Adv. Appl. Mech.18, 123–176 (1979). 2W. Rodi,Turbulence Model and Their Application in Hydraulics: A state of the Art Review (Int. Ass. for Hydraulic Research, Delft, The Netherlands, 1980). 3W. Rodi, “Examples of calculation methods for flow and mixing in stratified fluids,” J. Geophys. Res. Oce...

    work page 1979

  2. [2]

    On the accuracy of RANS simulations of 2D boundary layers with OpenFOAM,

    35S. Gomez, B. Graves, and S. Poroseva, “On the accuracy of RANS simulations of 2D boundary layers with OpenFOAM,” in44th AIAA Fluid Dynamics Conference(2014) p

    work page 2014

  3. [3]

    Prandtl number effects in convective turbulence,

    36R. Verzicco and R. Camussi, “Prandtl number effects in convective turbulence,” J. Fluid Mech.383, 55–73 (1999). 37A. Pandey, D. Krasnov, K. R. Sreenivasan, and J. Schumacher, “Convective mesoscale turbulence at very low Prandtl numbers,” J. Fluid Mech.948, A23 (2022). 38X.-M. Li, J.-D. He, Y. Tian, P. Hao, and S.-D. Huang, “Effects of Prandtl number in ...

    work page doi:10.5281/zenodo.17101542(2025 1999