Relaminarization in highly accelerated turbulent boundary layers

The mean flow development in an initially turbulent boundary layer subjected to a large favourable pressure gradient beginning at a point x0 is examined through analyses expected a priori to be valid on either side of relaminarization. The ‘quasi-laminar’ flow in the later stages of reversion, where the Reynolds stresses have by definition no significant effect on the mean flow, is described by an asymptotic theory constructed for large values of a pressure-gradient parameter Λ, scaled on a characteristic Reynolds stress gradient. The limiting flow consists of an inner laminar boundary layer and a matching inviscid (but rotational) outer layer. There is consequently no entrainment to lowest order in Λ−1, and the boundary layer thins down to conserve outer vorticity. In fact, the predictions of the theory for the common measures of boundary-layer thickness are in excellent agreement with experimental results, almost all the way from x0. On the other hand the development of wall parameters like the skin friction suggests the presence of a short bubble-shaped reverse-transitional region on the wall, where neither turbulent nor quasi-laminar calculations are valid. The random velocity fluctuations inherited from the original turbulence decay with distance, in the inner layer, according to inverse-power laws characteristic of quasi-steady perturbations on a laminar flow. In the outer layer, there is evidence that the dominant physical mechanism is a rapid distortion of the turbulence, with viscous and inertia forces playing a secondary role. All the observations available suggest that final retransition to turbulence quickly follows the onset of instability in the inner layer.It is concluded that reversion in highly accelerated flows is essentially due to domination of pressure forces over the slowly responding Reynolds stresses in an originally turbulent flow, accompanied by the generation of a new laminar boundary layer stabilized by the favourable pressure gradient.

Pulsatile flow in tubes of varying cross-sections

The pulsatile flow of an incompressible viscous fluid in a cylindrical tube of varying cross section is investigated for small Reynolds numbers. The solutions consist of a stedy and an oscillatory part. The shear stress distribution on the wall is evaluated and discussed in detail for special geometries like tapered tubes, locally constricted tubes and peristaltic tubes. The existence of separation in the flow field is noticed.

Numerical studies on quasilinear and linear elliptic equations

We explore here the acceleration of convergence of iterative methods for the solution of a class of quasilinear and linear algebraic equations. The specific systems are the finite difference form of the Navier-Stokes equations and the energy equation for recirculating flows. The acceleration procedures considered are: the successive over relaxation scheme; several implicit methods; and a second-order procedure. A new implicit method—the alternating direction line iterative method—is proposed in this paper. The method combines the advantages of the line successive over relaxation and alternating direction implicit methods. The various methods are tested for their computational economy and accuracy on a typical recirculating flow situation. The numerical experiments show that the alternating direction line iterative method is the most economical method of solving the Navier-Stokes equations for all Reynolds numbers in the laminar regime. The usual ADI method is shown to be not so attractive for large Reynolds numbers because of the loss of diagonal dominance. This loss can however be restored by a suitable choice of the relaxation parameter, but at the cost of accuracy. The accuracy of the new procedure is comparable to that of the well-tested successive overrelaxation method and to the available results in the literature. The second-order procedure turns out to be the most efficient method for the solution of the linear energy equation.

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.

Mechanics of transition in an axisymmetric laminar boundary layer on a circular cylinder

Measurements of the growth of artificially generated turbulent spots and intermittency distribution in the transition region on a circular cylinder in axial flow show that the instability Reynolds number of 11,000 has a marked effect on the properties. In particular, it is found that the spot production in the initial region when a single turbulent spot has not yet wrapped around the cylinder and the propagation is essentially two-dimensional, is significantly altered. But the transition in the downstream or latter region, where most of the turbulent spots propagate onedimensionally (like the turbulent plugs in a pipe), is not affected. When the radius Reynolds number is more than 11,000, the intermittency law in the initial region is essentially the same as in twodimensional flow on a flat plate and in the latter region it is the one-dimensional flow in a pipe, the demarcation between the two regions being quite sharp.

Zero-crossings in turbulent signals

