Average and extremal properties of heat transfer and shear stress on a wall surface in Rayleigh–Bénard convection

In this study surface-averaged and extremal properties of heat transfer and shear stress on the upper wall surface of Rayleigh–Bénard convection are numerically examined. The Prandtl number was raised up to 10³, and the Rayleigh number was changed between 10⁴ and 10⁷. As a result, average Nusselt number Nu and shear rate τ/Pr depends on Pr, Ra, and the entire numerical results are distributed between two correlation equations corresponding to small and large Pr. The small and large Pr equations are closely related to steady and unsteady flow regimes, respectively. Nevertheless, a single relation τ/Pr ~ Nu ^3.0 exists to explain the entire results. Similarly the change of local maximal properties Nu _max and τ _max/Pr depends on Pr, Ra, and these values are also distributed between two correlation equations corresponding to small and large Pr cases. Despite such complicated dependence we can obtain a correlation equation as a form of τ _max/Pr ~ Nu _max^2.6, which has not been obtained theoretically.     

Onset of Buoyancy-Driven Convection in Porous Media Saturated with Cold Water Cooled from Above

When porous media saturated with initially stagnant cold water around the density maximum temperature are cooled from above, convection may be induced in an unstable lower layer. In this study, the onset of buoyancy-driven convection during time-dependent cooling is investigated using the propagation theory, which transforms disturbance equations similarly, and also considering the density inversion effect. The critical Darcy–Rayleigh number Ra _D,c is found as a function of the dimensionless density maximum temperature θ _max. For Ra _D > Ra _D,c the dimensionless critical time τ _c to mark the onset of instability is presented as a function of Ra _D and θ _max. These critical conditions are compared with previous theoretical results.     

Nonsimilarity solutions for non-Darcy mixed convection from horizontal surfaces in a porous medium

Non-Darcy mixed convection in a porous medium from horizontal surfaces with variable surface heat flux of the power-law distribution is analyzed. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ζ _f =Ra^* _x /Pe² _x is found to characterize the effect of buoyancy forces on the forced convection with K ^′ U _∞/ν characterizing the effect of inertia resistance. The second region covers the natural convection dominated regime where the dimensionless parameter ζ _n =Pe _x /Ra^*1/2 _x is found to characterize the effect of the forced flow on the natural convection, with (K ^′ U _∞/ν)Ra^*1/2 _x /Pe _x characterizing the effect of inertia resistance. To obtain the solution that covers the entire mixed convection regime the solution of the first regime is carried out for ζ _f =0, the pure forced convection limit, to ζ _f =1 and the solution of the second is carried out for ζ _n =0, the pure natural convection limit, to ζ _n =1. The two solutions meet and match at ζ _f =ζ _n =1, and R ^* _h =G ^* _h . Also a non-Darcy model was used to analyze mixed convection in a porous medium from horizontal surfaces with variable wall temperature of the power-law form. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ξ _f =Ra _x /Pe _x ^3/2 is found to measure the buoyancy effects on mixed convection with Da _x Pe _x /ɛ as the wall effects. The second region covers the natural convection dominated region where ξ _n =Pe _x /Ra _x ^2/3 is found to measure the force effects on mixed convection with Da _x Ra _x ^2/3/ɛ as the wall effects. Numerical results for different inertia, wall, variable surface heat flux and variable wall temperature exponents are presented. Received on 8 July 1996     

Nusselt-Rayleigh correlations for free convection in 2D air-filled parallelogrammic enclosures with isothermal active walls

