Turbulent spot growth in favorable pressure gradients

Experimental results are presented on the lateral growth of turbulent spots in a series of flows with favorable pressure gradients. It is shown that the wedge angle increases slowly with the Reynolds number and that a favorable pressure gradient inhibits the growth of turbulent spots and, in general, results in a nonlinear turbulent wedge. As soon as the pressure gradient decreases to the point where the flow becomes supercritical, however, spot growth increases rapidly and the associated turbulent wedge becomes linear.

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.

Laminar and turbulent natural convection in rectangular cavities

A computer program has been developed for the prediction of buoyancy-driven laminar and turbulent flow in rectangular air-filled two-dimensional cavities with differentially heated side walls. Laminar flow predictions for a square cavity and Rayleigh numbers from Ra = 10^3 up to the onset of unsteady flow have been obtained. Accurate solutions for Ra = 5 x 10^6, 10^7, 5 x 10^7 and 10^8 are presented and an estimate for the critical Rayleigh number at which the steady laminar flow becomes unsteady is given for this geometry. Numerical predictions of turbulent flow have been obtained for RaH~0(10^9 -10^11 ) and compared with existing experimental data. A previously developed second moment closure model (Behnia et al. 1987) has been used to model the turbulence. Results indicate that a second moment closure model is capable of predicting the observed flow features.

Laminar and turbulent natural convection in rectangular cavities

"This chapter discusses laminar and turbulent natural convection in rectangular cavities. Natural convection in rectangular two-dimensional cavities has become a standard problem in numerical heat transfer because of its relevance in understanding a number of problems in engineering. Current research identified a number of difficulties with regard to the numerical methods and the turbulence modeling for this class of flows. Obtaining numerical predictions at high Rayleigh numbers proved computationally expensive such that results beyond Ra ∼ 1014 are rarely reported. The chapter discusses a study in which it was found that turbulent computations in square cavities can't be extended beyond Ra ∼ O (1012) despite having developed a code that proved very efficient for the high Ra laminar regime. As the Rayleigh number increased, thin boundary layers began to form next to the vertical walls, and the central region became progressively more stagnant and highly stratified. Results obtained for the high Ra laminar regime were in good agreement with existing studies. Turbulence computations, although of a preliminary nature, indicated that a second moment closure model was capable of predicting the experimentally observed flow features."--Publisher Summary

Computation of turbulent flow in an IFRF quarl burner

A computational model for isothermal axisymmetric turbulent flow in a quarl burner is set up using the CFD package FLUENT, and numerical solutions obtained from the model are compared with available experimental data. A standard k-e model and and two versions of the RNG k-e model are used to model the turbulence. One of the aims of the computational study is to investigate whether the RNG based k-e turbulence models are capable of yielding improved flow predictions compared with the standard k-e turbulence model. A difficulty is that the flow considered here features a confined vortex breakdown which can be highly sensitive to flow behaviour both upstream and downstream of the breakdown zone. Nevertheless, the relatively simple confining geometry allows us to undertake a systematic study so that both grid-independent and domain-independent results can be reported. The systematic study includes a detailed investigation of the effects of upstream and downstream conditions on the predictions, in addition to grid refinement and other tests to ensure that numerical error is not significant. Another important aim is to determine to what extent the turbulence model predictions can provide us with new insights into the physics of confined vortex breakdown flows. To this end, the computations are discussed in detail with reference to known vortex breakdown phenomena and existing theories. A major conclusion is that one of the RNG k-e models investigated here is able to correctly capture the complex forward flow region inside the recirculating breakdown zone. This apparently pathological result is in stark contrast to the findings of previous studies, most of which have concluded that either algebraic or differential Reynolds stress modelling is needed to correctly predict the observed flow features. Arguments are given as to why an isotropic eddy-viscosity turbulence model may well be able to capture the complex flow structure within the recirculating zone for this flow setup. With regard to the flow physics, a major finding is that the results obtained here are more consistent with the view that confined vortex breakdown is a type of axisymmetric boundary layer separation, rather than a manifestation of a subcritical flow state.

Analysis of turbulent confined highly swirling flow

We report here a CFD model of highly swirling flow in a quarl burner using three versions of the k-epsilon model. Results for the recirculating zone, the bounding shear layer and the downstream flow are presented. We discuss, with suitable qualifications, how the model predictions can inform our understanding of this class of flows.

