Experimental and numerical investigation of turbulent convection in a rotating cylinder

2009 ◽  
Vol 642 ◽  
pp. 445-476 ◽  
Author(s):  
R. P. J. KUNNEN ◽  
B. J. GEURTS ◽  
H. J. H. CLERCX

The effects of an axial rotation on the turbulent convective flow because of an adverse temperature gradient in a water-filled upright cylindrical vessel are investigated. Both direct numerical simulations and experiments applying stereoscopic particle image velocimetry are performed. The focus is on the gathering of turbulence statistics that describe the effects of rotation on turbulent Rayleigh–Bénard convection. Rotation is an important addition, which is relevant in many geophysical and astrophysical flow phenomena.A constant Rayleigh number (dimensionless strength of the destabilizing temperature gradient) Ra = 109 and Prandtl number (describing the diffusive fluid properties) σ = 6.4 are applied. The rotation rate, given by the convective Rossby number Ro (ratio of buoyancy and Coriolis force), takes values in the range 0.045 ≤ Ro ≤ ∞, i.e. between rotation-dominated flow and zero rotation. Generally, rotation attenuates the intensity of the turbulence and promotes the formation of slender vertical tube-like vortices rather than the global circulation cell observed without rotation. Above Ro ≈ 3 there is hardly any effect of the rotation on the flow. The root-mean-square (r.m.s.) values of vertical velocity and vertical vorticity show an increase when Ro is lowered below Ro ≈ 3, which may be an indication of the activation of the Ekman pumping mechanism in the boundary layers at the bottom and top plates. The r.m.s. fluctuations of horizontal and vertical velocity, in both experiment and simulation, decrease with decreasing Ro and show an approximate power-law behaviour of the shape Ro0.2 in the range 0.1 ≲ Ro ≲ 2. In the same Ro range the temperature r.m.s. fluctuations show an opposite trend, with an approximate negative power-law exponent Ro−0.32. In this Rossby number range the r.m.s. vorticity has hardly any dependence on Ro, apart from an increase close to the plates for Ro approaching 0.1. Below Ro ≈ 0.1 there is strong damping of turbulence by rotation, as the r.m.s. velocities and vorticities as well as the turbulent heat transfer are strongly diminished. The active Ekman boundary layers near the bottom and top plates cause a bias towards cyclonic vorticity in the flow, as is shown with probability density functions of vorticity. Rotation induces a correlation between vertical vorticity and vertical velocity close to the top and bottom plates: near the top plate downward velocity is correlated with positive/cyclonic vorticity and vice versa (close to the bottom plate upward velocity is correlated with positive vorticity), pointing to the vortical plumes. In contrast with the well-mixed mean isothermal bulk of non-rotating convection, rotation causes a mean bulk temperature gradient. The viscous boundary layers scale as the theoretical Ekman and Stewartson layers with rotation, while the thermal boundary layer is unaffected by rotation. Rotation enhances differences in local anisotropy, quantified using the invariants of the anisotropy tensor: under rotation there is strong turbulence anisotropy in the centre, while near the plates a near-isotropic state is found.

2018 ◽  
Vol 856 ◽  
pp. 135-168 ◽  
Author(s):  
S. T. Salesky ◽  
W. Anderson

