CMC simulations of lifted turbulent jet flame in a vitiated coflow

Lifted turbulent jet diffusion flame is simulated using Conditional Moment Closure (CMC). Specifically, the burner configuration of Cabra et al. [R. Cabra, T. Myhrvold, J.Y. Chen. R.W. Dibble, A.N. Karpetis, R.S. Barlow, Proc. Combust. Inst. 29 (2002) 1881-1887] is chosen to investigate H-2/N-2 jet flame supported by a vitiated coflow of products of lean H-2/air combustion. A 2D, axisymmetric flow-model fully coupled with the scalar fields, is employed. A detailed chemical kinetic scheme is included, and first order CIVIC is applied. Simulations are carried out for different jet velocities and coflow temperatures (T-c) The predicted liftoff generally agrees with experimental data, as well as joint-PDF results. Profiles of mean scalar fluxes in the mixture fraction space, for T-c = 1025 and 1080 K reveal that (1) Inside the flame zone, the chemical term balances the molecular diffusion term, and hence the Structure is of a diffusion flamelet for both cases. (2) In the pre-flame zone, the structure depends on the coflow temperature: for the 1025 K case, the chemical term being small, the advective term balances the axial turbulent diffusion term. However, for the 1080 K case. the chemical term is large and balances the advective term, the axial turbulent diffusion term being small. It is concluded that, lift-off is controlled (a) by turbulent premixed flame propagation for low coflow temperature while (b) by autoignition for high coflow temperature. (C) 2009 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Turbulent mixed convection in a shallow enclosure with a series of heat generating components

Turbulent mixed convection flow and heat transfer in a shallow enclosure with and without partitions and with a series of block-like heat generating components is studied numerically for a range of Reynolds and Grashof numbers with a time-dependent formulation. The flow and temperature distributions are taken to be two-dimensional. Regions with the same velocity and temperature distributions can be identified assuming repeated placement of the blocks and fluid entry and exit openings at regular distances, neglecting the end wall effects. One half of such module is chosen as the computational domain taking into account the symmetry about the vertical centreline. The mixed convection inlet velocity is treated as the sum of forced and natural convection components, with the individual components delineated based on pressure drop across the enclosure. The Reynolds number is based on forced convection velocity. Turbulence computations are performed using the standard k– model and the Launder–Sharma low-Reynolds number k– model. The results show that higher Reynolds numbers tend to create a recirculation region of increasing strength in the core region and that the effect of buoyancy becomes insignificant beyond a Reynolds number of typically 5×105. The Euler number in turbulent flows is higher by about 30 per cent than that in the laminar regime. The dimensionless inlet velocity in pure natural convection varies as Gr1/3. Results are also presented for a number of quantities of interest such as the flow and temperature distributions, Nusselt number, pressure drop and the maximum dimensionless temperature in the block, along with correlations.

Unsteady Incompressible Two-Dimensional and Axisymmetric Turbulent Boundary Layer Flows

The unsteady turbulent incompressible boundary-layer flow over two-dimensional and axisymmetric bodies with pressure gradient has been studied. An eddy-viscosity model has been used to model the Reynolds shear stress. The unsteadiness is due to variations in the free stream velocity with time. The nonlinear partial differential equation with three independent variables governing the flow has been solved using Keller's Box method. The results indicate that the free stram velocity distribution exerts strong influence on the boundary-layer characteristics. The point of zero skin friction is found to move upstream as time increases.

Wall Shear Fluctuations in a Turbulent Boundary Layer

IN the last two decades, the instantaneous structure of a turbulent boundary layer has been examined by many in an effort to understand the dynamics of the flow. Distinct and well-defined flow patterns that seem to have great relevance to the turbulence production mechanism have been observed in the wall region.1'2 The flow near the wall is intermittent with periodic eruptions of the fluid, a phenomenon generally termed "bursting process." Earlier investigations in this field were limited to liquid flows at low speeds and the entire flowpattern was observed using flow visualization techniques.Study was later extended to boundary-layer flows in windtunnels at higher speeds and Reynolds numbers using hot-wiresignals for the analysis of the bursting phenomenon.

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.

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

Evaluation of k-epsilon and RNG turbulence models for confined swirling and non-swirling flow

