pith. sign in

arxiv: 2604.25362 · v2 · submitted 2026-04-28 · ⚛️ physics.flu-dyn · math-ph· math.MP

The Wooding problem revisited

A. Barletta , D. A. S. Rees This is my paper

Pith reviewed 2026-05-07 14:59 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn math-phmath.MP
keywords porous mediaconvective instabilityBiot numberRayleigh numberWooding problemlinear stabilitysuction boundary layergeothermal flow
0
0 comments X

The pith

The heat-flux-based Rayleigh number has a finite limit at neutral stability for both infinite and zero Biot numbers.

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

The paper reconsiders the classic Wooding setup for convective instability onset in a semi-infinite porous layer with steady suction at the boundary. It introduces a Biot number to account for imperfect heat transfer across the top surface, recovering the original Wooding temperature condition only in the infinite-Biot limit. Two Rayleigh-number definitions are compared: the usual temperature-difference version and a heat-flux version. The heat-flux version remains finite at the stability threshold in both the perfect-heat-transfer and no-heat-transfer extremes, while the temperature version diverges as the Biot number approaches zero. A reader would care because this supplies a consistent nondimensional measure for predicting when convection begins under a wider range of realistic boundary conditions in geothermal and porous-media flows.

Core claim

Linear stability analysis of the stationary suction-driven boundary layer shows that the heat-flux-based Rayleigh number attains finite limiting values at the onset of convective instability as the Biot number tends to either zero or infinity, while the temperature-difference-based Rayleigh number diverges in the zero-Biot limit.

What carries the argument

The heat-flux-based Rayleigh number, formed from the imposed boundary heat flux rather than a temperature difference, which remains bounded at neutral stability across the full range of Biot numbers.

Load-bearing premise

The analysis assumes the standard Darcy-Boussinesq equations hold and that a steady suction-driven base state exists whose linear stability can be examined.

What would settle it

A laboratory experiment or numerical simulation that measures the critical suction velocity or heat flux at which convection onsets in a porous layer for both very large and very small Biot numbers would directly test whether the heat-flux Rayleigh number remains finite at those limits.

Figures

Figures reproduced from arXiv: 2604.25362 by A. Barletta, D. A. S. Rees.

