Large eddy simulation of a confined square cavity with natural convection based on compressible flow equations

Large eddy simulation of natural convection in a confined square cavity is described. The use of a complex compressible code with an artificial acoustic stiffness correction method, allows the use of higher time steps for a faster time and statistical convergence. We consider a broadly studied experimental case, consisting of a natural convective flow in a confined square cavity, with vertical walls heated at different rates (active walls), set at Ra = 1.58 × 10⁹. Turbulent boundary layers developing on the active walls and a vertical stable stratification characterize the mean flow. It is shown here that the results of this study match the experimental results reported in literature; for instance, mean velocity results. Although results for rms velocity fluctuations are barely over-predicted, the peak region is properly represented, while the greatest disagreements are found in the turbulent heat flow rate (velocity–temperature correlations). Turbulent structures were identified using different visualization methods and statistical studies. The authors found that the boundary layers on the active walls almost reach the fully turbulent regime, tending toward the laminar regime along the horizontal walls.     

Manipulation of the reverse-flow region downstream of a fence by spanwise vortices

The experimental investigation of a turbulent separated flow over a fence is presented. By introducing a periodic disturbance upstream of the separation region in front of the fence, the time averaged length of the separation region downstream of the fence was reduced by as much as 40%. Two types of flow manipulation were applied: an oscillating cross-flow with zero net mass-flux through a spanwise slot in the floor of the test section and a spanwise oriented, oscillating spoiler. The cross-flow was generated by a loudspeaker system connected to a chamber underneath the spanwise slot. Both types of flow manipulation generate spanwise vortices at the fence that convect into the region downstream of the fence where they enhance the mixing in the shear layer and reduce the time mean length of the reverse-flow region downstream of the fence. Velocity profiles phase averaged with respect to the forcing frequency and phase triggered flow visualisations show that the spanwise vortices cause the long reverse-flow region of the unmanipulated flow to break up into separate smaller regions. While the time mean length of the reverse-flow region is reduced in the manipulated case, the length of the region where instantaneous reverse-flow occurs is not changed. The data presented include wall pulsed-wire measurements of the wall shear-stress and its turbulent fluctuations, and LDA measurements of the streamwise and the wall-normal velocity components and turbulent stresses.     

Experimental characterization of non-premixed turbulent jet propane flames

This paper reports an experimental study conducted on turbulent jet propane flames aiming at further understanding of turbulent structure in non-premixed slow-chemistry combustion systems. Measurements of mean and fluctuating velocity and temperature fields, mean concentration of major chemical species, correlation between velocity and temperature fluctuations, and dissipation of temperature fluctuations are reported in a turbulent round jet non-premixed propane flame, Re=20 400 and 37 600, issuing vertically in still air. The experimental conditions were designed to provide a complete definition of the upstream boundary conditions in the measurement domain for the purpose of validating computational models. The measured data depicts useful flow field information for describing turbulent non-premixed slow-chemistry flames. Velocity–temperature correlation measurements show turbulent heat fluxes tended to be restricted to the mixing layer where large temperature gradients occurred. Observations of non-gradient diffusion of heat at x/D=10 were verified. Temperature fluctuation dissipation, χ, showed the highest values in the shear layer, where the variance of temperature fluctuations was maximum and combustion occurred. The isotropy between the temperature dissipation in the radial and tangential directions was confirmed. By contrast, the observed anisotropy between axial and radial directions of dissipation suggests the influence of large structures in the entrainment shear layer on the production of temperature fluctuations in the flame region. The value of the normalized scalar dissipation at the stoichiometric mixture fraction surface, χ_st, was calculated, and ranges between 2 and 4 s⁻¹. The measured data were used to estimate the budgets in the balance equations for turbulent kinetic energy, Reynolds shear stresses, turbulent heat flux and temperature variance, quantifying the mechanisms involved in the generation of turbulence as well as in the transport of the temperature.     

Higher-order and length-scale statistics of velocity and temperature fluctuations in turbulent boundary layer along a heated vertical flat plate

