Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter=9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 108. Under these conditions the convection is turbulent, and the time-averaged velocity at any point is ‘zero’. The Reynolds number based on the Taylor microscale, Reλ, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6–7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as Ra1/2Sc1/2, and the Reynolds number would scale as Ra1/2Sc−1/2. The velocity and the flux measurements appear to be consistent with the Ra1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are much higher compared to what would be obtained in Rayleigh–Bénard (R–B) convection for similar density differences.

Numerical prediction of nozzle flow separation: Issue of turbulence modeling

Numerical simulation of separated flows in rocket nozzles is challenging because existing turbulence models are unable to predict it correctly. This paper addresses this issue with the Spalart-Allmaras and Shear Stress Transport (SST) eddy-viscosity models, which predict flow separation with moderate success. Their performances have been compared against experimental data for a conical and two contoured subscale nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the nozzle pressure ratio (NPR) and nozzle type. A careful assessment indicated how the model had to be tuned for better, consistent prediction. It is learnt that SST model's failure is caused by limiting of the shear stress inside boundary layer according to Bradshaw's assumption, and by over prediction of jet spreading rate. Accordingly, SST's coefficients were empirically modified to match the experimental wall pressure data. Results confirm that accurate RANS prediction of separation depends on the correct capture of the jet spreading rate, and that it is feasible over a wide range of NPRs by modified values of the diffusion coefficients in the turbulence model. (C) 2015 Elsevier Masson SAS. All rights reserved.

Turbulence measurements in an equilibrium axisymmetric wall jet

This paper reports measurements of turbulent quantities in an axisymmetric wall jet subjected to an adverse pressure gradient in a conical diffuser, in such a way that a suitably defined pressure-gradient parameter is everywhere small. Self-similarity is observed in the mean velocity profile, as well as the profiles of many turbulent quantities at sufficiently large distances from the injection slot. Autocorrelation measurements indicate that, in the region of turbulent production, the time scale of ν fluctuations is very much smaller than the time scale of u fluctuations. Based on the data on these time scales, a possible model is proposed for the Reynolds stress. One-dimensional energy spectra are obtained for the u, v and w components at several points in the wall jet. It is found that self-similarity is exhibited by the one-dimensional wavenumber spectrum of $\overline{q^2}(=\overline{u^2}+\overline{v^2}+\overline{w^2})$, if the half-width of the wall jet and the local mean velocity are used for forming the non-dimensional wavenumber. Both the autocorrelation curves and the spectra indicate the existence of periodicity in the flow. The rate of dissipation of turbulent energy is estimated from the $\overline{q^2}$ spectra, using a slightly modified version of a previously suggested method.

Mean Flow Measurements in Turbulent Boundary Layers along Mildly Curved Surfaces

An experimental investigation of the mean flow characteristics of two-dimensional turbulent boundary layers over surfaces of mild longitudinal curvature is reported. The study covered both convex and concave walls of \d/Rw I « 0.013 (d being the boundary-layer thickness and Rw being the wall radius). It was found that, whereas the region close to the wall was not affected significantly by wall curvature, the outer region was very sensitive to even mild wall curvature. A detailed study of the wake region using present and other available data suggests a systematic effect of b/Rw on the wake structure. The paper also discusses in detail the effect of mild wall curvature on the boundary-layer development with particular emphasis on the difference in behavior of the boundary layer at short and long distances from the leading edge of the curved wall, an aspect which has not received sufficient attention in previous experimental investigations. An attempt has been made to explain this behavior from a consideration of the structure of turbulence in boundary layers over curved surfaces taken into account.

Experiments on the fine structure of turbulence

Investigations have been carried out of some aspects of the fine-scale structure of turbulence in grid flows, in boundary layers in a zero pressure gradient and in a boundary layer in a strong favourable pressure gradient leading to relaminarization. Using a narrow-band filter with suitable mid-band frequencies, the properties of the fine-scale structure (appearing as high frequency pulses in the filtered signal) were analysed using the variable discriminator level technique employed earlier by Rao, Narasimha & Badri Narayanan (1971). It was found that, irrespective of the type of flow, the characteristic pulse frequency (say Np) defined by Rao et al. was about 0·6 times the frequency of the zero crossings. It was also found that, over the small range of Reynolds numbers tested, the ratio of the width of the fine-scale regions to the Kolmogorov scale increased linearly with Reynolds number in grid turbulence as well as in flat-plate boundarylayer flow. Nearly lognormal distributions were exhibited by this ratio as well as by the interval between successive zero crossings. The values of Np and of the zero-crossing rate were found to be nearly constant across the boundary layer, except towards its outer edge and very near the wall. In the zero-pressure-gradient boundary-layer flow, very near the wall the high frequency pulses were found to occur mostly when the longitudinal velocity fluctuation u was positive (i.e. above the mean), whereas in the outer part of the boundary layer the pulses more often occurred when u was negative. During acceleration this correlation between the fine-scale motion and the sign of u was less marked.

