Influence of viscous dissipation on Bénard convection

1974 ◽  
Vol 64 (2) ◽  
pp. 369-374 ◽  
Author(s):  
D. L. Turcotte ◽  
A. T. Hsui ◽  
K. E. Torrance ◽  
G. Schubert

The approximations implicit in Bénard convection have been modified to include viscous dissipation. It is shown that both the influence of an adiabatic temperature gradient and of viscous dissipation are governed by the same dimensionless parameter Di = αgh/cp. Numerical calculations of finite amplitude convection are given for finite values of Di. It is found that increasing Di decreases flow velocities and finally stabilizes the flow.

1971 ◽  
Vol 42 ◽  
pp. 130-135 ◽  
Author(s):  
K. H. Böhm ◽  
J. Cassinelli

Outer convection zones of white dwarfs in the range 5800 K ≤ Teff ≤ 30000 K have been studied assuming that they have the same chemical composition as determined by Weidemann (1960) for van Maanen 2. Convection is important in all these stars. In white dwarfs Teff < 8000 K the adiabatic temperature gradient is strongly influenced by the pressure ionization of H, HeI and HeII which occurs within the convection zone. Partial degeneracy is also important.Convective velocities are very small for cool white dwarfs but they reach considerable values for hotter objects. For a white dwarf of Teff = 30000 K a velocity of 6.05 km/sec and an acoustic flux (generated by the turbulent convection) of 1.5 × 1011 erg cm−2 sec−1 is reached. The formation of white dwarf coronae is briefly discussed.


2018 ◽  
Vol 841 ◽  
pp. 825-850 ◽  
Author(s):  
Chong Shen Ng ◽  
Andrew Ooi ◽  
Detlef Lohse ◽  
Daniel Chung

Previous numerical studies on homogeneous Rayleigh–Bénard convection, which is Rayleigh–Bénard convection (RBC) without walls, and therefore without boundary layers, have revealed a scaling regime that is consistent with theoretical predictions of bulk-dominated thermal convection. In this so-called asymptotic regime, previous studies have predicted that the Nusselt number ($\mathit{Nu}$) and the Reynolds number ($\mathit{Re}$) vary with the Rayleigh number ($\mathit{Ra}$) according to $\mathit{Nu}\sim \mathit{Ra}^{1/2}$ and $\mathit{Re}\sim \mathit{Ra}^{1/2}$ at small Prandtl numbers ($\mathit{Pr}$). In this study, we consider a flow that is similar to RBC but with the direction of temperature gradient perpendicular to gravity instead of parallel to it; we refer to this configuration as vertical natural convection (VC). Since the direction of the temperature gradient is different in VC, there is no exact relation for the average kinetic dissipation rate, which makes it necessary to explore alternative definitions for $\mathit{Nu}$, $\mathit{Re}$ and $\mathit{Ra}$ and to find physical arguments for closure, rather than making use of the exact relation between $\mathit{Nu}$ and the dissipation rates as in RBC. Once we remove the walls from VC to obtain the homogeneous set-up, we find that the aforementioned $1/2$-power-law scaling is present, similar to the case of homogeneous RBC. When focusing on the bulk, we find that the Nusselt and Reynolds numbers in the bulk of VC too exhibit the $1/2$-power-law scaling. These results suggest that the $1/2$-power-law scaling may even be found at lower Rayleigh numbers if the appropriate quantities in the turbulent bulk flow are employed for the definitions of $\mathit{Ra}$, $\mathit{Re}$ and $\mathit{Nu}$. From a stability perspective, at low- to moderate-$\mathit{Ra}$, we find that the time evolution of the Nusselt number for homogenous vertical natural convection is unsteady, which is consistent with the nature of the elevator modes reported in previous studies on homogeneous RBC.


2015 ◽  
Vol 773 ◽  
pp. 395-417 ◽  
Author(s):  
K. Petschel ◽  
S. Stellmach ◽  
M. Wilczek ◽  
J. Lülff ◽  
U. Hansen