Time-developing direct numerical simulation (DNS) was performed to clarify the higher-order turbulent behaviors in the thermally-driven boundary layers both in air and water along a heated vertical flat plate. The predicted statistics of the heat transfer rates and the higher-order turbulent behaviors such as skewness factors, flatness factors and spatial correlation coefficients of the velocity and temperature fluctuations in the natural-convection boundary layer correspond well with those obtained from experiments for space-developing flows. The numerical results reveal that the turbulent structures of the buoyancy-driven boundary layers are mainly controlled by the fluid motions in the outer region of the boundary layer, and these large-scale structures are strongly connected with the generation of turbulence in the thermally-driven boundary layers, in accordance with the actual observations for space-developing flows. Moreover, to specify the turbulence structures of the boundary layers, the cross-correlation coefficients and the characteristic length scales are examined for the velocity and thermal fields. Consequently, it is found that with a slight increase in freestream velocity, the cross-correlation coefficient for the Reynolds shear stress and turbulent heat flux increases for opposing flow and decreases for aiding flow, and the integral scales for the velocity and temperature fields become larger for opposing flow and smaller for aiding flow compared with those for the pure natural-convection boundary layer.     

Structure of turbulent two-dimensional channel flow with strongly heated wall

A laser Doppler velocimeter and a resistance thermometer were used to study velocity and temperature statistics in a strongly heated turbulent two-dimensional channel flow, with the wall temperature up to 700 °C and a Reynolds number of 14,000. Normalized mean velocity and mean temperature profiles were not significantly affected by the wall heating. Turbulent intensities of temperature fluctuation were also insensitive to the heat flux. However, turbulent intensities of velocity fluctuation were suppressed in the region away from the wall, whereas those near the wall were not changed noticeably by the wall heating. This phenomenon was explained by the balance of three parameters: turbulent production, viscous dissipation and intermittency.     

Near-field measurements and development of a new boundary layer over a flat plate with localized suction

Suction on a turbulent boundary layer is applied through a narrow strip in order to understand the effects suction can have on the boundary layer development and turbulent structures in the flow. Detailed two-component laser Doppler velocimetry (LDV) and laser-induced fluorescence (LIF) based measurements have been undertaken in regions close to the suction strip and further downstream. The region close to the strip involves a flow reversal accompanied by a change in sign for the Reynolds shear stress and strong gradients in the flow variables. The mean streamwise velocity after suction remains larger than its corresponding no-suction value. Relative to the no-suction case, the velocity fluctuations first decrease with suction followed by a slow recovery which may involve a slight overshoot. LIF visualizations indicate that compared to the no-suction case, the low-speeds streaks stay closer to the wall and exhibit a smaller amount of spanwise and wall-normal oscillations with suction. The visualization results are consistent with two-point velocity correlation measurements. The streamwise and spanwise correlation measurements indicate that the structures are disrupted or removed from the boundary layer due to suction suggesting that the original boundary layer has been strongly influenced by suction. The results are explained by the development of a new inner layer that forms downstream of the suction strip.     

Turbulent free convection in a vertical converging channel

