No Arabic abstract
Direct Numerical Simulations are used to solve turbulent flow and heat transfer over a variety of rough walls in a channel. The wall geometries are exactly resolved in the simulations. The aim is to understand the effect of roughness morphology and its scaling on the augmentation of heat transfer relative to that of skin friction. A number of realistic rough surface maps obtained from the scanning of gas turbine blades and internal combustion engines as well as several artificially generated rough surfaces are examined. In the first part of the paper, effects of statistical surface properties, namely surface slope and roughness density, at constant roughness height are systematically investigated, and it is shown that Reynolds analogy factor (two times Stanton number divided by skin friction coefficient) varies meaningfully but moderately with the surface parameters except for the case with extremely low slope or density where the Reynolds analogy factor grows significantly and tends to that of a smooth wall. In the second part of the paper, the roughness height is varied (independently in both inner and outer units) while the geometrical similarity is maintained. Considering all the simulated cases, it is concluded that Reynolds analogy factor correlates fairly well with the equivalent sand roughness scaled in inner units and asymptotically tends to a plateau.
We present a systematic investigation of the effects of roughness geometry on turbulent Rayleigh-Benard convection (RBC) over rough plates with pyramid-shaped and periodically distributed roughness elements. Using a parameter $lambda$ defined as the height of a roughness element over its base width, the heat transport, the flow dynamics and local temperatures are measured for the Rayleigh number range $7.50times 10^{7} leq Raleq 1.31times 10^{11}$, and the Prandtl number $Pr$ from 3.57 to 23.34 at four values of $lambda$. It is found that the heat transport scaling, i.e. $Nusim Ra^{alpha}$ where $Nu$ is the Nusselt number, may be classified into three regimes. In Regime I, the system is in a dynamically smooth state. The heat transport scaling is the same as that in a smooth cell. In Regimes II and III, the heat transport enhances. When $lambda$ is increased from 0.5 to 4.0, $alpha$ increases from 0.36 to 0.59 in Regime II, and it increases from 0.30 to 0.50 in Regime III. The experiment demonstrates the heat transport scaling in turbulent RBC can be manipulated using $lambda$. Previous studies suggest that the transition from Regime I to Regime II, occurs when the thermal boundary layer (BL) thickness becomes smaller than the roughness height $h$. Direct measurements of the viscous BL in the present study suggest that the transition from Regime II to Regime III is likely a result of the viscous BL thickness becoming smaller $h$. The scaling exponent of the Reynolds number $Re$ vs. $Ra$ changes from 0.471 to 0.551 when $lambda$ is increased from 0.5 to 4.0. It is also found that increasing $lambda$ increases the clustering of thermal plumes which effectively increases the plumes lifetime that are ultimately responsible for the enhanced heat transport.
We report heat transfer and temperature profile measurements in laboratory experiments of rapidly rotating convection in water under intense thermal forcing (Rayleigh number $Ra$ as high as $sim 10^{13}$) and unprecedentedly strong rotational influence (Ekman numbers $E$ as low as $10^{-8}$). Measurements of the mid-height vertical temperature gradient connect quantitatively to predictions from numerical models of asymptotically rapidly rotating convection, separating various flow phenomenologies. Past the limit of validity of the asymptotically-reduced models, we find novel behaviors in a regime we refer to as rotationally-influenced turbulence, where rotation is important but not as dominant as in the known geostrophic turbulence regime. The temperature gradients collapse to a Rayleigh-number scaling as $Ra^{-0.2}$ in this new regime. It is bounded from above by a critical convective Rossby number $Ro^*=0.06$ independent of domain aspect ratio $Gamma$, clearly distinguishing it from well-studied rotation-affected convection.
A Direct Numerical Simulation (DNS) of the incompressible flow around a rectangular cylinder with chord-to-thickness ratio 5:1 (also known as the BARC benchmark) is presented. The work replicates the first DNS of this kind recently presented by Cimarelli et al (2018), and intends to contribute to a solid numerical benchmark, albeit at a relatively low value of the Reynolds number. The study differentiates from previous work by using an in-house finite-differences solver instead of the finite-volumes toolbox OpenFOAM, and by employing finer spatial discretization and longer temporal average. The main features of the flow are described, and quantitative differences with the existing results are highlighted. The complete set of terms appearing in the budget equation for the components of the Reynolds stress tensor is provided for the first time. The different regions of the flow where production, redistribution and dissipation of each component take place are identified, and the anisotropic and inhomogeneous nature of the flow is discussed. Such information is valuable for the verification and fine-tuning of turbulence models in this complex separating and reattaching flow.
We numerically investigate turbulent Rayleigh-Benard convection within two immiscible fluid layers, aiming to understand how the layer thickness and fluid properties affect the heat transfer (characterized by the Nusselt number $Nu$) in two-layer systems. Both two- and three-dimensional simulations are performed at fixed global Rayleigh number $Ra=10^8$, Prandtl number $Pr=4.38$, and Weber number $We=5$. We vary the relative thickness of the upper layer between $0.01 le alpha le 0.99$ and the thermal conductivity coefficient ratio of the two liquids between $0.1 le lambda_k le 10$. Two flow regimes are observed: In the first regime at $0.04lealphale0.96$, convective flows appear in both layers and $Nu$ is not sensitive to $alpha$. In the second regime at $alphale0.02$ or $alphage0.98$, convective flow only exists in the thicker layer, while the thinner one is dominated by pure conduction. In this regime, $Nu$ is sensitive to $alpha$. To predict $Nu$ in the system in which the two layers are separated by a unique interface, we apply the Grossmann-Lohse theory for both individual layers and impose heat flux conservation at the interface. Without introducing any free parameter, the predictions for $Nu$ and for the temperature at the interface well agree with our numerical results and previous experimental data.
To understand how internal flow structures manifest themselves in the global heat transfer, we study the correlation between different flow modes and the instantaneous Nusselt number ($Nu$) in a two-dimensional square Rayleigh-Benard convection cell. High-resolution and long-time direct numerical simulations are carried out for Rayleigh numbers between $10^{7}$ and $10^{9}$ and a Prandtl number of 5.3. The investigated Nusselt numbers include the volume-averaged $Nu_{text{vol}}$, the wall-averaged $Nu_{text{wall}}$, the kinetic energy dissipation based $Nu_{text{kinetic}}$, and the thermal energy dissipation based $Nu_{text{thermal}}$. The Fourier mode decomposition and proper orthogonal decomposition are adopted to extract the coherent flow structure. Our results show that the single-roll mode, the horizontally stacked double-roll mode, and the quadrupolar flow mode are more efficient for heat transfer on average. In contrast, the vertically stacked double-roll mode is inefficient for heat transfer on average. The volume-averaged $Nu_{text{vol}}$ and the kinetic energy dissipation based $Nu_{text{kinetic}}$ can better reproduce the correlation of internal flow structures with heat transfer efficiency than that of the wall-averaged $Nu_{text{wall}}$ and the thermal energy dissipation based $Nu_{text{thermal}}$, even though these four Nusselt numbers give consistent time-averaged mean values. The ensemble-averaged time trace of $Nu$ during flow reversal shows that only the volume-averaged $Nu_{text{vol}}$ can reproduce the overshoot phenomena that is observed in the previous experimental study. Our results reveal that the proper choice of $Nu$ is critical to obtain a meaningful interpretation.