Numerical calculations of two-dimensional large Prandtl number convection in a box

2013 ◽  
Vol 729 ◽  
pp. 584-602 ◽  
Author(s):  
J. A. Whitehead ◽  
A. Cotel ◽  
S. Hart ◽  
C. Lithgow-Bertelloni ◽  
W. Newsome
Keyword(s):  
Heat Flux ◽  
Rayleigh Number ◽  
Heat Source ◽  
Prandtl Number ◽  
Steady Flow ◽  
Maximum Velocity ◽  
Numerical Calculations ◽  
Two Dimensional ◽  
Long Time ◽  
Constant Viscosity

AbstractConvection from an isolated heat source in a chamber has been previously studied numerically, experimentally and analytically. These have not covered long time spans for wide ranges of Rayleigh number Ra and Prandtl number Pr. Numerical calculations of constant viscosity convection partially fill the gap in the ranges $\mathit{Ra}= 1{0}^{3} {{\unicode{x2013}}}1{0}^{6} $ and $\mathit{Pr}= 1, 10, 100, 1000$ and $\infty $. Calculations begin with cold fluid everywhere and localized hot temperature at the centre of the bottom of a square two-dimensional chamber. For $\mathit{Ra}\gt 20\hspace{0.167em} 000$, temperature increases above the hot bottom and forms a rising plume head. The head has small internal recirculation and minor outward conduction of heat during ascent. The head approaches the top, flattens, splits and the two remnants are swept to the sidewalls and diffused away. The maximum velocity and the top centre heat flux climb to maxima during head ascent and then adjust toward constant values. Two steady cells are separated by a vertical thermal conduit. This sequence is followed for every value of $Pr$ number, although lower Pr convection lags in time. For $\mathit{Ra}\lt 20\hspace{0.167em} 000$ there is no plume head, and no streamfunction and heat flux maxima with time. For sufficiently large Ra and all values of Pr, an oscillation develops at roughly $t= 0. 2$, with the two cells alternately strengthening and weakening. This changes to a steady flow with two unequal cells that at roughly $t= 0. 5$ develops a second oscillation.

Roll-diameter dependence in Rayleigh convection and its effect upon the heat flux

1972 ◽  
Vol 54 (2) ◽  
pp. 351-367 ◽  
Author(s):  
G. E. Willis ◽  
J. W. Deardorff ◽  
R. C. J. Somerville
Keyword(s):  
Heat Flux ◽  
Rayleigh Number ◽  
Prandtl Number ◽  
Silicone Oil ◽  
Numerical Calculations ◽  
Two Dimensional ◽  
Numerical Integrations ◽  
Rayleigh Convection

The average roll diameter in Rayleigh convection for 2000 < R < 31000, where R is the Rayleigh number, has been measured from photographs of three convecting fluids: air, water and a silicone oil with a Prandtl number σ of 450. For air the average dimensionless roll diameter was found to depend uniquely upon R and to increase especially rapidly in the range 2000 < R < 8000. The fluids of larger σ exhibited strong hysteresis but also had average roll diameters tending to increase with R. The increase in average roll diameter with R tended to decrease with σ. Through use of two-dimensional numerical integrations for the case of air it was found that the increase in average roll diameter with R provides an explanation for the usual discrepancy in heat flux observed between experiment and two-dimensional numerical calculations which prescribe a fixed wavelength.

Natural Convection in a Micro Size Horizontal Cylindrical Annulus With Periodic Volumetric Heat Flux

2009 ◽  
Author(s):  
Kamyar Mansour
Keyword(s):  
Heat Flux ◽  
Natural Convection ◽  
Rayleigh Number ◽  
Viscous Fluid ◽  
Two Dimensional ◽  
Dimensional Problem ◽  
Similarity Parameter ◽  
Symbolic Calculation ◽  
Cylindrical Annulus ◽  
Gap Ratio