A primary motivation for this work arises from the contradictory results obtained in some recent measurements of the zero-crossing frequency of turbulent fluctuations in shear flows. A systematic study of the various factors involved in zero-crossing measurements shows that the dynamic range of the signal, the discriminator characteristics, filter frequency and noise contamination have a strong bearing on the results obtained. These effects are analysed, and explicit corrections for noise contamination have been worked out. New measurements of the zero-crossing frequency N0 have been made for the longitudinal velocity fluctuation in boundary layers and a wake, for wall shear stress in a channel, and for temperature derivatives in a heated boundary layer. All these measurements show that a zero-crossing microscale, defined as Λ = (2πN0)−1, is always nearly equal to the well-known Taylor microscale λ (in time). These measurements, as well as a brief analysis, show that even strong departures from Gaussianity do not necessarily yield values appreciably different from unity for the ratio Λ/λ. Further, the variation of N0/N0 max across the boundary layer is found to correlate with the familiar wall and outer coordinates; the outer scaling for N0 max is totally inappropriate, and the inner scaling shows only a weak Reynolds-number dependence. It is also found that the distribution of the interval between successive zero-crossings can be approximated by a combination of a lognormal and an exponential, or (if the shortest intervals are ignored) even of two exponentials, one of which characterizes crossings whose duration is of the order of the wall-variable timescale ν/U2*, while the other characterizes crossings whose duration is of the order of the large-eddy timescale δ/U[infty infinity]. The significance of these results is discussed, and it is particularly argued that the pulse frequency of Rao, Narasimha & Badri Narayanan (1971) is appreciably less than the zero-crossing rate.

Peristaltic motion of a micropolar fluid

Peristaltic motion of a micropolar fluid is studied for small amplitudes of peristalic waves under low Reynolds number analysis. The effect of pressure gradient on the secondary motion reveals many interesting and useful results. The critical value of the pressure gradient ensuing the reversal effect in both velocity field and microrotation is evaluated and discussed.

Rapid distortion of axisymmetric turbulence

A generalization of the isotropic theory of Batchelor & Proudman (1954) is developed to estimate the effect of sudden but arbitrary three-dimensional distortion on homogeneous, initially axisymmetric turbulence. The energy changes due to distortion are expressed in terms of the Fourier coefficients of an expansion in zonal harmonics of the two independent scalar functions that describe the axisymmetric spectral tensor. However, for two special but non-trivial forms of this tensor, which represent possibly the simplest kinds of non-isotropic turbulence and specify the angular distribution but not the wavenumber dependence, the energy ratios have been determined in closed form. The deviation of the ratio from its isotropic value is the product of a factor containing R, the initial value of the ratio of the longitudinal to the transverse energy component, and another factor depending only on the geometry of the distortion. It is found that, in axisymmetric and large two-dimensional contractions, the isotropic theory gives nearly the correct longitudinal energy, but (when R > 1) over-estimates the increase in the transverse energy; the product of the two intensities varies little unless the distortion is very large, thus accounting for the stress-freezing observed in rapidly accelerated shear flows.Comparisons with available experimental data for the spectra and for the energy ratios show reasonable agreement. The different ansatzes predict results in broad qualitative agreement with a simple strategem suggested by Reynolds & Tucker (1975), but the quantitative differences are not always negligible.

The structure of turbulent boundary layers along mildly curved surfaces

This paper describes a detailed study of the structure of turbulence in boundary layers along mildly curved convex and concave surfaces. The surface curvature studied corresponds to δ/Rw = ± 0·01, δ being the boundary-layer thickness and Rw the radius of curvature of the wall, taken as positive for convex and negative for concave curvature. Measurements of turbulent energy balance, autocorrelations, auto- and cross-power spectra, amplitude probability distributions and conditional correlations are reported. It is observed that even mild curvature has very strong effects on the various aspects of the turbulent structure. For example, convex curvature suppresses the diffusion of turbulent energy away from the wall, reduces drastically the integral time scales and shifts the spectral distributions of turbulent energy and Reynolds shear stress towards high wavenumbers. Exactly opposite effects, though generally of a smaller magnitude, are produced by concave wall curvature. It is also found that curvature of either sign affects the v fluctuations more strongly than the u fluctuations and that curvature effects are more significant in the outer region of the boundary layer than in the region close to the wall. The data on the conditional correlations are used to study, in detail, the mechanism of turbulent transport in curved boundary layers. (Published Online April 12 2006)

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.

Response of a turbulent boundary layer to a step change in surface heat flux