In the present study Nu-Ra-α correlations are proposed to calculate the steady-state natural convection heat transfer taking place in 2D air-filled cavities of parallelogrammic section. The thermal conditions and the dimensions of the enclosures permit to cover a large range of Rayleigh numbers, 1.7 × 10³  ≤ Ra ≤ 3.0 × 10⁹, suitable for diverse engineering applications. The two active walls of the cavities are kept vertical and isothermal at hot and cold temperatures T _h and T _c respectively. Separated by a horizontal distance H, they have the same height H and are connected by a closed adiabatic channel whose upper and lower walls can be inclined at an angle α with respect to the horizontal, varying between −60° to +60°. That gives rise to a conducting or insulating cavity, in the convective sense of the term (diode cavity). A computational model based on the finite volume method is used to solve the governing equations. The large number of treated configurations led to propose Nu-Ra-α correlations for large ranges of Ra and α which can be applied to many engineering areas. The results of this numerical study have been successfully compared with calculated and measured available data.     

Onset of Buoyancy-driven Instability in Porous Media Melted from Below

When a nonhomogeneous solid is melting from below, convection may be induced in a thermally–unstable melt layer. In this study, the onset of buoyancy-driven convection during time-dependent melting is investigated by using similarly transformed disturbance equations. The critical Darcy–Rayleigh numbers based on the melt-layer thickness, Ra _H,c, are found numerically for various conditions. For small superheats, the present predictions show that Ra _H,c is located between 27.1 and 4π ² and it approaches the well-known results of the original Horton–Rogers–Lapwood problem. However, for high superheats, it is dependent on the phase change rate λ and the relation of Ra _H,c λ = 25.89 is shown.     

Recoverable strains and retardation times of monodisperse suspensions of silicon dioxide spheres in poly(dimethylsiloxane)

Steady-state viscosities η, steady-state recoverable strains γ _rs and characteristic retardation time τ _1/2 were measured for suspensions of monodisperse silicon dioxide (SiO₂) spheres in poly(dimethylsiloxane) (PDMS) with various volume fractions Φ of the suspended spheres at various creep stresses σ ₀. Two different regions are found in plots of η/η _m vs γ _rs, where η/η _m denotes the relative viscosity of the suspensions. In one region, η/η _m is proportional to γ _rs, while γ _rs is independent of η/η _m in the other region. In both regions, τ _1/2 is the functions of the shear strain rate in the steady-state of creep test independently of Φ. The origin of the elasticity is related to the ‘maximally distorted’ cages recovered owing to the repulsive interaction between the SiO₂ spheres and recovery of the cages in the shear-induced clusters of the suspended spheres.     

The laminar hole pressure for Newtonian fluids

For flows with wall turbulence the hole pressure, P _H , was shown empirically by Franklin and Wallace (J Fluid Mech, 42, 33–48, 1970) to depend solely on R ₊, the Reynolds number constructed from the friction velocity and the hole diameter b. Here this dependence is extended to the laminar regime by numerical simulation of a Newtonian fluid flowing in a plane channel (gap H) with a deep tap hole on one wall. Calculated hole pressures are in good agreement with experimental values, and for two hole sizes are well represented by: (P _H −P _HS )/τ _w = √(k ² + c ² R ₊²)−k, where the Stokes hole pressure P _HS /τ _w = s (b/H)³, k, c, s are fitted constants, and τ _w is the wall shear stress.     

Stability of double-diffusive convection induced by selective absorption of radiation in a fluid layer

Linear and nonlinear stability analyses were performed on a fluid layer with a concentration-based internal heat source. Clear bimodal behaviour in the neutral curve (with stationary and oscillatory modes) is observed in the region of the onset of oscillatory convection, which is a previously unobserved phenomenon in radiation-induced convection. The numerical results for the linear instability analysis suggest a critical value γ _c of γ, a measure for the strength of the internal heat source, for which oscillatory convection is inhibited when γ > γ _c . Linear instability analyses on the effect of varying the ratio of the salt concentrations at the upper and lower boundaries conclude that the ratio has a significant effect on the stability boundary. A nonlinear analysis using an energy approach confirms that the linear theory describes the stability boundary most accurately when γ is such that the linear theory predicts the onset of mostly stationary convection. Nevertheless, the agreement between the linear and nonlinear stability thresholds deteriorates for larger values of the solute Rayleigh number for any value of γ.     

The laminar plume above a line heat source in a transverse magnetic field