We consider the two-dimensional problem of steady natural convection in a narrow (Micro size) Horizontal Cylindrical annulus filled with viscous fluid and periodic volumetric heat flux. The solution is expanded in powers of a single combined similarity parameter, which is the product of the Gap ratio to the power of four, and Rayleigh number and the series extended by means of symbolic calculation up to 16 terms. Analysis of these expansions allows the exact computation for arbitrarily accuracy up to 50000 figures. Although the range of the radius of convergence is almost zero but Pade approximation lead our result to be good even for much higher value of the similarity parameter.

scholarly journals Finite Prandtl number convection with nearly insulating boundaries

1981 ◽  
Vol 24 (3) ◽  
pp. 339-347
Author(s):  
N. Riahi
Keyword(s):  
Heat Flux ◽  
Prandtl Number ◽  
Maximum Amount ◽  
Two Dimensional ◽  
Thermal Conductivities

Arbitrary Prandtl number convection in a layer with nearly insulating boundaries is investigated. For B1/3 > 3.95P (P is the Prandtl number and B being the ratio between the thermal conductivities of the boundary and of the fluid) two-dimensional rolls are stable. The heat transported by the stable rolls reaches its peak at a critical P = B1/3/5.77 beyond which the heat flux decreases with increasing P. For B1/3 = 3.95P rolls become unstable to disturbances in the form of rolls oriented at a right angle to the original rolls. For B1/3 > 3.95P square pattern convection represents the preferred stable convection. The stable square pattern transports the maximum amount of heat at a critical P = B1/3/3.7.

The effect of distant side-walls on the evolution and stability of finite-amplitude Rayleigh-Bénard convection

1981 ◽  
Vol 378 (1775) ◽  
pp. 539-566 ◽  
Keyword(s):  
Rayleigh Number ◽  
Prandtl Number ◽  
Finite Amplitude ◽  
Two Dimensional ◽  
Bénard Convection ◽  
Benard Convection ◽  
Side Walls ◽  
Rayleigh Benard Convection ◽  
Rayleigh Numbers ◽  
Rayleigh Benard

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 ).

A study of Bénard convection with and without rotation

1969 ◽  
Vol 36 (2) ◽  
pp. 309-335 ◽  
Author(s):  
H. T. Rossby
Keyword(s):  
Boundary Layer ◽  
Heat Flux ◽  
Rayleigh Number ◽  
Prandtl Number ◽  
Linear Theory ◽  
Silicone Oil ◽  
Taylor Number ◽  
Subcritical Instability ◽  
Two Fluids ◽  
Wide Range

An experimental study of the response of a thin uniformly heated rotating layer of fluid is presented. It is shown that the stability of the fluid depends strongly upon the three parameters that described its state, namely the Rayleigh number, the Taylor number and the Prandtl number. For the two Prandtl numbers considered, 6·8 and 0·025 corresponding to water and mercury, linear theory is insufficient to fully describe their stability properties. For water, subcritical instability will occur for all Taylor numbers greater than 5 × 104, whereas mercury exhibits a subcritical instability only for finite Taylor numbers less than 105. At all other Taylor numbers there is good agreement between linear theory and experiment.The heat flux in these two fluids has been measured over a wide range of Rayleigh and Taylor numbers. Generally, much higher Nusselt numbers are found with water than with mercury. In water, at any Rayleigh number greater than 104, it is found that the Nusselt number will increase by about 10% as the Taylor number is increased from zero to a certain value, which depends on the Rayleigh number. It is suggested that this increase in the heat flux results from a perturbation of the velocity boundary layer with an ‘Ekman-layer-like’ profile in such a way that the scale of boundary layer is reduced. In mercury, on the other hand, the heat flux decreases monotonically with increasing Taylor number. Over a range of Rayleigh numbers (at large Taylor numbers) oscillatory convection is preferred although it is inefficient at transporting heat. Above a certain Rayleigh number, less than the critical value for steady convection according to linear theory, the heat flux increases more rapidly and the convection becomes increasingly irregular as is shown by the temperature fluctuations at a point in the fluid.Photographs of the convective flow in a silicone oil (Prandtl number = 100) at various rotation rates are shown. From these a rough estimate is obtained of the dominant horizontal convective scale as a function of the Rayleigh and Taylor numbers.

Natural convection flow in an open rectangular cavity with cold sidewalls and constant volumetric heat source