Measurements of both the velocity and the temperature field have been made in the thermal layer that grows inside a turbulent boundary layer which is subjected to a small step change in surface heat flux. Upstream of the step, the wall heat flux is zero and the velocity boundary layer is nearly self-preserving. The thermal-layer measurements are discussed in the context of a self-preserving analysis for the temperature disturbance which grows underneath a thick external turbulent boundary layer. A logarithmic mean temperature profile is established downstream of the step but the budget for the mean-square temperature fluctuations shows that, in the inner region of the thermal layer, the production and dissipation of temperature fluctuations are not quite equal at the furthest downstream measurement station. The measurements for both the mean and the fluctuating temperature field indicate that the relaxation distance for the thermal layer is quite large, of the order of 1000θ0, where θ0 is the momentum thickness of the boundary layer at the step. Statistics of the thermal-layer interface and conditionally sampled measurements with respect to this interface are presented. Measurements of the temperature intermittency factor indicate that the interface is normally distributed with respect to its mean position. Near the step, the passive heat contaminant acts as an effective marker of the organized turbulence structure that has been observed in the wall region of a boundary layer. Accordingly, conditional averages of Reynolds stresses and heat fluxes measured in the heated part of the flow are considerably larger than the conventional averages when the temperature intermittency factor is small.

Simulation of laminar flow in a three-dimensional lid-driven cavity by lattice Boltzmann method

Purpose - The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid-driven flow in a three-dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two-dimensional (2D) cavity. Design/methodology/approach - The second-order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25-4) and Reynolds numbers (0.01-1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier-Stokes (NS) equations and with available experimental data. Findings - The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side-walls in transverse direction (in 3D), due to which the vorticity in the core-region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry-plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid-motion are observed; these are weak at low. Reynolds numbers, but become quite strong at high Reynolds numbers. Originality/value - The findings will be useful in the study of variety of enclosed fluid flows.

Numerical Study of Heat Transfer From Pin-Fin Heat Sink Using Steady and Pulsated Impinging Jets

This paper investigates numerically the heat transfer characteristics of confined slot jet impingement on a pin-fin heat sink. A variety of pin-fin heat sinks is investigated, and the resulting enhancement of heat transfer studied. The distribution of heat transfer coefficient on the top surface of the base plate and that along the fin height are examined. Both steady and pulsated jets are studied. It is observed that for a steady jet impingement on a pin-fin heat sink, the effective heat transfer coefficient increases with fin height, leading to a corresponding decrease in base plate temperature for the same heat flux. In the case of pulsated jets, the influence of pulse frequency and the Reynolds number is examined, and their effect on the effective heat transfer coefficient is studied.

Liquid Sheet Breakup in Gas-Centered Swirl Coaxial Atomizers

The study deals with the breakup behavior of swirling liquid sheets discharging from gas-centered swirl coaxial atomizers with attention focused toward the understanding of the role of central gas jet on the liquid sheet breakup. Cold flow experiments on the liquid sheet breakup were carried out by employing custom fabricated gas-centered swirl coaxial atomizers using water and air as experimental fluids. Photographic techniques were employed to capture the flow behavior of liquid sheets at different flow conditions. Quantitative variation on the breakup length of the liquid sheet and spray width were obtained from the measurements deduced from the images of liquid sheets. The sheet breakup process is significantly influenced by the central air jet. It is observed that low inertia liquid sheets are more vulnerable to the presence of the central air jet and develop shorter breakup lengths at smaller values of the air jet Reynolds number Re-g. High inertia liquid sheets ignore the presence of the central air jet at smaller values of Re-g and eventually develop shorter breakup lengths at higher values of Re-g. The experimental evidences suggest that the central air jet causes corrugations on the liquid sheet surface, which may be promoting the production of thick liquid ligaments from the sheet surface. The level of surface corrugations on the liquid sheet increases with increasing Re-g. Qualitative analysis of experimental observations reveals that the entrainment process of air established between the inner surface of the liquid sheet and the central air jet is the primary trigger for the sheet breakup.

Metamodeling approach to predict friction factor of alluvial channel

The flow resistance of an alluvial channel flow is not only affected by the Reynolds number and the roughness conditions but also the Froude number. Froude number is the most basic parameter in the case of the alluvial channel, thus effect of Froude number on resistance to flow should be considered in the formulation of the friction factor, which is not in the case of present available resistance equations. At present, no generally acceptable quantitative description of the effects of the Froude number on hydraulic resistance has been developed. Metamodeling technique, which is particularly useful in modeling a complex processes or where knowledge of the physics is limited, is presented as a tool complimentary to modeling friction factor in alluvial channels. Present work uses, a radial basis metamodel, which is a type of neural network modeling, to find the effect of Froude number on the flow resistance. Based on the experimental data taken from different sources, it has been found that the predicting capability of the present model is on acceptable level. Present work also tries in formulating an empirical equation for resistance in alluvial channel comprising all the three majorm, parameters, namely, roughness parameter, Froude number and Reynolds number. (C) 2009 Elsevier B.V. All rights reserved.