This paper presents the results of experimental and numerical investigations of the problem of turbulent natural convection in a converging-plate vertical channel. The channel has two isothermally heated inclined walls and two adiabatic vertical side walls. The parameters involved in this study are the channel geometry represented by the channel width at exit, the inclination of the heated walls and the temperature difference between the heated walls and the ambient. The investigation covered modified Rayleigh numbers up to 10⁸ in the computational study and up to 9.3 × 10⁶ in the experimental work. The experimental measurements focused on the velocity field and were carried out using a PIV system and included measurements of the mean velocity profiles as well as the root-mean-square velocity and shear stress profiles. The experiments were conducted for an inclination angle of 30°, a gap width of 10 mm and two temperature differences (∆T=25.4°C and 49.8°C). The velocity profiles in the lower part of the channel indicated the presence of two distinct layers. The first layer is adjacent to the heated plate and driven by buoyancy forces while the second layer extends from the point of maximum velocity to the channel center plane and driven mainly by shear forces. The velocity profile at the upper portion of the channel has shown the merging of the two boundary layers growing over the two heated walls. The measured values of the Reynolds shear stress and root mean square of the horizontal and vertical velocity fluctuation components have reached their maximum near the wall while having smaller values in the core region. The computational results have shown that the average Nusselt number increases approximately linearly with the increase of the modified Rayleigh number when plotted on log–log scale. The variation of the local Nusselt number indicated infinite values at the channel inlet (leading edge effect) and high values at the channel exit (trailing edge effect). For a fixed value of the top channel opening, the increase of the inclination angle tended to reduce flow velocity at the inlet section while changing the flow structure near the heated plates in such a way to create boundary-layer type flow. The maximum value of the average Nusselt number occurs when θ = 0 and decreases with the increase of the inclination angle. On the other hand, the increase of the channel width at exit for the same inclination angle caused a monotonic increase in the flow velocity at the channel inlet.     

Near wall characterization of the flow over a two-dimensional steep smooth hill

Near-wall data for the strongly perturbed flow in a neutrally stable boundary layer encountering a steep, smooth, two-dimensional hill are presented. Observations were made on the centerplane of a water channel at thirteen stations relative to the hill by laser Doppler anemometry. The large reverse flow region that is formed on the lee of the hill was particularly scrutinized through seven measuring stations. Results are presented for the mean and turbulent properties of the flow. Wall shear stress was evaluated through fitting procedures that resorted to the near wall behavior of the velocity profile. Logarithmic fits as well as predictions through the Reynolds stress profiles are also presented.     

Turbulent mixed convection flow over a backward-facing step––the effect of the step heights

Effect of the backward-facing step heights on turbulent mixed convection flow along a vertical flat plate is examined experimentally. The step geometry consists of an adiabatic backward-facing step, an upstream wall and a downstream wall. Both the upstream and downstream walls are heated to a uniform and constant temperature. Laser–Doppler velocimeter and cold wire anemometer were used, respectively, to measure simultaneously the time-mean velocity and temperature distributions and their turbulent fluctuations. The experiment was carried out for step heights of 0, 11, and 22 mm, at a free stream air velocity, u_∞, of 0.41 m/s, and a temperature difference, ΔT, of 30 °C between the heated walls and the free stream air. The present results reveal that the turbulence intensity of the streamwise and transverse velocity fluctuations and the intensity of temperature fluctuations downstream of the step increase as the step height increases. Also, it was found that both the reattachment length and the heat transfer rate from the downstream heated wall increase with increasing step height.     

Shear flow disturbances due to interaction with convective rolls

The two-dimensional stationary problem of convection in a two-component stratified medium (for example, humid air) over a nonuniformly heated horizontal surface is solved in the linear approximation. In the case in which the surface temperature is harmonically dependent on the horizontal coordinate, convective rolls are formed above the surface. Depending on the background temperature stratification, situations characterized by temperature response sign reversal can occur (for example, a decrease in the temperature of the medium as a response to heating from beneath, or “negative” heat capacity of the medium). The disturbances arising in the background horizontal flow with vertical shear directed along the convective rolls are also calculated.     

Experiments in Turbulent Pipe Flow with Polymer Additives at Maximum Drag Reduction