2011 ◽  
Vol 225 (5) ◽  
pp. 1191-1201 ◽  
Author(s):  
M Saleem ◽  
S Asghar ◽  
M A Hossain
Keyword(s):  
Natural Convection ◽  
Rayleigh Number ◽  
Finite Difference ◽  
Heat Source ◽  
Prandtl Number ◽  
Rectangular Cavity ◽  
Maximum Temperature ◽  
Convection Flow ◽  
Natural Convection Flow ◽  
Volumetric Heat Source

The transient two-dimensional natural convection flow of Newtonian fluid in an open rectangular cavity has been studied numerically. The flow is induced due to constant internal heat generation. The alternate direct implicit (ADI) finite difference, together with upwind finite-difference scheme and successive over relaxation method, are used to solve the non-dimensional equations numerically. Effects of Rayleigh number, Ra, Prandtl number, Pr, and cavity aspect ratio, A, on the flow patterns and isotherms as well as on the heat transfer rate are studied graphically. The maximum temperature induced due to the constant volumetric heat source is found with the increase in cavity width, and to decrease with the increase in Prandtl number and Rayleigh number. The numerical model employed here is found to be in good agreement with the previous work.

Development of hybrid method for coupled conduction-radiation heat transfer in two-dimensional irregular enclosure considering thermo-radiative effects and varying thermal conductivity

2019 ◽  
Vol 30 (4) ◽  
pp. 1815-1837
Author(s):  
Mehdi Zare ◽  
Sadegh Sadeghi
Keyword(s):  
Heat Transfer ◽  
Thermal Conductivity ◽  
Heat Flux ◽  
Heat Source ◽  
Radiation Heat Transfer ◽  
Isotropic Scattering ◽  
Two Dimensional ◽  
Radiation Heat ◽  
Content Type ◽  
Irregular Enclosure

Purpose This study aims to perform a comprehensive investigation to model the thermal characteristics of a coupled conduction-radiation heat transfer in a two-dimensional irregular enclosure including a triangular-shaped heat source. Design/methodology/approach For this purpose, a promising hybrid technique based on the concepts of blocked-off method, FVM and DOM is developed. The enclosure consists of several horizontal, vertical and oblique walls, and thermal conductivity within the enclosure varies directly with temperature and indirectly with position. To simplify the complex geometry, a promising mathematical model is introduced using blocked-off method. Emitting, absorbing and non-isotropic scattering gray are assumed as the main radiative characteristics of the steady medium. Findings DOM and FVM are, respectively, applied for solving radiative transfer equation (RTE) and the energy equation, which includes conduction, radiation and heat source terms. The temperature and heat flux distributions are calculated inside the enclosure. For validation, results are compared with previous data reported in the literature under the same conditions. Results and comparisons show that this approach is highly efficient and reliable for complex geometries with coupled conduction-radiation heat transfer. Finally, the effects of thermo-radiative parameters including surface emissivity, extinction coefficient, scattering albedo, asymmetry factor and conduction-radiation parameter on temperature and heat flux distributions are studied. Originality/value In this paper, a hybrid numerical method is used to analyze coupled conduction-radiation heat transfer in an irregular geometry. Varying thermal conductivity is included in this analysis. By applying the method, results obtained for temperature and heat flux distributions are presented and also validated by the data provided by several previous papers.

scholarly journals NATURAL CONVECTION IN HIGH ASPECT RATIO THREE-DIMENSIONAL ENCLOSURES WITH UNIFORM HEAT FLUX ON THE HEATED WALL

2004 ◽  
Vol 3 (2) ◽  
pp. 100
Author(s):  
T. Dias Jr. ◽  
L. F. Milanez
Keyword(s):  
Heat Flux ◽  
Nusselt Number ◽  
Natural Convection ◽  
Rayleigh Number ◽  
Aspect Ratio ◽  
Prandtl Number ◽  
High Aspect Ratio ◽  
Three Dimensional ◽  
Uniform Heat Flux ◽  
Uniform Heat

