**Musterbildung, Energielandschaften und Skalierungsgesetze**

**Direktor: **

Felix Otto

**Kontakt: **Email

*Telefon:*

+49 (0) 341 - 9959 - 950

*Fax:*

+49 (0) 341 - 9959 - 658

**Anschrift:**

Inselstr. 22

04103 Leipzig

**Sekretariat:**

Katja Heid

Email, *Telefon:*

+49 (0) 341 - 9959 - 951

# Rayleigh-Bénard Convection

## A model of thermal convection

**Rayleigh-Bénard convection** (RBC) is the buoyancy-driven ﬂow of a ﬂuid heated from below and cooled from above.
This model of thermal convection is a paradigm for nonlinear and chaotic dynamics, pattern formation and fully developed
turbulence (Kadanoﬀ 2001 [1]). RBC plays an important
role in a large range of phenomena in geophysics, astrophysics, meteorology,
oceanography and engineering.
The problem under investigation is:

Given an incompressible ﬂuid enclosed in a container heated from below and
cooled from above, what are the ﬂow dynamics? In particular, what is the
heat transfer from the bottom to the top?

The mechanisms that determine the dynamics are:

**Thermal diﬀusion:**heat flux, due to the temperature gradient in the box.**Buoyancy-driven convection:**bulk fluid movement due to temperature differences which affect the viscosity and thus the relative buoyancy.**Inner friction:**force that resists the deformation of the fluid parcels. It's effect is very strong near the plates, where the ﬂuid is at rest.

In the numerical simulation below, the ﬂuid is at rest in the bulk until **convection rolls** appear which in turn are destabilzed by
ﬁnger-like structures that detach from the boundary layers, called plumes.
Eventually these formations trigger the breakdown of the convection rolls and the
motion becomes chaotic/turbulent.

*Video by M. Zimmermann. Click to play.*

In the Oberbeck (1879) -Boussinesq (1903) approximation, where the density \(\rho\) is assumed to depend linearly on the temperature \(T\) , the dimensionless equations of motion are :

Figure 3: The geometry of the setup and boundary conditions. Periodic sidewall conditions are imposed.

\( \begin{cases} \frac{\partial }{\partial t}T+\mathbf{u}\cdot \nabla T=\Delta T & (1)\\ \frac{1}{Pr}\left(\frac{\partial }{\partial t}\mathbf{u}+\mathbf{u}\cdot \nabla \mathbf{u}\right)-\Delta \mathbf{u}+\nabla p=Ra T\mathbf{e_z} & (2)\\ \nabla\cdot \mathbf{u}=0 & (3)\\ \end{cases} \)

where \(\mathbf{u}\) is the velocity ﬁeld and \(p\) is the pressure.
On the two plates at height \(z = 0\) and \(z = 1\), respectively, the velocity ﬁeld
satisﬁes *no-slip* boundary conditions and the temperature is \(1\) at \(z = 0\) and
\(0\) at \(z = 1\).

All the quantities \((\textbf{u},T,p)\) are assumed to be periodic in the the horizontal
variables \((x,y)\).

From the non-dimensionalization, only two control parameters are left:

the Prandtl number, \(Pr\), and the Rayleigh number, \(Ra\).

By deﬁnition

\(Pr=\frac{\nu}{\kappa}\)

where \(\nu\) is the fluid's kinematic viscosity and \(\kappa\) is the thermal diffusivity and

\(Ra\propto (T_{hot}-T_{cold}) H^3\)

where \(H\) is the hight of the box.

## The Nusselt number

In the applications, it is of interest to measure the effectiveness of the motion.
The most natural and accepted measure that quantifies the enhancement of vertical heat flux due to convection
is the **Nusselt number**, \(Nu\).

Integrating in \(z\in (0,1)\) the long-time and horizontal-space averaged heat flux \(uT-\nabla T\) in the vertical direction,
we obtain the mathematical definition of the Nusselt number

\(Nu=\int\limits_0^1 \langle (\mathbf{u} T-\nabla T)\cdot \mathbf{e_z}\rangle dz.\)

Since the convective fluid flow increases vertical heat transport beyond the purely conductive flux, our challange is to determine the relationship

\(Nu=Nu(Ra,Pr) \)

from the equations of motion.

In 1954 W.V.R Malkus [2] predicted the scaling law

\(Nu\sim Ra^{\frac{1}{3}}\)

by a *marginally stable boundary layer* argument,
based on the concept that
the boundary layer thickness \(\delta\) adjust itself so as to be,
as a convection layer, marginally stable.