In this paper we report on (two-component) LDV experiments in a fully developed turbulent pipe flow with a drag-reducing polymer (partially hydrolyzed polyacrylamide) dissolved in water. The Reynolds number based on the mean velocity, the pipe diameter and the local viscosity at the wall is approximately 10000. We have used polymer solutions with three different concentrations which have been chosen such that maximum drag reduction occurs. The amount of drag reduction found is 60–70%. Our experimental results are compared with results obtained with water and with a very dilute solution which exhibits only a small amount of drag reduction. We have focused on the observation of turbulence statistics (mean velocities and turbulence intensities) and on the various contributions to the total shear stress. The latter consists of a turbulent, a solvent (viscous) and a polymeric part. The polymers are found to contribute significantly to the total stress. With respect to the mean velocity profile we find a thickening of the buffer layer and an increase in the slope of the logarithmic profile. With respect to the turbulence statistics we find for the streamwise velocity fluctuations an increase of the root mean square at low polymer concentration but a return to values comparable to those for water at higher concentrations. The root mean square of the normal velocity fluctuations shows a strong decrease. Also the Reynolds (turbulent) shear stress and the correlation coefficient between the stream wise and the normal components are drastically reduced over the entire pipe diameter. In all cases the Reynolds stress stays definitely non-zero at maximum drag reduction. The consequence of the drop of the Reynolds stress is a large polymer stress, which can be 60% of the total stress. The kinetic-energy balance of the mean flow shows a large transfer of energy directly to the polymers instead of the route by turbulence. The kinetic energy of the turbulence suggests a possibly negative polymeric dissipation of turbulent energy. This revised version was published online in July 2006 with corrections to the Cover Date.     

Turbulent thermal convection in wide horizontal fluid layers

Turbulence in thermal convection is investigated for flows in which the production of turbulence energy is due solely to buoyancy, and the statistics of the flow are homogeneous in horizontal planes. New experimental results for high Rayleigh number unsteady turbulent convection in a horizontal layer heated from below and insulated from above are presented and compared to turbulent Rayleigh convection, convection in the planetary boundary layer, and laboratory penetrative convection. Mean temperature fields are correlated in terms of wall layer scales and convection scales. Joint statistics of turbulent temperature and horizontal velocity and vertical velocity through fourth order are presented for the core region of the convection layer.This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

EXPERIMENTAL STUDY OF MEASUREMENT FOR DISSIPATION RATE SCALING EXPONENT IN HEATED WALL TURBULENCE

Experimental investigations have been devoted to the study of scaling law of coarse-grained dissipation rate structure function for velocity and temperature fluctuation of non-isotropic and inhomogeneous turbulent flows at moderate Reynolds number. Much attention has been paid to the case of turbulent boundary layer, which is typically the non-istropic and inhomogeneous trubulence because of the dynamically important existence of organized coherent structure burst process in the near wall region . Longitudinal velocity and temperature have been measured at different vertical positions in turbulent boundary layer over a heated and unheated flat plate in a wind tunnel using hot wire anemometer. The influence of non-isotropy and inhomogeneity and heating the wall on the scaling law of the dissipation rate structure function is studied because of the existence of organized coherent structure burst process in the near wall region . The scaling law of coarse-grained dissipation rate structure function is foun     

Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation(DNS) of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed.The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers(P r = 0.71,1.5,and 3.0) and a shear Reynolds number(Reτ = 180).Some typical thermal statistics,including normalized mean temperature and their fluctuations,turbulent heat fluxes,Nusselt number and so on,are analyzed.The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number.Two reasons can explain this.First,the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases,and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field.Second,the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number,and the thermal feedback of the particles to turbulence becomes weak.     

Nonlinear Dynamics of Turbulent Thermals in Shear Flow

The nonlinear integral model of a turbulent thermal is extended to the case of the horizontal component of its motion relative to the medium (e.g., thermal floating-up in shear flow). In contrast to traditional models, the possibility of a heat source in the thermal is taken into account. For a piecewise constant vertical profile of the horizontal velocity of the medium and a constant vertical velocity shear, analytical solutions are obtained which describe different modes of dynamics of thermals. The nonlinear interaction between the horizontal and vertical components of thermal motion is studied because each of the components influences the rate of entrainment of the surrounding medium, i.e., the growth rate of the thermal size and, hence, its mobility. It is shown that the enhancement of the entrainment of the medium due to the interaction between the thermal and the cross flow can lead to a significant decrease in the mobility of the thermal.     

Velocity-temperature correlation in strongly heated channel flow