In this work, the laminar natural convection in high aspect ratio three-dimensional enclosures has been numerically studied. The enclosures studied here were heated with uniform heat flux on a vertical wall and cooled at constant temperature on the opposite wall. The remaining walls were considered adiabatic. Fluid properties were assumed constant except for the density change with temperature on the buoyancy term. The governing equations were solved using the finite volumes method and the dimensionless form of these equations has the Prandtl number and the modified Rayleigh number as parameters. The influences of the Rayleigh number and of the cavity aspect ratio on the Nusselt number, for a Prandtl number of 0.7, were analyzed. Results were obtained for values of the modified Rayleigh number up to 106 and for aspect ratios ranging from 1 to 20. The results were compared with two-dimensional results available in the literature and the variation of the average Nusselt number with the parameters studied were discussed.

Horizontal convection is non-turbulent

2002 ◽  
Vol 466 ◽  
pp. 205-214 ◽  
Author(s):  
F. PAPARELLA ◽  
W. R. YOUNG
Keyword(s):  
Thermal Diffusivity ◽  
Energy Dissipation ◽  
Rayleigh Number ◽  
Prandtl Number ◽  
Numerical Solutions ◽  
Inviscid Limit ◽  
Two Dimensional ◽  
Uniform Temperature ◽  
Horizontal Convection ◽  
Boussinesq Fluid

Consider the problem of horizontal convection: a Boussinesq fluid, forced by applying a non-uniform temperature at its top surface, with all other boundaries insulating. We prove that if the viscosity, ν, and thermal diffusivity, κ, are lowered to zero, with σ ≡ ν/κ fixed, then the energy dissipation per unit mass, κ, also vanishes in this limit. Numerical solutions of the two-dimensional case show that despite this anti-turbulence theorem, horizontal convection exhibits a transition to eddying flow, provided that the Rayleigh number is sufficiently high, or the Prandtl number σ sufficiently small. We speculate that horizontal convection is an example of a flow with a large number of active modes which is nonetheless not ‘truly turbulent’ because ε→0 in the inviscid limit.

Boundary-layer solutions of single-mode convection equations

1981 ◽  
Vol 102 ◽  
pp. 211-219 ◽  
Author(s):  
N. Riahi
Keyword(s):  
Boundary Layer ◽  
Heat Flux ◽  
Rayleigh Number ◽  
Prandtl Number ◽  
Thermal Convection ◽  
Single Mode ◽  
Layer Method ◽  
Boundary Layer Method ◽  
Inertial Terms ◽  
Large Rayleigh Number

Nonlinear thermal convection between two stress-free horizontal boundaries is studied using the modal equations for cellular convection. Assuming a large Rayleigh number R the boundary-layer method is used for different ranges of the Prandtl number σ. The heat flux F is determined for the values of the horizontal wavenumber a which maximizes F. For a large Prandtl number, σ [Gt ] R⅙(log R)−1, inertial terms are insignificant, a is either of order one (for $\sigma \geqslant R^{\frac{2}{3}}$) or proportional to $R^{\frac{1}{3}}\sigma^{-\frac{1}{2}}$ (for $\sigma \ll R^{\frac{2}{3}}$) and F is proportional to $R^{\frac{1}{3}}$. For a moderate Prandtl number, \[ (R^{-1}\log R)^{\frac{1}{9}} \ll \sigma \ll R^{\frac{1}{6}}(\log R)^{-1}, \] inertial terms first become significant in an inertial layer adjacent to the viscous buoyancy-dominated interior, and a and F are proportional to R¼ and \[ R^{\frac{3}{10}}\sigma^{\frac{1}{5}} (\log\sigma R^{\frac{1}{4}})^{\frac{1}{10}}, \] respectively. For a small Prandtl number, $R^{-1} \ll \sigma \ll (R^{-1} \log R)^{\frac{1}{9}}$, inertial terms are significant both in the interior and the boundary layers, and a and F are proportional to ($R \sigma)^{\frac{9}{32}} (\log R\sigma)^{-\frac{1}{32}}$ and ($R \sigma)^{\frac{5}{16}} (\log R \sigma)^{\frac{3}{16}}$, respectively.

Sign in / Sign up

Export Citation Format

Share Document