A number of recent studies have demonstrated the existence of so-called large- and very-large-scale motions (LSM, VLSM) that occur in the logarithmic region of inertia-dominated wall-bounded turbulent flows. These regions exhibit significant streamwise coherence, and have been shown to modulate the amplitude and frequency of small-scale inner-layer fluctuations in smooth-wall turbulent boundary layers. In contrast, the extent to which analogous modulation occurs in inertia-dominated flows subjected to convective thermal stratification (low Richardson number) and Coriolis forcing (low Rossby number), has not been considered. And yet, these parameter values encompass a wide range of important environmental flows. In this article, we present evidence of amplitude modulation (AM) phenomena in the unstably stratified (i.e. convective) atmospheric boundary layer, and link changes in AM to changes in the topology of coherent structures with increasing instability. We perform a suite of large eddy simulations spanning weakly ($-z_{i}/L=3.1$) to highly convective ($-z_{i}/L=1082$) conditions (where$-z_{i}/L$is the bulk stability parameter formed from the boundary-layer depth$z_{i}$and the Obukhov length $L$) to investigate how AM is affected by buoyancy. Results demonstrate that as unstable stratification increases, the inclination angle of surface layer structures (as determined from the two-point correlation of streamwise velocity) increases from$\unicode[STIX]{x1D6FE}\approx 15^{\circ }$for weakly convective conditions to nearly vertical for highly convective conditions. As$-z_{i}/L$increases, LSMs in the streamwise velocity field transition from long, linear updrafts (or horizontal convective rolls) to open cellular patterns, analogous to turbulent Rayleigh–Bénard convection. These changes in the instantaneous velocity field are accompanied by a shift in the outer peak in the streamwise and vertical velocity spectra to smaller dimensionless wavelengths until the energy is concentrated at a single peak. The decoupling procedure proposed by Mathiset al.(J. Fluid Mech., vol. 628, 2009a, pp. 311–337) is used to investigate the extent to which amplitude modulation of small-scale turbulence occurs due to large-scale streamwise and vertical velocity fluctuations. As the spatial attributes of flow structures change from streamwise to vertically dominated, modulation by the large-scale streamwise velocity decreases monotonically. However, the modulating influence of the large-scale vertical velocity remains significant across the stability range considered. We report, finally, that amplitude modulation correlations are insensitive to the computational mesh resolution for flows forced by shear, buoyancy and Coriolis accelerations.


2014 ◽  
Vol 71 (2) ◽  
pp. 615-639 ◽  
Author(s):  
Minoru Chikira

Abstract The eastward-propagating intraseasonal oscillation represented by the Chikira–Sugiyama cumulus scheme in a general circulation model was investigated focusing on the variation of the free-tropospheric humidity. The net effect of the vertical advection and cloud process amplifies the positive moisture anomaly in the mature phase, supporting the moisture-mode theory. The horizontal advection causes the eastward propagation of the field. The variation of the moisture profile is accurately understood by using environmental vertical velocity outside cumuli. The velocity is regulated by a thermodynamic balance under a weak temperature gradient. A nondimensional parameter α plays an important role in the moisture variation, which characterizes the efficiency of moistening (drying) induced by external heating (cooling). In the middle and lower troposphere, the major moistening factor is the radiative warming anomaly, which induces the upward environmental vertical velocity anomaly. The reevaporation of the precipitation works as drying, since its cooling effect induces the downward environmental vertical velocity anomaly. Snow melting significantly cools and thereby dries the midtroposphere. The moistening of the midtroposphere is important for moistening the lower troposphere through the reduction of α. The efficiency of moistening depends on the heating profile, and congestus clouds play an important role in it. The heating profile, which maximizes the moistening of the free troposphere, is realized in the mature phase. The atmosphere is marginally unstable even in the mature phase, which is a favorable condition for the congestus clouds to occur.


Author(s):  
Sepideh Kavousi ◽  
Dorel Moldovan

Using phase field modeling simulation approach we investigate the effect of various parameters on the primary and secondary dendrite arm spacing during directional solidification in a single component system. In previous studies the effect of temperature gradient was assumed to be negligible in the transversal directions with a temperature rate equal to the product of thermal gradient and solidification rate. In our study the temperature field is obtained from energy conservation equation by considering the balance of latent heat released in the regions where solidification occurs and energy dissipation due to directional temperature gradient as boundary condition. In our simulations, we implemented a numerical method that enables the investigation of solidification in larger domains. Specifically, the temperature and the order parameter equations are solved only in the domains close to the solidification front; approach that reduces the computational costs significantly. We investigate the interplay and the effect of thermal gradient, solidification rate, undercooling temperature, and the cooling heat flux on arm spacing. By using a well-established power law relation the primary and secondary arm spacing are calculated for various solidification parameters. We also show that, for large heat fluxes, the secondary arm spacing is almost constant for different undercooling temperatures; behavior that demonstrates the need for correction of the power law relation by including the effect of heat flux.


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.


2019 ◽  
Vol 870 ◽  
pp. 460-491 ◽  
Author(s):  
Jérémy Basley ◽  
Laurent Perret ◽  
Romain Mathis

