Hostname: page-component-7c8c6479df-7qhmt Total loading time: 0 Render date: 2024-03-28T11:26:54.408Z Has data issue: false hasContentIssue false

High-Rayleigh-number convection in a fluid-saturated porous layer

Published online by Cambridge University Press:  03 February 2004

JESSE OTERO
Affiliation:
Department of Mathematics, Ohio State University, Columbus, OH 43210, USA
LUBOMIRA A. DONTCHEVA
Affiliation:
Department of Computer Science & Engineering, University of Washington, Seattle, WA 98195-2350, USA
HANS JOHNSTON
Affiliation:
Department of Mathematics, University of Michigan, Ann Arbor, MI 48109-1109, USA
RODNEY A. WORTHING
Affiliation:
Breasco LLC, Ann Arbor, MI 48197, USA
ALEXANDER KURGANOV
Affiliation:
Department of Mathematics, Tulane University, New Orleans, LA 70118, USA
GUERGANA PETROVA
Affiliation:
Department of Mathematics, Texas A&M University, College Station, TX 77843-3368, USA
CHARLES R. DOERING
Affiliation:
Department of Mathematics, University of Michigan, Ann Arbor, MI 48109-1109, USA Michigan Center for Theoretical Physics, Ann Arbor, MI 48109-1120, USA

Abstract

The Darcy–Boussinesq equations at infinite Darcy–Prandtl number are used to study convection and heat transport in a basic model of porous-medium convection over a broad range of Rayleigh number $Ra$. High-resolution direct numerical simulations are performed to explore the modes of convection and measure the heat transport, i.e. the Nusselt number Nu, from onset at $Ra \,{=}\, 4\pi^2$ up to $Ra\,{=}\,10^4$. Over an intermediate range of increasing Rayleigh numbers, the simulations display the ‘classical’ heat transport $\hbox{\it Nu} \,{\sim}\, Ra$ scaling. As the Rayleigh number is increased beyond $Ra \,{=}\, 1255$, we observe a sharp crossover to a form fitted by $\hbox{\it Nu} \,{\approx}\, 0.0174 \times Ra^{0.9}$ over nearly a decade up to the highest $Ra$. New rigorous upper bounds on the high-Rayleigh-number heat transport are derived, quantitatively improving the most recent available results. The upper bounds are of the classical scaling form with an explicit prefactor: $\hbox{\it Nu} \,{\le}\, 0.0297 \times Ra$. The bounds are compared directly to the results of the simulations. We also report various dynamical transitions for intermediate values of $Ra$, including hysteretic effects observed in the simulations as the Rayleigh number is decreased from $1255$ back down to onset.

Type
Papers
Copyright
© 2004 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)