The dynamics of a buoyant plume rising above a horizontal line heat source in a transverse, horizontal magnetic field is investigated. Similarity is shown to occur when the magnetic field strength varies as the −2/5 power of vertical distance from the source. The plume depends on two parameters — the Prandtl number (Pr) and the Lykoudis number (Z _L). Families of exact closed form solutions are derived for Pr=5/9 and Pr≥2. A family of numerical integrations for Pr=0.01 (typical of liquid metals) is also reported. The magnetic field is shown to affect the profiles of velocity and temperature by altering the similarity functions, the coefficients, and the value of the independent similarity variable corresponding to a fixed physical position. An approximate closed form solution valid for low Pr and high Z _L is presented. Possible experimental tests of the theory are proposed. Research sponsored by the U.S. Energy Research and Development Administration under interagency agreement with Union Carbide Corporation.     

Further Study on Pure Shear

It is known that the Cauchy stress tensor T is a pure shear when trT = 0. An elementary derivation is given for a coordinate system such that, when referred to this coordinate system, the diagonal elements of T vanish while the off-diagonal elements τ ₁, τ ₂, τ ₃, are the pure shears. The structure of τ _i (i = 1, 2, 3) depends on one non-dimensional parameter q = 54(detT)² / [tr(T ²)]³, 0 ≤ q ≤ 1. When q = 0, one of the three τ _i vanishes. A coordinate system can be chosen such that the remaining two have the same magnitude or one of the remaining two also vanishes. When q = 1, all three τ _i have the same magnitude. However, there is a one-parameter family of coordinate systems that gives the same three τ _i . For q ≠ 0 or 1, none of the three τ _i vanishes and the three τ _i in general have different magnitudes. Nevertheless, a coordinate system can be chosen such that two of the three τ _i have the same magnitude. ^★Professor Emeritus of University of Illinois at Chicago and Consulting Professor of Stanford University.     

Control of Flow and Heat Transfer in a Porous Enclosure due to an Adiabatic Thin fin on the Hot Wall

Natural convection flow in a differentially heated square enclosure filled with porous matrix with a solid adiabatic thin fin attached at the hot left wall is studied numerically. The Brinkman–Forchheimer-extended Darcy model is used to solve the momentum equations, in the porous medium. The numerical investigation is done through streamlines, isotherms, and heat transfer rates. A parametric study is carried out using the following parameters: Darcy number (Da) from 10⁻⁴ to 10⁻², dimensionless thin fin lengths (L _p) 0.3, 0.5, and 0.7, dimensionless positions (S _p) 0.25, 0.5, and 0.75 with Prandtl numbers (Pr) 0.7 and 100 for Ra = 10⁶. For Da = 10⁻³ and Pr = 0.7, it is observed that there is a counter clock-wise secondary flow formation around the tip of the fin for S _p = 0.5 for all lengths of L _p. Moreover when Da = 10⁻² the secondary circulation behavior has been observed for S _p = 0.25 and 0.75 and there is another circulation between the top wall and the fin that is separated from the primary circulation. However, these secondary circulations features are not observed for Pr = 100. It is also found that the average Nusselt number decreases as the length of the fin increases for all locations. However, the rate of decrease of average Nusselt number becomes slower as the location of fin moves from the bottom wall to the top wall. The overall heat transfer rate can be controlled with a suitable selection of the fin location and length.     

Prediction of natural convection flow using network model and numerical simulations inside enclosure with distributed solid blocks

Steady state natural convection of a fluid with Pr ≈ 1 within a square enclosure containing uniformly distributed, conducting square solid blocks is investigated. The side walls are subjected to differential heating, while the top and bottom ones are kept adiabatic. The natural convection flow is predicted employing the nondimensional volumetric flow rate (Q_max^* Q_{\max }^{*} ) by using a network model and also using numerical simulations. For identical solid and fluid thermal conductivities (i.e. k _s  = k _f ), a parametric study of the effect of number of blocks (N ²), gap size (δ) and enclosure Rayleigh number (Ra) on Q_max^* Q_{\max }^{*} is performed using the two approaches. Network model predictions are observed to agree well with that from the simulations until Raδ³ ~ 12. Considering the enclosure with blocks as a porous medium, for a fixed enclosure Ra number, increasing the number of blocks for a fixed volumetric porosity leads to a decrease in enclosure permeability, which in turn reduces the flow rate. When the number of blocks is fixed, and for a given Ra number, the flow rate increases as the porosity increases by widening the gap between the blocks.     