Determination of skin friction in strong pressure-gradient equilibrium and near-equilibrium turbulent boundary layers

The conventional Clauser-chart method for determination of local skin friction in zero or weak pressure-gradient turbulent boundary layer flows fails entirely in strong pressure-gradient situations. This failure occurs due to the large departure of the mean velocity profile from the universal logarithmic law upon which the conventional Clauser-chart method is based. It is possible to extend this method,even for strong pressure-gradient situations involving equilibrium or near-equilibrium turbulent boundary layers by making use of the so-called non-universal logarithmic laws. These non-universal log laws depend on the local strength of the pressure gradient and may be regarded as perturbations of the universal log law.The present paper shows that the modified Clauser-chart method, so developed, yields quit satisfactory results in terms of estimation of local skin friction in strongly accelerated or retarded equilibrium and near-equilibrium turbulent boundary layers that are not very close to relaminarization or separation.

Transition of turbulent pipe flow

The velocity profile in turbulent pipe flow is usually divided into two regions, a wall or inner region and a core or outer region. For the inner region, the viscosity and wall shear stress are the important parameters governing the velocity distribution whereas for the outer region, the wall reduces the velocity below the maximum velocity independent of viscosity. In the present work, a velocity model is proposed for turbulent flow in the wall region of a pipe covering the entire transition from smooth to rough flows. Coupling this model for the wall region with the power law velocity model for the core region, an equation for the friction factor is obtained. The model constants are evaluated by using Nikuradse's experiments in the fully smooth and rough turbulent flows. The model shows good agreement with the friction factor and the velocity profiles obtained by Nikuradse for the transition region of turbulent flow.

Parametric Analysis Of Turbulent Wall Jets

The concept of a fully-developed flow based on the hypothesis of selective memory is here applied to general wall-jet type flows. In the presence of a (constant) external stream, the free-stream velocity and the jet momentum flux are taken to be the chief quantities governing the development of the wall jet: two additional nondimensional parameters, representing a momentum flux Reynolds number and the relative momentum defect in the initial boundary layer, are shown to have only a secondary effect on the fully-developed flow. The standard correlations so determined are also found to predict quite well the flow development in Gartshore and Newman's experiments on wall jets in adverse pressure gradients; possible reasons for this somewhat surprising result are discussed. Finally it is shown, by application to the still-air case, that the parameters discovered in incompressible flow are, with appropriate but straightforward modification, successful in describing compressible wall jets also.

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.

Turbulent flow computations on a hybrid cartesian point distribution using meshless solver LSFD-U

This paper may be considered as a sequel to one of our earlier works pertaining to the development of an upwind algorithm for meshless solvers. While the earlier work dealt with the development of an inviscid solution procedure, the present work focuses on its extension to viscous flows. A robust viscous discretization strategy is chosen based on positivity of a discrete Laplacian. This work projects meshless solver as a viable cartesian grid methodology. The point distribution required for the meshless solver is obtained from a hybrid cartesian gridding strategy. Particularly considering the importance of an hybrid cartesian mesh for RANS computations, the difficulties encountered in a conventional least squares based discretization strategy are highlighted. In this context, importance of discretization strategies which exploit the local structure in the grid is presented, along with a suitable point sorting strategy. Of particular interest is the proposed discretization strategies (both inviscid and viscous) within the structured grid block; a rotated update for the inviscid part and a Green-Gauss procedure based positive update for the viscous part. Both these procedures conveniently avoid the ill-conditioning associated with a conventional least squares procedure in the critical region of structured grid block. The robustness and accuracy of such a strategy is demonstrated on a number of standard test cases including a case of a multi-element airfoil. The computational efficiency of the proposed meshless solver is also demonstrated. (C) 2010 Elsevier Ltd. All rights reserved.

A possible origin of viscosity in Keplerian accretion disks due to secondary perturbation: Turbulent transport without magnetic fields