Boundary-layer effects on internal flow choking in dual-thrust solid rocket motors

Theoretical studies have been carried out to examine internal flow choking in the inert simulators of a dual-thrust motor. Using a two-dimensional k-omega turbulence model, detailed parametric studies have been carried out to examine aerodynamic choking and the existence of a fluid throat at the transition region during the startup transient of dual-thrust motors. This code solves standard k-omega turbulence equations with shear flow corrections using a coupled second-order-implicit unsteady formulation. In the numerical study, a fully implicit finite volume scheme of the compressible, Reynolds-averaged, Navier-Stokes equations is employed. It was observed that, at the subsonic inflow conditions, there is a possibility of the occurrence of internal flow choking in dual-thrust motors due to the formation of a fluid throat at the beginning of the transition region induced by area blockage caused by boundary-layer-displacement thickness. It has been observed that a 55% increase in the upstream port area of the dual-thrust motor contributes to a 25% reduction in blockage factor at the transition region, which could negate the internal How choking and supplement with an early choking of the dual-thrust motor nozzle. If the height of the upstream port relative to the motor length is too small, the developing boundary layers from either side of the port can interact, leading to a choked,flow. On the other hand, if the developing boundary layers are far enough apart, then choking does not occur. The blockage factor is greater in magnitude for the choked case than for the unchoked case. More tangible explanations are presented in this paper for the boundary-layer blockage and the internal flow choking in dual-thrust motors, which hitherto has been unexplored.

Turbulence-induced melting of a nonequilibrium vortex crystal in a forced thin fluid film

To gain a better understanding of recent experiments on the turbulence-induced melting of a periodic array of vortices in a thin fluid film, we perform a direct numerical simulation of the two-dimensional Navier-Stokes equations forced such that, at low Reynolds numbers, the steady state of the film is a square lattice of vortices. We find that as we increase the Reynolds number, this lattice undergoes a series of nonequilibrium phase transitions, first to a crystal with a different reciprocal lattice and then to a sequence of crystals that oscillate in time. Initially, the temporal oscillations are periodic; this periodic behaviour becoming more and more complicated with increasing Reynolds number until the film enters a spatially disordered nonequilibrium statistical steady state that is turbulent. We study this sequence of transitions using fluid-dynamics measures, such as the Okubo-Weiss parameter that distinguishes between vortical and extensional regions in the flow, ideas from nonlinear dynamics, e.g. Poincare maps, and theoretical methods that have been developed to study the melting of an equilibrium crystal or the freezing of a liquid and that lead to a natural set of order parameters for the crystalline phases and spatial autocorrelation functions that characterize short- and long-range order in the turbulent and crystalline phases, respectively.

Large-scale structures in turbulent and reverse-transitional sink flow boundary layers

Aspects of large-scale organized structures in sink flow turbulent and reverse-transitional boundary layers are studied experimentally using hot-wire anemometry. Each of the present sink flow boundary layers is in a state of 'perfect equilibrium' or 'exact self-preservation' in the sense of Townsend (The Structure of Turbulent Shear Flow, 1st and 2nd edns, 1956, 1976, Cambridge University Press) and Rotta (Progr. Aeronaut. Sci., vol. 2, 1962, pp. 1-220) and conforms to the notion of 'pure wall-flow' (Coles, J. Aerosp. Sci., vol. 24, 1957, pp. 495-506), at least for the turbulent cases. It is found that the characteristic inclination angle of the structure undergoes a systematic decrease with the increase in strength of the streamwise favourable pressure gradient. Detectable wall-normal extent of the structure is found to be typically half of the boundary layer thickness. Streamwise extent of the structure shows marked increase as the favourable pressure gradient is made progressively severe. Proposals for the typical eddy forms in sink flow turbulent and reverse-transitional flows are presented, and the possibility of structural self-organization (i.e. individual hairpin vortices forming streamwise coherent hairpin packets) in these flows is also discussed. It is further indicated that these structural ideas may be used to explain, from a structural viewpoint, the phenomenon of soft relaminarization or reverse transition of turbulent boundary layers when subjected to strong streamwise favourable pressure gradients. Taylor's 'frozen turbulence' hypothesis is experimentally shown to be valid for flows in the present study even though large streamwise accelerations are involved, the flow being even reverse transitional in some cases. Possible conditions, which are required to be satisfied for the safe use of Taylor's hypothesis in pressure-gradient-driven flows, are also outlined. Measured convection velocities are found to be fairly close to the local mean velocities (typically 90% or more) suggesting that the structure gets convected downstream almost along with the mean flow.