In a very recent study [1] the Renormalisation Group (RNG) turbulence model was used to obtain flow predictions in a strongly swirling quarl burner, and was found to perform well in predicting certain features that are not well captured using less sophisticated models of turbulence. The implication is that the RNG approach should provide an economical and reliable tool for the prediction of swirling flows in combustor and furnace geometries commonly encountered in technological applications. To test this hypothesis the present work considers flow in a model furnace for which experimental data is available [2]. The essential features of the flow which differentiate it from the previous study [1] are that the annular air jet entry is relatively narrow and the base wall of the cylindrical furnace is at 90 degrees to the inlet pipe. For swirl numbers of order 1 the resulting flow is highly complex with significant inner and outer recirculation regions. The RNG and standard k-epsilon models are used to model the flow for both swirling and non-swirling entry jets and the results compared with experimental data [2]. Near wall viscous effects are accounted for in both models via the standard wall function formulation [3]. For the RNG model, additional computations with grid placement extending well inside the near wall viscous-affected sublayer are performed in order to assess the low Reynolds number capabilities of the model.

Computation of confined swirling flow using RNG and k-epsilon turbulence models

In this work we numerically model isothermal turbulent swirling flow in a cylindrical burner. Three versions of the RNG k-epsilon model are assessed against performance of the standard k-epsilon model. Sensitivity of numerical predictions to grid refinement, differing convective differencing schemes and choice of (unknown) inlet dissipation rate, were closely scrutinised to ensure accuracy. Particular attention is paid to modelling the inlet conditions to within the range of uncertainty of the experimental data, as model predictions proved to be significantly sensitive to relatively small changes in upstream flow conditions. We also examine the characteristics of the swirl--induced recirculation zone predicted by the models over an extended range of inlet conditions. Our main findings are: - (i) the standard k-epsilon model performed best compared with experiment; - (ii) no one inlet specification can simultaneously optimize the performance of the models considered; - (iii) the RNG models predict both single-cell and double-cell IRZ characteristics, the latter both with and without additional internal stagnation points. The first finding indicates that the examined RNG modifications to the standard k-e model do not result in an improved eddy viscosity based model for the prediction of swirl flows. The second finding suggests that tuning established models for optimal performance in swirl flows a priori is not straightforward. The third finding indicates that the RNG based models exhibit a greater variety of structural behaviour, despite being of the same level of complexity as the standard k-e model. The plausibility of the predicted IRZ features are discussed in terms of known vortex breakdown phenomena.

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.

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.

Stability and transition in the flow of polymer solutions

The linear stability analysis of a plane Couette flow of viscoelastic fluid have been studied with the emphasis on two dimensional disturbances with wave number k similar to Re-1/2, where Re is Reynolds number based on maximum velocity and channel width. We employ three models to represent the dilute polymer solution: the classical Oldroyd-B model, the Oldroyd-B model with artificial diffusivity and the non-homogeneous polymer model. The result of the linear stability analysis is found to be sensitive to the polymer model used. While the plane Couette flow is found to be stable to infinitesimal disturbances for the first two models, the last one exhibits a linear instability.

Effect of surfactants on drop breakage in turbulent liquid dispersions

The existingm odels of drop breakage in stirred dispersions grossly overpredict the maximum drop size when surface active agents are present inspite of using the lowered value of interfacial tension. It is shown that the difference in the values of dynamic and static interfacial tension, aids the turbulent stresses in drop breakage. When the difference is zero, e.g. for pure liquids and for high concentration of surfactants, the influence of the addition of surfactant is merely to reduce the interfacial tension and can be accounted for by existingm odels. A modified model has been developed, where the drop breakage is assumed to be represented by a Voigt element. The deforming stresses are due to turbulence and the difference between dynamic and static interfacial tensions. The resisting stresses arise due to interfacial tension and the viscous flow inside the drop. The model yields the existing expressions for dmax as special cases. The model has been found to be satisfactory when tested against experimental results using the styrene-water-teepol system.

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.

Turbulent near-wake development behind a flat plate

An asymptotic analysis of the two-dimensional turbulent near-wake flow behind a Rat plate with sharp trailing edge has been formulated, The feature that the near-wake, which is dominated by the mixing of the oncoming turbulent boundary layers retains, to a large extent, the memory of the turbulent structure of the upstream boundary layer has been exploited to develop the analysis. This analysis leads to two regions of the near-wake flow (the inner near-wake and the outer near-wake) for which the governing equations are derived. The matching conditions among these regions lead to a logarithmic variation in the normal direction in the overlapping region surrounding the inner near-wake. These features are validated by the available experimental data. Similarity solutions for the velocity distribution (which satisfy the required matching conditions) in the inner near-wake and outer near-wake regions have been obtained by making the appropriate eddy-viscosity assumptions, Uniformly valid solutions for velocity distribution have been constructed for the near-wake. The solutions show good agreement with available experimental data. (C) Elsevier, Paris.

Constant pressure laminar, transitional and turbulent flows - An approximate unified treatment

A nondimensional number that is constant in two-dimensional, incompressible and constant pressure laminar and fully turbulent boundary, layer flows has been proposed. An extension of this to constant pressure transitional flow is discussed.