skip to main content


Title: The Nusselt numbers of horizontal convection
In the problem of horizontal convection a non-uniform buoyancy, $b_{s}(x,y)$ , is imposed on the top surface of a container and all other surfaces are insulating. Horizontal convection produces a net horizontal flux of buoyancy, $\boldsymbol{J}$ , defined by vertically and temporally averaging the interior horizontal flux of buoyancy. We show that $\overline{\boldsymbol{J}\boldsymbol{\cdot }\unicode[STIX]{x1D735}b_{s}}=-\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ ; the overbar denotes a space–time average over the top surface, angle brackets denote a volume–time average and $\unicode[STIX]{x1D705}$ is the molecular diffusivity of buoyancy  $b$ . This connection between $\boldsymbol{J}$ and $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ justifies the definition of the horizontal-convective Nusselt number, $Nu$ , as the ratio of $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ to the corresponding quantity produced by molecular diffusion alone. We discuss the advantages of this definition of $Nu$ over other definitions of horizontal-convective Nusselt number. We investigate transient effects and show that $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ equilibrates more rapidly than other global averages, such as the averaged kinetic energy and bottom buoyancy. We show that $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ is the volume-averaged rate of Boussinesq entropy production within the enclosure. In statistical steady state, the interior entropy production is balanced by a flux through the top surface. This leads to an equivalent ‘surface Nusselt number’, defined as the surface average of vertical buoyancy flux through the top surface times the imposed surface buoyancy $b_{s}(x,y)$ . In experimental situations it is easier to evaluate the surface entropy flux, rather than the volume integral of $|\unicode[STIX]{x1D735}b|^{2}$ demanded by $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ .  more » « less
Award ID(s):
1657041
NSF-PAR ID:
10169194
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
894
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The central open question about Rayleigh–Bénard convection – buoyancy-driven flow in a fluid layer heated from below and cooled from above – is how vertical heat flux depends on the imposed temperature gradient in the strongly nonlinear regime where the flows are typically turbulent. The quantitative challenge is to determine how the Nusselt number $Nu$ depends on the Rayleigh number $Ra$ in the $Ra\to \infty$ limit for fluids of fixed finite Prandtl number $Pr$ in fixed spatial domains. Laboratory experiments, numerical simulations and analysis of Rayleigh's mathematical model have yet to rule out either of the proposed ‘classical’ $Nu \sim Ra^{1/3}$ or ‘ultimate’ $Nu \sim Ra^{1/2}$ asymptotic scaling theories. Among the many solutions of the equations of motion at high $Ra$ are steady convection rolls that are dynamically unstable but share features of the turbulent attractor. We have computed these steady solutions for $Ra$ up to $10^{14}$ with $Pr=1$ and various horizontal periods. By choosing the horizontal period of these rolls at each $Ra$ to maximize $Nu$ , we find that steady convection rolls achieve classical asymptotic scaling. Moreover, they transport more heat than turbulent convection in experiments or simulations at comparable parameters. If heat transport in turbulent convection continues to be dominated by heat transport in steady rolls as $Ra\to \infty$ , it cannot achieve the ultimate scaling. 
    more » « less
  2. A numerical investigation of an asymptotically reduced model for quasigeostrophic Rayleigh-Bénard convection is conducted in which the depth-averaged flows are numerically suppressed by modifying the governing equations. At the largest accessible values of the Rayleigh number Ra, the Reynolds number and Nusselt number show evidence of approaching the diffusion-free scalings of Re ∼ RaE/Pr and Nu ∼ Pr−1/2Ra3/2E2, respectively, where E is the Ekman number and Pr is the Prandtl number. For large Ra, the presence of depth-invariant flows, such as large-scale vortices, yield heat and momentum transport scalings that exceed those of the diffusion-free scaling laws. The Taylor microscale does not vary significantly with increasing Ra, whereas the integral length scale grows weakly. The computed length scales remain O(1) with respect to the linearly unstable critical wave number; we therefore conclude that these scales remain viscously controlled. We do not find a point-wise Coriolis-inertia-Archimedean (CIA) force balance in the turbulent regime; interior dynamics are instead dominated by horizontal advection (inertia), vortex stretching (Coriolis) and the vertical pressure gradient. A secondary, subdominant balance between the Archimedean buoyancy force and the viscous force occurs in the interior and the ratio of the root mean square (rms) of these two forces is found to approach unity with increasing Ra. This secondary balance is attributed to the turbulent fluid interior acting as the dominant control on the heat transport. These findings indicate that a pointwise CIA balance does not occur in the high Rayleigh number regime of quasigeostrophic convection in the plane layer geometry. Instead, simulations are characterized by what may be termed a nonlocal CIA balance in which the buoyancy force is dominant within the thermal boundary layers and is spatially separated from the interior Coriolis and inertial forces. 
    more » « less
  3. This work employs single-mode equations to study convection and double-diffusive convection in a porous medium where the Darcy law provides large-scale damping. We first consider thermal convection with salinity as a passive scalar. The single-mode solutions resembling steady convection rolls reproduce the qualitative behavior of root-mean-square and mean temperature profiles of time-dependent states at high Rayleigh numbers from direct numerical simulations (DNS). We also show that the single-mode solutions are consistent with the heat-exchanger model that describes well the mean temperature gradient in the interior. The Nusselt number predicted from the single-mode solutions exhibits a scaling law with Rayleigh number close to that followed by exact 2D steady convection rolls, although large aspect ratio DNS results indicate a faster increase. However, the single-mode solutions at a high wavenumber predict Nusselt numbers close to the DNS results in narrow domains. We also employ the single-mode equations to analyze the influence of active salinity, introducing a salinity contribution to the buoyancy, but with a smaller diffusivity than the temperature. The single-mode solutions are able to capture the stabilizing effect of an imposed salinity gradient and describe the standing and traveling wave behaviors observed in DNS. The Sherwood numbers obtained from single-mode solutions show a scaling law with the Lewis number that is close to the DNS computations with passive or active salinity. This work demonstrates that single-mode solutions can be successfully applied to this system whenever periodic or no-flux boundary conditions apply in the horizontal. 
    more » « less
  4. null (Ed.)
    This study explores thermal convection in suspensions of neutrally buoyant, non-colloidal suspensions confined between horizontal plates. A constitutive diffusion equation is used to model the dynamics of the particles suspended in a viscous fluid and it is coupled with the flow equations. We employ a simple model that was proposed by Metzger, Rahli & Yin ( J. Fluid Mech. , vol. 724, 2013, pp. 527–552) for the effective thermal diffusivity of suspensions. This model considers the effect of shear-induced diffusion and gives the thermal diffusivity increasing linearly with the thermal Péclet number ( Pe ) and the particle volume fraction ( ϕ ). Both linear stability analysis and numerical simulation based on the mathematical models are performed for various bulk particle volume fractions $({\phi _b})$ ranging from 0 to 0.3. The critical Rayleigh number $(R{a_c})$ grows gradually by increasing ${\phi _b}$ from the critical value $(R{a_c} = 1708)$ for a pure Newtonian fluid, while the critical wavenumber $({k_c})$ remains constant at 3.12. The transition from the conduction state of suspensions is subcritical, whereas it is supercritical for the convection in a pure Newtonian fluid $({\phi _b} = 0)$ . The heat transfer in moderately dense suspensions $({\phi _b} = 0.2\text{--}0.3)$ is significantly enhanced by convection rolls for small Rayleigh number ( Ra ) close to $R{a_c}$ . We also found a power-law increase of the Nusselt number ( Nu ) with Ra , namely, $Nu\sim R{a^b}$ for relatively large values of Ra where the scaling exponent b decreases with ${\phi _b}$ . Finally, it turns out that the shear-induced migration of particles can modify the heat transfer. 
    more » « less
  5. Turbulence and mixing in a near-bottom convectively driven flow are examined by numerical simulations of a model problem: a statically unstable disturbance at a slope with inclination $\unicode[STIX]{x1D6FD}$ in a stable background with buoyancy frequency $N$ . The influence of slope angle and initial disturbance amplitude are quantified in a parametric study. The flow evolution involves energy exchange between four energy reservoirs, namely the mean and turbulent components of kinetic energy (KE) and available potential energy (APE). In contrast to the zero-slope case where the mean flow is negligible, the presence of a slope leads to a current that oscillates with $\unicode[STIX]{x1D714}=N\sin \unicode[STIX]{x1D6FD}$ and qualitatively changes the subsequent evolution of the initial density disturbance. The frequency, $N\sin \unicode[STIX]{x1D6FD}$ , and the initial speed of the current are predicted using linear theory. The energy transfer in the sloping cases is dominated by an oscillatory exchange between mean APE and mean KE with a transfer to turbulence at specific phases. In all simulated cases, the positive buoyancy flux during episodes of convective instability at the zero-velocity phase is the dominant contributor to turbulent kinetic energy (TKE) although the shear production becomes increasingly important with increasing  $\unicode[STIX]{x1D6FD}$ . Energy that initially resides wholly in mean available potential energy is lost through conversion to turbulence and the subsequent dissipation of TKE and turbulent available potential energy. A key result is that, in contrast to the explosive loss of energy during the initial convective instability in the non-sloping case, the sloping cases exhibit a more gradual energy loss that is sustained over a long time interval. The slope-parallel oscillation introduces a new flow time scale $T=2\unicode[STIX]{x03C0}/(N\sin \unicode[STIX]{x1D6FD})$ and, consequently, the fraction of initial APE that is converted to turbulence during convective instability progressively decreases with increasing $\unicode[STIX]{x1D6FD}$ . For moderate slopes with $\unicode[STIX]{x1D6FD}<10^{\circ }$ , most of the net energy loss takes place during an initial, short ( $Nt\approx 20$ ) interval with periodic convective overturns. For steeper slopes, most of the energy loss takes place during a later, long ( $Nt>100$ ) interval when both shear and convective instability occur, and the energy loss rate is approximately constant. The mixing efficiency during the initial period dominated by convectively driven turbulence is found to be substantially higher (exceeds 0.5) than the widely used value of 0.2. The mixing efficiency at long time in the present problem of a convective overturn at a boundary varies between 0.24 and 0.3. 
    more » « less