The origin of hydrodynamic turbulence in rotating shear flow is a long standing puzzle. Resolving it is especially important in astrophysics when the flow's angular momentum profile is Keplerian which forms an accretion disk having negligible molecular viscosity. Hence, any viscosity in such systems must be due to turbulence, arguably governed by magnetorotational instability, especially when temperature T greater than or similar to 10(5). However, such disks around quiescent cataclysmic variables, protoplanetary and star-forming disks, and the outer regions of disks in active galactic nuclei are practically neutral in charge because of their low temperature, and thus are not expected to be coupled with magnetic fields enough to generate any transport due to the magnetorotational instability. This flow is similar to plane Couette flow including the Coriolis force, at least locally. What drives their turbulence and then transport, when such flows do not exhibit any unstable mode under linear hydrodynamic perturbation? We demonstrate that the three-dimensional secondary disturbance to the primarily perturbed flow that triggers elliptical instability may generate significant turbulent viscosity in the range 0.0001 less than or similar to nu(t) less than or similar to 0.1, which can explain transport in accretion flows.

Studies on Cavitation Inception Process in Separated Flows

Studies related to cavitation inception process in separated flows are reported. Experimental observations of bubble appearance in grooves with laminar or turbulent boundary layer over them have clearly shown that gaseous diffusion process is significantly enhanced in turbulent flow. This process can lead to local nuclei size modification in environment similar to that of flow over a groove, like laminar separation "bubbles." Cavitation inception modeling including this aspect is carried out for predicting inception conditions associated with "bubble-ring" cavitation commonly observed on hemispherically nosed axisymmetric body. Qualitative dependence of predicted inception numbers with velocity is found to agree very well with experimental observations of Carroll (1981).

An experimental study of reverse transition in two-dimensional channel flow

An experimental investigation on reverse transition from turbulent to laminar flow in a two-dimensional channel was carried out. The reverse transition occurred when Reynolds number of an initially turbulent flow was reduced below a certain value by widening the duct in the lateral direction. The experiments were conducted at Reynolds numbers of 625, 865, 980 and 1250 based on half the height of the channel and the average of the mean velocity. At all these Reynolds numbers the initially turbulent mean velocity profiles tend to become parabolic. The longitudinal and vertical velocity fluctuations ($\overline{u^{\prime 2}}$ and $\overline{v^{\prime 2}}$) averaged over the height of the channel decrease exponentially with distance downstream, but $\overline{u^{\prime}v^{\prime}} $ tends to become zero at a reasonably well-defined point. During reverse transition $\overline{u^{\prime}}\overline{v^{\prime}}/\sqrt{\overline{u^{\prime 2}}}\sqrt{\overline{v^{\prime 2}}}$ also decreases as the flow moves downstream and Lissajous figures taken with u’ and v’ signals confirm this trend. There is approximate similarly between $\overline{u^{\prime 2}} $ profiles if the value of $\overline{u^{\prime 2}_{\max}} $ and the distance from the wall at which it occurs are taken as the reference scales. The spectrum of $\overline{u^{\prime 2}} $ is almost similar at all stations and the non-dimensional spectrum is exponential in wave-number. All the turbulent quantities, when plotted in appropriate co-ordinates, indicate that there is a definite critical Reynolds number of 1400±50 for reverse transition.

Two-dimensional moist stratified turbulence and the emergence of vertically sheared horizontal flows

Moist stratified turbulence is studied in a two-dimensional Boussinesq system influenced by condensation and evaporation. The problem is set in a periodic domain and employs simple evaporation and condensation schemes, wherein both the processes push parcels towards saturation. Numerical simulations demonstrate the emergence of a moist turbulent state consisting of ordered structures with a clear power-law type spectral scaling from initially spatially uncorrelated conditions. An asymptotic analysis in the limit of rapid condensation and strong stratification shows that, for initial conditions with enough water substance to saturate the domain, the equations support a straightforward state of moist balance characterized by a hydrostatic, saturated, vertically sheared horizontal flow (VSHF). For such initial conditions, by means of long time numerical simulations, the emergence of moist balance is verified. Specifically, starting from uncorrelated data, subsequent to the development of a moist turbulent state, the system experiences a rather abrupt transition to a regime which is close to saturation and dominated by a strong VSHF. On the other hand, initial conditions which do not have enough water substance to saturate the domain, do not attain moist balance. Rather, the system is observed to remain in a turbulent state and oscillates about moist balance. Even though balance is not achieved with these general initial conditions, the time scale of oscillation about moist balance is much larger than the imposed time scale of condensation and evaporation, thus indicating a distinct dominant slow component in the moist stratified two-dimensional turbulent system. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3694805]