Thermophoresis particle deposition in a non-Darcy porous medium under the influence of Soret, Dufour effects

Thermophoresis particle deposition in free convection on a vertical plate embedded in a fluid saturated non-Darcy porous medium is studied using similarity solution technique. The effect of Soret and Dufour parameters on concentration distribution, wall thermophoretic deposition velocity, heat transfer and mass transfer is discussed in detail for different values of dispersion parameters (Ra _γ, Ra _ξ) inertial parameter F and Lewis number Le. The result indicates that the Soret effect is more influential in increasing the concentration distribution in both aiding as well as opposing buoyancies. Also, the non-dimensional heat transfer coefficient and non-dimensional mass transfer coefficient changes according to different values of thermophoretic coefficient k.     

Scaling law in saddle-node bifurcations for one-dimensional maps: a complex variable approach

The study of transient dynamical phenomena near bifurcation thresholds has attracted the interest of many researchers due to the relevance of bifurcations in different physical or biological systems. In the context of saddle-node bifurcations, where two or more fixed points collide annihilating each other, it is known that the dynamics can suffer the so-called delayed transition. This phenomenon emerges when the system spends a lot of time before reaching the remaining stable equilibrium, found after the bifurcation, because of the presence of a saddle-remnant in phase space. Some works have analytically tackled this phenomenon, especially in time-continuous dynamical systems, showing that the time delay, τ, scales according to an inverse square-root power law, τ∼(μ−μ _c )^−1/2, as the bifurcation parameter μ, is driven further away from its critical value, μ _c . In this work, we first characterize analytically this scaling law using complex variable techniques for a family of one-dimensional maps, called the normal form for the saddle-node bifurcation. We then apply our general analytic results to a single-species ecological model with harvesting given by a unimodal map, characterizing the delayed transition and the scaling law arising due to the constant of harvesting. For both analyzed systems, we show that the numerical results are in perfect agreement with the analytical solutions we are providing. The procedure presented in this work can be used to characterize the scaling laws of one-dimensional discrete dynamical systems with saddle-node bifurcations.     

Numerical study of transient laminar natural convection heat transfer over a sphere subjected to a constant heat flux

Transient laminar natural convection over a sphere which is subjected to a constant heat flux has been studied numerically for high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and a wide range of Prandtl numbers (Pr = 0.02, 0.7, 7, and 100). A plume with a mushroom-shaped cap forms above the sphere and drifts upward continuously with time. The size and the level of temperature of the transient cap and plume stem decrease with increasing Gr and Pr. Flow separation and an associated vortex may appear in the wake of the sphere depending on the magnitude of Gr and Pr. A recirculation vortex which appears and grows until “steady state” is attained was found only for the very high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and the lowest Prandtl number considered (Pr = 0.02). The appearance and subsequent disappearance of a vortex was observed for Gr = 10⁹ and Pr = 0.7. Over the lower hemisphere, the thickness of both the hydrodynamic (δ_H) and the thermal (δ_T) boundary layers remain nearly constant and the sphere surface is nearly isothermal. The surface temperature presents a local maximum in the wake of the sphere whenever a vortex is established in the wake of the sphere. The surface pressure recovery in the wake of the sphere increases with decreasing Pr and with increasing Gr. For very small Pr, unlike forced convection, the ratio δ_T/δ_H remains close to unity. The results are in good agreement with experimental data and in excellent agreement with numerical results available in the literature. A correlation has also been presented for the overall Nusselt number as a function of Gr and Pr.     