The scaling

\(Nu\sim Pr^{\frac{1}{2}} Ra^{\frac{1}{2}}\)

has been postulated (Kraichnan (1962) [3] and Spiegel (1971) [4]) as an asymptotic regime in which the heat transfer and the strength of turbulence become independent of the kinematic viscosity and the thermal diffusivity.

## Infinite Prandtl-number limit

In many situations in nature, when the fluid is very viscous (e.g, Earth's mantel and engine oils), the Prandtl number is very big. This motivates the interest in studying the limiting case when the inertia of the fluid can be neglected:

\(\begin{cases} \frac{\partial }{\partial t}T+\mathbf{u}\cdot \nabla T=\Delta T\\ -\Delta \mathbf{u}+\nabla p=Ra T \mathbf{e_z}\\ \nabla\cdot \mathbf{u}=0\\ \end{cases}\)

In this case of infinite-Prandtl-number fluids the scaling

\(Nu\sim Ra^{\frac{1}{3}}\)

is believed to be valid in the asymptotically high-Rayleigh number regime (Grossmann & Lohse 2000 [5]).

The derivation of physically relevant upper bounds has a long history, that goes
back to the sixties (Howard (1963) [6] and Busse (1969) [7]). Several decades after,
the introduction of the *background field method* (Doering & Constantin [11]) has produced significant
results. This method consists of decomposing the temperature field into a background profile
and a perturbation term,

\(T(x,y,z,t)=\tau(z)+\theta(x,y,z,t)\) , where

\(\tau(0)=1\), \(\tau(1)=0\quad\) and \(\quad\theta(x,y,0,t)=0=\theta(x,y,1,t)\).

By the background field method (BFM) the problem of finding upper bounds for \(Nu\) is
reduced to the problem of constructing a background profile \(\tau\) that satisfies a *marginal* stability constraint.
Then the Dirichlet integral of a marginally stable
background temperature profile produces an upper bound
for the Nusselt number.

Here, the possibility to construct a "good" stable background profile is possible
thanks to the instantaneous *linear slaving* of the velocity
to the temperature field.

Within the BFM (via marginal stability)
many important results have been obtained.

The first rigorous upper bound, optimal up to logarithmic correction,
was proved by P. Constantin and C. Doering in [8]:

\(Nu\lesssim \big(\ln^{\frac{2}{3}}(Ra)\big) Ra^{\frac{1}{3}}.\)

Their proof relies on an \(L^{\infty}\) maximal regularity argument for the Stokes equation together
with bounds on the average Laplacian squared of the temperature.

The construction of a non-monotone background profile with a "log-layer" in the bulk
is the central idea in Doering, Otto & Westdickenberg [9] to derive

\(Nu\lesssim \big(\ln^{\frac{1}{3}}(Ra)\big) Ra^{\frac{1}{3}}.\)

By refinement of the argument in [9], F. Otto and C. Seis in [10] improved the last bound obtaining

\(Nu\lesssim \big(\ln^{\frac{1}{15}}(Ra)\big) Ra^{\frac{1}{3}}.\)

Despite the fact that the background field method has produced many significant upper bounds, it is not optimal. Indeed C. Nobili and F. Otto [in preparation] showed that the Nusselt number produced by the background field method

\( \widetilde{Nu}:=\inf\limits_{\substack{\tau:(0,1)\rightarrow \mathbb{R} \\ \tau(0)=1, \tau(1)=0}}\left\{\int_0^1\left(\frac{d\tau}{dz}\right)^2dz\quad \vert \;\;\tau \text{ marginally stable }\right\}\)

can be bounded from below by \(\big(\text{ln}^\frac{1}{15}\left(Ra\right)\big) Ra^\frac{1}{3}\).
In particular this implies that \(\widetilde{Nu} \sim \big(\text{ln}^\frac{1}{15}\left(Ra\right)\big) Ra^\frac{1}{3} \) and
no better bound can be produced by the background field method.

Nevertheless, the combination of the background stability with other methods
can be implemented to produce better bounds.

Indeed by the additional use of \(L^{\infty}\)
maximal regularity for Stokes equation, F. Otto and C. Seis [10] obtained the upper bound

\(Nu\lesssim \Big(\ln^{\frac{1}{3}}\big(\ln(Ra)\big)\Big) Ra^{\frac{1}{3}}.\)

which is, up to now, the best known.

## Finite Prandtl numbers

When \(Pr\) is finite, the equation for the temperature is coupled with the full Navier-Stokes equation for the velocity field, i.e

\(\frac{1}{Pr}\left(\frac{\partial \mathbf{u}}{\partial t}+\mathbf{u}\cdot \nabla \mathbf{u}\right)-\Delta \mathbf{u}+\nabla p=Ra T \mathbf{e_z},\)

\(\nabla\cdot \mathbf{u}=0.\)

As opposed to the case \(Pr=\infty\), here a major difficulty comes from the fact that the velocity and the temperature field are not
instantaneously slaved to each other.
Our current interest is in deriving rigorous upper bounds for \(Nu\)
that reproduce both physical scalings \(Nu\sim Ra^{\frac{1}{3}}\) and
\(Nu\sim Ra^{\frac{1}{2}}\) in some parameter regimes, up to logarithms.

In 1996 C. Doering and P. Constantin [11] applied the background field method
to the problem, finding

\(Nu\lesssim c Ra^{\frac{1}{2}}\)

for no-slip (or stress-free) boundary conditions.

Up to now, the best rigorous upper bound for large but finite Prandtl number
is proved by X. Wang in 2006 [12].
With a perturbation argument on the Stokes equation the author proves

\(Nu\lesssim \big(\ln^{\frac{2}{3}}(Ra)\big) Ra^{\frac{1}{3}} \;\; \mbox{ for }\;\; Pr\geq c_0 Ra\)

where \(c_0\) depends only on the aspect ratio of the domain.

Combining (logarithmically failing) maximal regularity estimates in \(L^{\infty}\) and \(L^1\) for the
nonstationary Stokes equation with force terms given by the bouyancy term and the nonlinear term respectively, A. Choffrut, C.Nobili and F.Otto
proved [in preparation] that

\( Nu \lesssim \begin{cases} \big( \text{ln}^\frac{1}{3}(Ra) \big) Ra^\frac{1}{3} &\;\; \mbox{ for }\;\; Pr \gtrsim \big( \text{ln}^\frac{1}{3}(Ra) \big) Ra^\frac{1}{3}\;, \\ \big( \text{ln}^\frac{1}{2}(Ra) \big) Pr^{-\frac{1}{2}} Ra^\frac{1}{2} &\;\; \mbox{ for }\;\; Pr \lesssim \big( \text{ln}^\frac{1}{3}(Ra) \big) Ra^\frac{1}{3}\;. \\ \end{cases} \)

The result on one hand improves the bound in [12] and on the other hand captures a regime where the inertial effects are so strong that the convection unlinearity is not negligible. In particular, the notion of solution for the Navier-Stokes equation is that of Leray.

## References

- Kadanoff, L. P.

Turbulent heat flow: structures and scaling.

Phys. Today 54 (8) (2001), 34-39. - Malkus, W.V.R.

The Heat Transport and Spectrum of Thermal Turbulence.

Proc. R: Soc. Lond. A 225, 196-2012 (1954) - Kraichnan, R.H.

Turbulent Thermal Convection at Arbitrary Prandtl Number.

Phys. Fluids 5 1374 (1962) - Spiegel, E.A.

Convection in Stars I. Basic Boussinesq Convection.

Ann. Rev Astron. Astrophys. 98, 323 (1971) - Grossmann, S. and Lohse, D.

Scaling in thermal convection: a unifying theory.

J. Fluid Mech. 407, 27 (2000). - Howard, L.

Heat transport by turbulent convection.

J. Fluid Mech. 17, 405 (1963) - Busse, F.H.

Fundamentals of Thermal Convection.

Peltier, W., eds.: Mantle Convections, Plate Tectonics and Global Dynamics , Gordon and Breach Publishing 23 (1989) - Doering, C.R and Constantin, P.

Infinite Prandtl number convection.

J. Stat. Phys. 94 (1-2) (1999), 159. - Charles R. Doering ; Felix Otto and Maria G. Westdickenberg.

Bounds on vertical heat transport for infinite-Prandtl-number Raleigh-Benard convection.

Journal of fluid mechanics, 560 (2006), 229 - Felix Otto and Christian Seis.

Rayleigh-Bénard convection : improved bounds on the Nusselt number.

Journal of mathematical physics, 52 (8) (2011) - Doering, C.R and Constantin, P.

Variational bounds on energy dissipation in incompressible flows: III. Convection.

Phys. Rev. E 53 5957 (1996) - Wang, X.

Bound on vertical heat transport at large Prandtl number.

Comm.Pure 61, 789 (2008)