The influence of a cube-based canopy on coherent structures of the flow was investigated in a high Reynolds number boundary layer (thickness $\unicode[STIX]{x1D6FF}\sim 30\,000$ wall units). Wind tunnel experiments were conducted considering wall configurations that represent three idealised urban terrains. Stereoscopic particle image velocimetry was employed using a large field of view in a streamwise–spanwise plane ($0.55\unicode[STIX]{x1D6FF}\times 0.5\unicode[STIX]{x1D6FF}$) combined to two-point hot-wire measurements. The analysis of the flow within the inertial layer highlights the independence of its characteristics from the wall configuration. The population of coherent structures is in agreement with that of smooth-wall boundary layers, i.e. consisting of large- and very-large-scale motions, sweeps and ejections, as well as smaller-scale vortical structures. The characteristics of vortices appear to be independent of the roughness configuration while their spatial distribution is closely linked to large meandering motions of the boundary layer. The canopy geometry only significantly impacts the wall-normal exchanges within the roughness sublayer. Bi-dimensional spectral analysis demonstrates that wall-normal velocity fluctuations are constrained by the presence of the canopy for the densest investigated configurations. This threshold in plan area density above which large scales from the overlying boundary layer can penetrate the roughness sublayer is consistent with the change of the flow regime reported in the literature and constitutes a major difference with flows over vegetation canopies.


Author(s):  
H. Hattori ◽  
Y. Nagano

Direct numerical simulations (DNS) of boundary layers with various thermal stratifications are carried out to investigate the turbulent structures of these flows. The present DNSs quantitatively provide the characteristics of thermally stratified turbulent boundary layers. In particular, the counter gradient diffusion phenomenon is found in a strong, stable stratified boundary layer. On the other hand, in order to adequately predict turbulent boundary layers with various thermal stratifications, an appropriate turbulence model should be employed in the calculation. Thus, using a database obtained by DNS, the strict assessment of turbulent heat transfer model is made so as to construct a reliable advanced turbulence model. The results of in-depth turbulent model evaluation are indicated, in which we have explored the prediction potential of the proposed nonlinear eddy diffusivity models for momentum and heat in both stable and unstable stratified boundary layers.


2014 ◽  
Vol 753 ◽  
pp. 360-401 ◽  
Author(s):  
R. Vicquelin ◽  
Y. F. Zhang ◽  
O. Gicquel ◽  
J. Taine

AbstractThe role of radiative energy transfer in turbulent boundary layers is carefully analysed, focusing on the effect on temperature fluctuations and turbulent heat flux. The study is based on direct numerical simulations (DNS) of channel flows with hot and cold walls coupled to a Monte-Carlo method to compute the field of radiative power. In the conditions studied, the structure of the boundary layers is strongly modified by radiation. Temperature fluctuations and turbulent heat flux are reduced, and new radiative terms appear in their respective balance equations. It is shown that they counteract turbulence production terms. These effects are analysed under different conditions of Reynolds number and wall temperature. It is shown that collapsing of wall-scaled profiles is not efficient when radiation is considered. This drawback is corrected by the introduction of a radiation-based scaling. Finally, the significant impact of radiation on turbulent heat transfer is studied in terms of the turbulent Prandtl number. A model for this quantity, based on the new proposed scaling, is developed and validated.


1979 ◽  
Vol 95 (1) ◽  
pp. 97-118 ◽  
Author(s):  
Fritz H. Bark ◽  
Lennart S. Hultgren

The effects of thermally insulating boundaries on rapidly and almost rigidly rotating gas flows are examined. It is shown that, on a thermally insulating boundary, all boundary layers disappear to zeroth order and that the geostrophic flow alone satisfies the kinematical boundary condition on such a boundary. The temperature gradient of the geostrophic flow is on a horizontal thermally insulating boundary corrected by a weak Ekman layer of strength E½ where E is the Ekman number. On a vertical thermally insulating boundary, the temperature gradient of the geostrophic flow is in the general case corrected by E¼ and $E^{\frac{1}{3}}$ Stewartson layers of strengths E¼ and $E^{\frac{1}{3}}$ respectively.


Sign in / Sign up

Export Citation Format

Share Document