On the Reynolds number scaling of vorticity production at no-slip walls during vortex-wall collisions

Recently, numerical studies revealed two different scaling regimes of the peak enstrophy Z and palinstrophy P during the collision of a dipole with a no-slip wall [Clercx and van Heijst, Phys. Rev. E 65, 066305, 2002]: Z μ Re^0.8{Z\propto{\rm Re}^{0.8}} and P μ Re^2.25{P\propto {\rm Re}^{2.25}} for 5 × 10² ≤ Re ≤ 2 × 10⁴ and Z μ Re^0.5{Z\propto{\rm Re}^{0.5}} and P μ Re^1.5{P\propto{\rm Re}^{1.5}} for Re ≥ 2 × 10⁴ (with Re based on the velocity and size of the dipole). A critical Reynolds number Re _c (here, Re_c ? 2×10⁴{{\rm Re}_c\approx 2\times 10^4}) is identified below which the interaction time of the dipole with the boundary layer depends on the kinematic viscosity ν. The oscillating plate as a boundary-layer problem can then be used to mimick the vortex-wall interaction and the following scaling relations are obtained: Z μ Re^3/4, P μ Re^9/4{Z\propto{\rm Re}^{3/4}, P\propto {\rm Re}^{9/4}} , and dP/dt μ Re^11/4{\propto {\rm Re}^{11/4}} in agreement with the numerically obtained scaling laws. For Re ≥ Re _c the interaction time of the dipole with the boundary layer becomes independent of the kinematic viscosity and, applying flat-plate boundary-layer theory, this yields: Z μ Re^1/2{Z\propto{\rm Re}^{1/2}} and P μ Re^3/2{P\propto {\rm Re}^{3/2}}.     

Rheo-dielectric behavior of low molecular weight liquid crystal 2. Behavior of 8CB in nematic and smectic states

Rheo-dielectric behavior was examined for 4−4^′−n-octyl-cyanobiphenyl (8CB) having large dipoles parallel to its principal axis (in the direction of the C≡N bond). In the quiescent state at all temperatures (T) examined, orientational fluctuation of the 8CB molecules was observed as dielectric dispersions at characteristic frequencies ω_c>10⁶ s⁻¹. In the isotropic state at high T, no detectable changes of the complex dielectric constant ɛ^*(ω) were found under slow flow at shear rates ˙γ≫ω_c. In the nematic state at intermediate T, the terminal relaxation intensity of ɛ^*(ω) was decreased under such slow flow. In the smectic state at lower T, the flow effect became much less significant. These results were related to the flow-induced changes of the liquid crystalline textures in the nematic and smectic states, and the differences of the rheo-dielectric behavior in these states are discussed in relation to a difference of the symmetry of molecular arrangements in the nematic and smectic textures. Received: 1 October 1998 Accepted: 13 January 1999     

Soret effect and slow mass diffusion as a catalyst for overstability in Marangoni-Bénard flows

We study the onset of time dependent Marangoni-Bénard convection in binary mixtures subject to Soret effect by numerical computation of linear instability thresholds in infinite fluid layers and two-dimensional boxes. The calculations are done for positive Marangoni numbers (Ma > 0) and negative Marangoni Soret parameters S _M = –(D _S γ _c )/(Dγ _T ) where D _S and D are the Soret and mass diffusion coefficients, respectively, and γ _T , γ _c are the first derivatives of the surface tension with respect to temperature and concentration. Our purpose is to understand why for particular choices of Prandtl and Schmidt numbers, the increase of the stabilizing solutal contribution leads to a decrease of the critical temperature difference, a phenomenon already reported by Chen & Chen [5] and Skarda et al. [12] For various choices of Prandtl and Schmidt numbers we analyze the evolution of the critical Marangoni number Ma _c , critical wavenumber k _c and angular frequency ω _c with S _M and compute the corresponding eigenvectors. We next propose a physical mechanism which explains how the stabilizing solutal contribution acts as a catalyst for overstability. Finally, we extend our results to two dimensional boxes of small aspect ratio.