Velocity-temperature correlations in a strongly heated channel flow were investigated experimentally by a LDV and a resistance thermometer. The wall heat flux is varied up to 50,000 W/m² with reference mean-velocity of 15 m/s, and then, the wall temperature reaches up to 1,000 K. The results show that the ejection fluid motion is intensified by the strong heating near the wall increasing the turbulent heat flux from the wall. The intensification of ejection motion balances the destruction of turbulent heat flux. Then, the overall turbulent heat transfer does not change clearly. Part of this work was supported by Tanikawa Fund, Promotion of Thermal Technology     

A Direct Numerical Study of Transport and Anisotropy in a Stably Stratified Turbulent Flow with Uniform Horizontal Shear

Direct numerical simulations of homogeneous turbulence in stably stratified shear flow have been performed to aid the understanding of turbulence and turbulent mixing in geophysical flow. Two cases are compared. In the first case, which has been studied in the past, the mean velocity has vertical shear and the mean density is vertically stably stratified. In the second case, which has not been studied systematically before, the mean velocity has horizontal shear and the mean density is again vertically stably stratified. The critical value of the gradient Richardson number, for which a constant turbulence level is obtained, is found to be an order of magnitude larger in the horizontal shear case. The turbulent transport coefficients of momentum and vertical mass transfer are also an order of magnitude larger in the horizontal shear case. The anisotropy of the turbulence intensities are found to be in the range expected of flows with mean shear with no major qualitative change in the range of Richardson numbers studied here. However, the anisotropy of the turbulent dissipation rate is strongly affected by stratification with the vertical component dominating the others. This revised version was published online in July 2006 with corrections to the Cover Date.     

Temperature and velocity field characteristics of turbulent natural convection in a vertical parallel-plate channel with asymmetric heating

Turbulent natural convection in an asymmetrically heated vertical parallel-plate channel has been studied experimentally and numerically using LDA and CFD. Simultaneous velocity and temperature measurements across the channel at different elevations have been carried out. Three different Ra(b/h) values of 1.91 × 10⁷, 2.74 × 10⁷ and 3.19 × 10⁷ are considered with the channel aspect ratio of b/h = 1/20. Experimental and numerical data are presented in the form of streamwise direction heated wall surface temperature, mean velocity, mean temperature, Reynolds shear stress and turbulent kinetic energy profiles along the channel for one case. These profiles exhibit the flow field development along the channel emphatically. The numerical technique used predicts temperature field fairly well, considerably over-estimating velocity field in the core region.     

Turbulent heat transfer characteristics of water flow in a rotating pipe

Large-Eddy-Simulation of turbulent heat transfer for water flow in rotating pipe is performed, for various rotation ratios (0 ≤ N ≤ 14). The value of the Reynolds number, based on the bulk velocity and pipe diameter, is Re = 5,500. The aim of this study is to examine the effect of the rotating pipe on the turbulent heat transfer for water flow, as well as the reliability of the LES approach for predicting turbulent heat transfer in water flow. Some predictions for the case of non-rotating pipe are compared to the available results of literature for validation. To depict the influence of the rotation ratio on turbulent heat transfer, many statistical quantities are analyzed (distributions of mean temperature, rms of fluctuating temperature, turbulent heat fluxes, higher-order statistics). Some contours of instantaneous temperature fluctuations are examined.     

A method for separating shock wave motion and turbulence in LDV measurements

 Two-component laser Doppler velocimetry (LDV) measurements were made in a planar, two-dimensional flow containing an unsteady oblique shock wave formed by the convergence of two supersonic streams past a thick plate. High-speed wall pressure measurements locate the shock wave and, consequently, allow separation of the effects of shock wave motion from the turbulence fluctuations in the LDV measurements of the shock-separated free shear layer. In the current flow isolating the large-scale changes in the position of the shock from the turbulence primarily reduces the experimental scatter rather than significantly changing the shapes or magnitudes of the turbulent stress profiles. Changes in the direction of shock motion do not significantly affect the mean velocity, but do affect the turbulent stresses. Received: 11 August 1997/Accepted: 30 September 1998