The kinetic energy balance in Rayleigh–Bénard convection is investigated by means of direct numerical simulations for the Prandtl number range $0.01\leqslant \mathit{Pr}\leqslant 150$ and for fixed Rayleigh number $\mathit{Ra}=5\times 10^{6}$. The kinetic energy balance is divided into a dissipation, a production and a flux term. We discuss the profiles of all the terms and find that the different contributions to the energy balance can be spatially separated into regions where kinetic energy is produced and where kinetic energy is dissipated. By analysing the Prandtl number dependence of the kinetic energy balance, we show that the height dependence of the mean viscous dissipation is closely related to the flux of kinetic energy. We show that the flux of kinetic energy can be divided into four additive contributions, each representing a different elementary physical process (advection, buoyancy, normal viscous stresses and viscous shear stresses). The behaviour of these individual flux contributions is found to be surprisingly rich and exhibits a pronounced Prandtl number dependence. Different flux contributions dominate the kinetic energy transport at different depths, such that a comprehensive discussion requires a decomposition of the domain into a considerable number of sublayers. On a less detailed level, our results reveal that advective kinetic energy fluxes play a key role in balancing the near-wall dissipation at low Prandtl number, whereas normal viscous stresses are particularly important at high Prandtl number. Finally, our work reveals that classical velocity boundary layers are deeply connected to the kinetic energy transport, but fail to correctly represent regions of enhanced viscous dissipation.


A recent study by Cross et al . (1980) has described a class of finite-amplitude phase-winding solutions of the problem of two-dimensional Rayleigh-Bénard convection in a shallow fluid layer of aspect ratio 2 L (≫ 1) confined laterally by rigid side-walls. These solutions arise at Rayleigh numbers R = R 0 + O ( L -1 ) where R 0 is the critical Rayleigh number for the corresponding infinite layer. Nonlinear solutions of constant phase exist for Rayleigh numbers R = R 0 + O ( L -2 ) but of these only the two that bifurcate at the lowest value of R are stable to two-dimensional linearized disturbances in this range (Daniels 1978). In the present paper one set of the class of phase-winding solutions is found to be stable to two-dimensional disturbances. For certain values of the Prandtl number of the fluid and for stress-free horizontal boundaries the results predict that to preserve stability there must be a continual readjustment of the roll pattern as the Rayleigh number is raised, with a corresponding increase in wavelength proportional to R - R 0 . These solutions also exhibit hysteresis as the Rayleigh number is raised and lowered. For other values of the Prandtl number the number of rolls remains unchanged as the Rayleigh number is raised, and the wavelength remains close to its critical value. It is proposed that the complete evolution of the flow pattern from a static state must take place on a number of different time scales of which t = O(( R - R 0 ) -1 ) and t = O(( R - R 0 ) -2 ) are the most significant. When t = O(( R - R 0 ) -1 ) the amplitude of convection rises from zero to its steady-state value, but the final lateral positioning of the rolls is only completed on the much longer time scale t = O(( R - R 0 ) -2 ).


1980 ◽  
Vol 96 (3) ◽  
pp. 515-583 ◽  
Author(s):  
Gary T. Jarvis ◽  
Dan P. Mckenzie

An approximate set of equations is derived for a compressible liquid of infinite Prandtl number. These are referred to as the anelastic-liquid equations. The approximation requires the product of absolute temperature and volume coefficient of thermal expansion to be small compared to one. A single parameter defined as the ratio of the depth of the convecting layer,d, to the temperature scale height of the liquid,HT, governs the importance of the non-Boussinesq effects of compressibility, viscous dissipation, variable adiabatic temperature gradients and non-hydrostatic pressure gradients. Whend/HT[Lt ] 1 the Boussinesq equations result, but whend/HTisO(1) the non-Boussinesq terms become important. Using a time-dependent numerical model, the anelastic-liquid equations are solved in two dimensions and a systematic investigation of compressible convection is presented in whichd/HTis varied from 0·1 to 1·5. Both marginal stability and finite-amplitude convection are studied. Ford/HT[les ] 1·0 the effect of density variations is primarily geometric; descending parcels of liquid contract and ascending parcels expand, resulting in an increase in vorticity with depth. Whend/HT> 1·0 the density stratification significantly stabilizes the lower regions of the marginal state solutions. At all values ofd/HT[ges ] 0·25, an adiabatic temperature gradient proportional to temperature has a noticeable stabilizing effect on the lower regions. Ford/HT[ges ] 0·5, marginal solutions are completely stabilized at the bottom of the layer and penetrative convection occurs for a finite range of supercritical Rayleigh numbers. In the finite-amplitude solutions adiabatic heating and cooling produces an isentropic central region. Viscous dissipation acts to redistribute buoyancy sources and intense frictional heating influences flow solutions locally in a time-dependent manner. The ratio of the total viscous heating in the convecting system, ϕ, to the heat flux across the upper surface,Fu, has an upper limit equal tod/HT. This limit is achieved at high Rayleigh numbers, when heating is entirely from below, and, for sufficiently large values ofd/HT, Φ/Fuis greater than 1·00.


Sign in / Sign up

Export Citation Format

Share Document