Flow Generated By Pitched Blade Turbines Ii: Simulation Using Κ-Ε Model

Experimental data on average velocity and turbulence intensity generated by pitched blade downflow turbines (PTD) were presented in Part I of this paper. Part II presents the results of the simulation of flow generated by PTD The standard κ-ε model along with the boundary conditions developed in the Part 1 have been employed to predict the flow generated by PTD in cylindrical baffled vessel. This part describes the new software FIAT (Flow In Agitated Tanks) for the prediction of three dimensional flow in stirred tanks. The basis of this software has been described adequately. The influence of grid size, impeller boundary conditions and values of model parameters on the predicted flow have been analysed. The model predictions successfully reproduce the three dimensionality and the other essential characteristics of the flow. The model can be used to improve the overall understanding about the relative distribution of turbulence by PTD in the agitated tank

A study on boundary-layer transition induced by free-stream turbulence

Boundary-layer transition at different free-stream turbulence levels has been investigated using the particle-image velocimetry technique. The measurements show organized positive and negative fluctuations of the streamwise fluctuating velocity component, which resemble the forward and backward jet-like structures reported in the direct numerical simulation of bypass transition. These fluctuations are associated with unsteady streaky structures. Large inclined high shear-layer regions are also observed and the organized negative fluctuations are found to appear consistently with these inclined shear layers, along with highly inflectional instantaneous streamwise velocity profiles. These inflectional velocity profiles are similar to those in the ribbon-induced boundary-layer transition. An oscillating-inclined shear layer appears to be the turbulent spot-precursor. The measurements also enabled to compare the actual turbulent spot in bypass transition with the simulated one. A proper orthogonal decomposition analysis of the fluctuating velocity field is carried out. The dominant flow structures of the organized positive and negative fluctuations are captured by the first few eigenfunction modes carrying most of the fluctuating energy. The similarity in the dominant eigenfunctions at different Reynolds numbers suggests that the flow prevails its structural identity even in intermittent flows. This analysis also indicates the possibility of the existence of a spatio-temporal symmetry associated with a travelling wave in the flow.

A numerical study of division of flow in open channels

A two-dimensional numerical model which employs the depth-averaged forms of continuity and momentum equations along with k-e turbulence closure scheme is used to simulate the flow at the open channel divisions. The model is generalised to flows of arbitrary geometries and MacCormack finite volume method is used for solving governing equations. Application of cartesian version of the model to analyse the flow at right-angled junction is presented. The numerical predictions are compared with experimental data of earlier investigators and measurements made as part of the present study. Performance of the model in predicting discharge distribution, surface profiles, separation zone parameters and energy losses is evaluated and discussed in detail. To illustrate the application of the numerical model to analyse the flow in acute angled offtakes and streamlined branch entries, a few computational results are presented.

Turbulence measurements in an ‘equilibrium’ axisymmetric wall jet Journal

This paper reports measurements of turbulent quantities in an axisymmetric wall jet subjected to an adverse pressure gradient in a conical diffuser, in such a way that a suitably defined pressure-gradient parameter is everywhere small. Self-similarity is observed in the mean velocity profile, as well as the profiles of many turbulent quantities at sufficiently large distances from the injection slot. Autocorrelation measurements indicate that, in the region of turbulent production, the time scale of ν fluctuations is very much smaller than the time scale of u fluctuations. Based on the data on these time scales, a possible model is proposed for the Reynolds stress. One-dimensional energy spectra are obtained for the u, v and w components at several points in the wall jet. It is found that self-similarity is exhibited by the one-dimensional wavenumber spectrum of $\overline{q^2}(=\overline{u^2}+\overline{v^2}+\overline{w^2})$, if the half-width of the wall jet and the local mean velocity are used for forming the non-dimensional wavenumber. Both the autocorrelation curves and the spectra indicate the existence of periodicity in the flow. The rate of dissipation of turbulent energy is estimated from the $\overline{q^2}$ spectra, using a slightly modified version of a previously suggested method.