Figure 1
Figure 1. Figure 1: A sketch of the semi-infinite porous medium and its boundary conditions. view at source ↗
Figure 2
Figure 2. Figure 2: Neutral stability values of Ra versus ymax for k = 0.5 (a) and for k = 2 (b), with different Biot numbers. [18]. A description of the shooting method applied to a similar problem, i.e. the Wooding problem for an inclined boundary, is also available in Rees and Bassom [5]. As pointed out by these authors [5], a special feature of the Wooding problem is that it is relative to a semi-infinite domain. Then, un… view at source ↗
Figure 3
Figure 3. Figure 3: Neutral stability curves in the (k, Ra) plane (a) and in the (k, Rm) plane (b), with different Biot numbers. On the other hand, with Bi → 0, linear instability is possible provided that the value of Rm is suffi￾ciently large. The reason of this dissimilarity is that Rm is a Rayleigh number based on the dimensional temperature difference ∆Tm, defined by (24). As already mentioned in Section 3.4, ∆Tm is a he… view at source ↗
Figure 4
Figure 4. Figure 4: Critical values of k, Ra and Rm versus Bi. The dashed lines (red and magenta) show the small-Bi and large-Bi approximations given by (33) and (40). together with the functions Ψ1(y), Θˆ 1(y), Ψk1(y) and Θˆ k1(y). Obviously, one may develop also the 2nd - order differential problem and so on. The linear approximations of kc and Rmc when Bi ≪ 1 yield kc ≈ 0.709207 + 0.072203 Bi, Rmc ≈ 10.492375 + 2.215147 Bi… view at source ↗
Figure 5
Figure 5. Figure 5: Perturbation streamlines (black) and isotherms (turquoise) at critical conditions ( view at source ↗
read the original abstract

The threshold conditions to convective instability in a semi-infinite porous layer saturated by a fluid are determined. The classical setup for this problem in geothermal fluid dynamics was originally modelled by Wooding in 1960. Its formulation is here reconsidered to allow for an imperfect heat transfer across the boundary, parametrised through the Biot number. The temperature boundary condition considered by Wooding is here recovered as the limit of an infinite Biot number. The linear stability analysis of the stationary boundary layer which establishes in the porous medium when a boundary steady suction occurs is carried out. Two different versions of the Rayleigh number are considered, namely, a temperature-difference-based version and a heat-flux-based version. While the former is the classical Rayleigh number for flow in porous media, the latter is a variant definition which displays a finite limit at neutral stability in both the opposite limiting cases of an infinite or of a zero Biot number.

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

0 major / 2 minor

Summary. The manuscript revisits the Wooding problem of convective instability in a semi-infinite porous layer with steady suction, extending the classical setup to a Robin thermal boundary condition parametrized by the Biot number Bi. Linear stability analysis is performed on the resulting steady base state. Two Rayleigh numbers are defined and analyzed: the conventional temperature-difference-based Rayleigh number and a heat-flux-based variant. The central claim is that the critical value of the heat-flux-based Rayleigh number remains finite at neutral stability in both the Bi → ∞ limit (recovering the original Wooding case) and the Bi → 0 limit.

Significance. If the linear stability results are correct, the work supplies a useful alternative Rayleigh number definition that remains well-behaved across the full range of Biot numbers. This is relevant to geothermal and porous-media convection modeling where boundary heat transfer is imperfect. The analysis stays within the standard Darcy-Boussinesq framework and linear stability theory. The stress-test concern that the Bi → 0 limit violates the Boussinesq approximation does not land: the paper presents the finite limit as a mathematical property of the nondimensionalization inside the model and makes no claim that the extreme limit is physically realizable with finite heat flux.

minor comments (2)
  1. The abstract states the finite-limit claim for the heat-flux-based Rayleigh number but supplies no explicit base-state solution, perturbation equations, or eigenvalue problem, making independent assessment of the result difficult from the summary alone.
  2. The introduction or methods section should explicitly contrast the definitions of the two Rayleigh numbers (including how the heat-flux version incorporates the realized temperature jump C) to clarify why only one remains finite as Bi → 0.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive and constructive assessment of our manuscript on the finite-Biot-number generalization of the Wooding problem. The recommendation of minor revision is noted; we will incorporate any editorial or presentational improvements in the revised version. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No circularity: critical values obtained from independent linear stability eigenvalue problem

full rationale

The paper performs a standard linear stability analysis of the steady suction base flow under a Robin thermal boundary condition parameterized by the Biot number Bi. The temperature-based Rayleigh number is the classical definition Ra_T = g β K ΔT δ / (ν α), while the heat-flux-based version Ra_q is introduced as a rescaled variant using the imposed flux q. Critical thresholds are computed by solving the resulting eigenvalue problem for neutral stability; the Bi-dependence of these thresholds is obtained directly from the dispersion relation and boundary conditions, approaching finite limits as Bi → 0 and Bi → ∞ without any redefinition of the target quantity in terms of the output, without parameter fitting to data, and without load-bearing self-citations. The derivation chain is self-contained within the Darcy-Boussinesq model and the stated base-state assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract relies on standard Darcy-Boussinesq porous-media equations and a steady base state but lists no explicit free parameters or new entities.

axioms (2)
  • domain assumption Darcy-Boussinesq equations govern flow and heat transport
    Implicit in reference to the Wooding setup for porous convection.
  • domain assumption Steady suction-driven base state exists and is linearly stable below threshold
    Required for the linear stability analysis described.

pith-pipeline@v0.9.0 · 8185 in / 1116 out tokens · 93130 ms · 2026-05-07T14:59:01.122714+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

19 extracted references · 19 canonical work pages

  1. [1]

    R. A. Wooding, Rayleigh instability of a thermal boundary layer in flow through a porous medium, Journal of Fluid Mechanics 9 (1960) 183–192

    work page 1960

  2. [2]

    C. J. Van Duijn, R. A. Wooding, G. J. M. Pieters, A. Van Der Ploeg, Stability criteria for the boundary layer formed by throughflow at a horizontal surface of a porous medium, Technical Report, Reports on Applied and Numerical Analysis, Department of Mathematics and Computing Science, Eindhoven University of Technology, Netherlands, 2001. 15

    work page 2001

  3. [3]

    D. A. S. Rees, The onset and nonlinear development of vortex instabilities in a horizontal forced convection boundary layer with uniform surface suction, Transport in Porous Media 77 (2009) 243–265

    work page 2009

  4. [4]

    P. M. Patil, D. A. S. Rees, Linear instability of a horizontal thermal boundary layer formed by vertical throughflow in a porous medium: the effect of local thermal nonequilibrium, Transport in Porous Media 99 (2013) 207–227

    work page 2013

  5. [5]

    D. A. S. Rees, A. P. Bassom, The inclined Wooding problem, Transport in Porous Media 125 (2018) 465–482

    work page 2018

  6. [6]

    Islam, M

    S. Islam, M. Ferdows, D. A. S. Rees, A. P. Bassom, The onset of convection in the Wooding prob- lem with anisotropic permeability, Special Topics & Reviews in Porous Media: An International Journal 15 (2024) 61–75

    work page 2024

  7. [7]

    F. M. Sutton, Onset of convection in a porous channel with net through flow, Physics of Fluids 13 (1970) 1931–1934

    work page 1970

  8. [8]

    G. M. Homsy, A. E. Sherwood, Convective instabilities in porous media with through flow, AIChE Journal 22 (1976) 168–174

    work page 1976

  9. [9]

    M. C. Jones, J. M. Persichetti, Convective instability in packed beds with throughflow, AIChE journal 32 (1986) 1555–1557

    work page 1986

  10. [10]

    Capone, J

    F. Capone, J. A. Gianfrani, G. Massa, D. A. S. Rees, A weakly nonlinear analysis of the effect of vertical throughflow on Darcy-Bénard convection, Physics of Fluids 35 (2023) 014107

    work page 2023

  11. [11]

    J. A. Gianfrani, D. A. S. Rees, Eckhaus instability in the Sutton problem: Dynamics of two interacting modes, Physics of Fluids 37 (2025) 054114

    work page 2025

  12. [12]

    Deepika, P

    N. Deepika, P. A. L. Narayana, A. A. Hill, Onset of Darcy-Brinkman convection with a uniform internal heat source and vertical throughflow, International Journal of Thermal Sciences 117 (2017) 136–144

    work page 2017

  13. [13]

    Turkyilmazoglu, Darcy-Bénard convection through a uniformly permeable porous slab, Com- munications in Nonlinear Science and Numerical Simulation 125 (2023) 107392

    M. Turkyilmazoglu, Darcy-Bénard convection through a uniformly permeable porous slab, Com- munications in Nonlinear Science and Numerical Simulation 125 (2023) 107392

    work page 2023

  14. [14]

    Capone, R

    F. Capone, R. De Luca, G. Massa, Throughflow effect on bi-disperse convection, Ricerche di Matematica 73 (2024) 67–84

    work page 2024

  15. [15]

    Hitchen, A

    J. Hitchen, A. J. Wells, The impact of imperfect heat transfer on the convective instability of a thermal boundary layer in a porous media, Journal of Fluid Mechanics 794 (2016) 154–174

    work page 2016

  16. [16]

    D. A. Nield, A. Bejan, Convection in Porous Media, 5th ed., Springer, New York, 2017

    work page 2017

  17. [17]

    Straughan, Stability and Wave Motion in Porous Media, Springer, New York, 2008

    B. Straughan, Stability and Wave Motion in Porous Media, Springer, New York, 2008

    work page 2008

  18. [18]

    Barletta, Routes to Absolute Instability in Porous Media, Springer, New York, 2019

    A. Barletta, Routes to Absolute Instability in Porous Media, Springer, New York, 2019

    work page 2019

  19. [19]

    J. W. Eaton, D. Bateman, S. Hauberg, R. Wehbring, GNU Octave, A high-level interactive language for numerical computations, 11th ed.,https://docs.octave.org/octave.pdf, 2026. 16

    work page 2026