NUMERICAL PREDICTION OF SEMI-CONFINED JET IMPINGEMENT AND COMPARISON WITH EXPERIMENTAL DATA

The standard k–ε eddy viscosity model of turbulence in conjunction with the logarithmic law of the wall has been applied to the prediction of a fully developed turbulent axisymmetric jet impinging within a semi-confined space. A single geometry with a Reynolds number of 20,000 and a nozzle-to-plate spacing of two diameters has been considered with inlet boundary conditions based on measured profiles of velocity and turbulence. Velocity, turbulence and heat transfer data have been obtained using laser–Doppler anemometry and liquid crystal thermography respectively. In the developing wall jet, numerical results of heat transfer compare to within 20% of experiment where isotropy prevails and the trends in turbulent kinetic energy are predicted. However, stagnation point heat transfer is overpredicted by about 300%, which is attributed directly to the turbulence model and inapplicability of the wall function.     

Computations for a jet impinging obliquely on a flat surface

A SIMPLE-C algorithm and Jones-Launder k-ε two-equation turbulence model are used to simulate a two-dimensional jet impinging obliquely on a flat surface. Both the continuity and momentum equations for the unsteady state are cast into suitable finite difference equations. The pressure, velocity, turbulent kinetic energy and turbulent energy dissipation rate distributions are solved and show good agreement with various experimental data. The calculations show that the flow field structure of the jet impinging obliquely on a flat surface is strongly affected by the oblique impingement angle. The maximum pressure zone of the obliquely impinging jet flow field moves towards the left as the oblique impingement angle is decreased.     

On the analysis of an impinging jet on ground effects

Laser Doppler measurements and flow visualization are presented for a turbulent circular jet emerging into a low-velocity cross-stream and, then, impinging on a flat surface perpendicular to the jet-nozzle axis. The experiments were performed for a Reynolds number based on the jet-exit conditions of 6 × 10⁴, a jet-to-crossflow velocity ratio of 30 and for an impingement height of 5 jet diameters and include mean and turbulent velocity characteristics along the two normal directions contained in planes parallel to the nozzle axis. The results, which have relevance to flows found beneath VSTOL aircraft in ground effect, show the presence of a complex 3-D scarf vortex formed around the impinging jet. In zones where measurement data are not available, the flow details are numerically-visualized using a solution of the finite difference form of the fully threedimensional Reynolds-averaged Navier-Stokes equations, incorporating the turbulence viscosity concept. The turbulent structure of the flow is affected by flow distortion at the impinging zone, which results in an unconventional behaviour of the dimensionless structure parameters that determine the empirical constants in engineering models of turbulence. The relative magnitude of the terms involved in the transport equations for the turbulent stresses is quantified from the experimental data in order to assess the importance of these effects and show the extent to which the turbulent structure of the impingement zone is affected by extra rates of strain.     

Unsteady flow and heat transfer analysis of an impinging synthetic jet

The present paper focuses on the analysis of unsteady flow and heat transfer regarding an axisymmetric impinging synthetic jet on a constant heat flux disc. Synthetic jet is a zero net mass flux jet that provides an unsteady flow without any external source of fluid. Present results are validated against the available experimental data showing that the SST/k − ω turbulence model is more accurate and reliable than the standard and low-Re k − ε models for predicting heat transfer from an impinging synthetic jet. It is found that the time-averaged Nusselt number enhances as the nozzle-to-plate distance is increased. As the oscillation frequency in the range of 16–400 Hz is increased, the heat transfer is enhanced. It is shown that the instantaneous Nu distribution along the wall is influenced mainly by the interaction of produced vortex ring and wall boundary layer. Also, the fluctuation level of Nu decreases as the frequency is raised.     

SOME VALIDATION OF STANDARD,MODIFIED AND NON-LINEAR k–ε TURBULENCE MODELS

Standard, modified and non-linear k–ε: turbulence models are validated against three axisymmetric flow problems—flow through a pipe expansion, flow through a pipe constriction and an impinging jet problem—to underpin knowledge about the solution quality obtained from two-equation turbulence models. The extended models improve the prediction of turbulence as a flow approaches a stagnation point and the non-linear model allows for the prediction of anisotropic turbulence. Significantly different values for the non-linear model coefficients are proposed in comparison with values found in the literature. Nevertheless, current turbulence models are still unable to accurately predict the spreading rate of shear layers. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids, 24: 965–986, 1997.     

On the performance of sudden‐expansion pipe without/with cross‐flow injection: experimental and numerical studies

The paper explores the possibilities that different turbulence closures offer, for in‐depth analysis of a complex flow. The case under investigation is steady, turbulent flow in a pipe with sudden expansion without/with normal‐to‐wall injection through jets. This is a typical geometry where generation of turbulence energy takes place, due to sudden change in boundary conditions. This study is aimed at investigating the capability of a developed computational program by the present authors with three different turbulence models to calculate the mean flow variables. Three two‐equation models are implemented, namely the standard linear k ? ε model, the low Reynolds number k ? ε model and the cubic nonlinear eddy viscosity (NLEV) k ? ε model. The performance of the chosen turbulence models is investigated with regard to the available data in the literature including velocity profiles, turbulent kinetic energy and reattachment position in a pipe expansion. In order to further assess the reliability of the turbulence models, an experimental program was conducted by the present authors also at the fluid mechanics laboratory of Menoufiya University. Preliminary measurements, including the surface pressure along the two walls of the expansion pipe and the pressure drop without and with the presence of different arrangements of wall jets produced by symmetrical or asymmetrical fluid cross‐flow injection, are introduced. The results of the present studies demonstrate the superiority of the cubic NLEV k ? ε model in predicting the flow characteristics over the entire domain. The simple low Reynolds number k ? ε model also gives good prediction, especially when the reattachment point is concerned. The evaluation of the reattachment point and the pressure‐loss coefficient is numerically addressed in the paper using the cubic NLEV k ? ε model. The results show that the injection location can control the performance of the pipe‐expansion system. It is concluded that the introduction of flow injection can increase the energy loss in the pipe expansion. The near‐field turbulence structure is also considered in the present study and it is noticed that the turbulence level is strongly affected by the cross‐flow injection and the jet location. Copyright © 2011 John Wiley & Sons, Ltd.     

A numerical study of turbulent flow over a two‐dimensional hill

Calculations of mean velocities and Reynolds stresses are reported for the recirculating flow established in the wake of two‐dimensional polynomial‐shaped obstacles that are symmetrical about a vertical axis and mounted in the water channel downstream of a fully developed channel flow for Re=6×10⁴. The study involves calculations of mean and fluctuating flow properties in the streamwise and spanwise directions and include comparisons with experimental data [Almeida GP, Durão DFG, Heitor MV. Wake flows behind two‐dimensional model hills. Experimental Thermal and Fluid Science 1993; 7: 87–101] for flow around a single obstacle with data resulting from the interaction of consecutive obstacles, using two versions of the low‐Reynolds number differential second‐moment (DSM) closure model. The results include analysis of the turbulent stresses in local flow co‐ordinates and reveal flow structure qualitatively similar to that found in other turbulent flows with a reattachment zone. It is found that the standard isotropization of production model (IPM), based on that proposed by Gibson and Launder [Ground effects on pressure fluctuations in the atmospheric boundary layer. Journal of Fluid Mechanics 1978; 86(3): 191–511], with the incorporation of the wall reflection model of Craft and Launder [New wall‐reflection model applied to the turbulent impinging jet. AIAA Journal 1992; 32(12): 2970–2972] predicts the mean velocities quite well, but underestimates the size of the recirculation region and turbulent quantities in the shear layer. These inadequacies are circumvented by adopting a new cubic Reynolds stress closure scheme based on that more recently developed by Craft and Launder [A Reynolds stress closure designed for complex geometries. International Journal of Heat and Fluid Flow 1996; 17: 245–254] which satisfies the two component limit (TCL) of turbulence. In this model the geometry‐specific quantities, such as the wall‐normal vector or wall distance, are replaced by invariant dimensionless gradient indicators. Also, the model captures the diverse behaviour of the different components of the stress dissipation, ε_ij, near the wall and uses a novel decomposition for the fluctuating pressure terms. Copyright © 2001 John Wiley & Sons, Ltd.     

Control of a submerged jet in a thin rectangular cavity

An innovative method is presented for control of an oscillatory turbulent jet in a thin rectangular cavity with a thickness to width ratio of 0.16. Jet flow control is achieved by mass injection of a secondary jet into the region above the submerged primary jet nozzle exit and perpendicular to the primary nozzle axis. An experimental model, a 2-D and a 3-D computational fluid dynamics (CFD) model are used to investigate the flow characteristics under various secondary injection mass flow rates and injection positions. Two-dimensional laser Doppler anemometry (LDA) measurements are compared with results from the CFD models, which incorporate a standard k–ε turbulence model or a 2-D and 3-D realisable k–ε model. Experimental results show deflection angles up to 23.3° for 24.6% of relative secondary mass flow are possible. The key to high jet control sensitivity is found to be lateral jet momentum with the optimum injection position at 12% of cavity width (31.6% of the primary nozzle length) above the primary nozzle exit. CFD results also show that a standard k–ε turbulence closure with nonequilibrium wall functions provides the best predictions of the flow.     

Numerical computation of an impinging jet with uniform wall suction

The paper presents numerical predictions of a turbulent axisymmetric jet impinging onto a porous plate, based on a finite volume method of solving the Navier-Stokes equations for an incompressible air jet with the K–ε turbulence model. The velocity and pressure terms of the momentum equations are solved by the SIMPLE (semi-implicit method for pressure-linked equation) method. In this study, non-uniform staggered grids are used. The parameters of interest include the nozzle-to-wall distance and the suction velocity. The results of the present calculations are compared with available data reported in the literature. It is found that suction effects reduce the boundary layer thickness and increase the velocity gradient near the wall.     

Progress in the Use of Non-Linear Two-Equation Models in the Computation of Convective Heat-Transfer in Impinging and Separated Flows

The paper considers the application of the Craft et al. [6]non-linear eddy-viscosity model to separating and impinging flows. The original formulation was found to lead to numerical instabilities when applied to flow separating from a sharp corner. An alternative formulation for the variation of the turbulent viscosity parameterc _μ with strain rate is proposed which, together with a proposed improvement in the implementation of the non-linear model, removes this weakness. It does, however, lead to worse predictions in an impinging jet, and a further modification in the expression for c _μ is proposed, which both retains the stability enhancements and improves the prediction of the stagnating flow. The Yap [24] algebraic length-scale correction term, included in the original model, is replaced with a differential form, developed from that proposed by Iacovides and Raisee [10]. This removes the need to prescribe the wall-distance, and is shown to lead to superior heat-transfer predictions in both an abrupt pipe flow and the axisymmetric impinging jet. One predictive weakness still, however, remains. The proposed model, in common with other near-wall models tested for the abrupt pipe expansion, returns a stronger dependence of Nusselt number on the Reynolds number than that indicated by the experimental data. This revised version was published online in July 2006 with corrections to the Cover Date.     

Heat and mass transfer study of impinging turbulent premixed flames

 Impinging jet combusting flows on granite plates are studied. A mathematical model for calculating heat release in turbulent impinging premixed flames is developed. The combustion including radiative heat transfer and local extinction effects, and flow characteristics are modeled using a finite volume computational approach. Two different eddy viscosity turbulence models, namely the standard k–ɛ and the RNG k–ɛ model with and without radiation (discrete transfer model) are assessed. The heat released predictions are compared with experimental data and the agreement is satisfactory only when both radiative heat transfer and local extinction modeling are taken into account. The results indicate that the main effect of radiation is the decrease of temperature values near the jet stagnation point and along the plate surface. Radiation increases temperature gradients and affects predicted turbulence levels independently of the closure model used. Also, the RNG k–ɛ predicts higher temperatures close the solid plate, with and without radiative heat transfer. Received on 13 November 2000 / Published online: 29 November 2001     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

Stereo-PIV measurements of vapor-induced flow modifications in confined jet impingement boiling

A single subcooled jet of water which undergoes boiling upon impingement on a discrete heat source is studied experimentally using time-resolved stereo particle image velocimetry (PIV). The impinging jet issues from a 3.75 mm diameter sharp-edged orifice in a confining orifice plate positioned 4 orifice diameters from the target surface. The behavior at jet Reynolds numbers of 5,000 and 15,000 is compared for a constant jet inlet subcooling of 10 °C. Fluorescent illumination allows for simultaneous imaging of both the flow tracers and the vapor bubbles in the flow. Flow structure, time-averaged velocities, and turbulence statistics are reported for the liquid regions within the confinement gap for a range of heat inputs at both Reynolds numbers, and the effect of the vapor generation on the flow is discussed. Vapor generation from boiling is found to modify the liquid velocities and turbulence fluctuations in the confinement gap. Flow in the confinement gap is dominated by vapor flow, and the vapor bubbles disrupt both the vertical impinging jet and horizontal wall jet flow. Moreover, vapor bubbles are a significant source of turbulence kinetic energy and dissipation, with the bubbly regions above the heated surface experiencing the most intense turbulence modification. Spectral analysis indicates that a Strouhal number of 0.023 is characteristic of the interaction between bubbles and turbulent liquid jets.     

湍流冲击射流流动与传热的数值研究进展

湍流冲击射流由于其冲击表面时具有很高的局部传热率和冲击力,被广泛应用于如表面的加热、电子元件的冷却、纸张的干燥和材料的切割等工程应用和工业过程中.由于其流动的复杂性,也常被作为一种理想的测试实例来评价湍流模型的性能.此外,湍升力射流与地面之间的空气动力作用对V/STOL (垂直或短距离起落)飞机的性能具有很大的影响.长期以来,人们从理论分析、实验测量和数值模拟方面对冲击射流进行了广泛而系统的研究,积累了丰富的资料.本文在分析了湍流冲击射流的数值研究现状的基础上,对近年来有关湍流冲击射流流动与传热的数值研究方面的文献有选择地进行了综述,重点评述了不同湍流模型对冲击射流流动与传热的预测能力,讨论了存在的问题并对该领域今后的研究方向进行了展望.      

Large Eddy Simulation of a plane impinging jet

Large Eddy Simulations of a plane turbulent impinging jet have been carried out using the dynamic Smagorinsky model. The statistical results are first validated with the measurements from the literature: mean and turbulent quantities along the jet axis and at different vertical locations are presented. This study is completed by the analysis of the wall shear stress at the impingement wall. The effect of the jet Reynolds number (3000Re13500) on the kinematic development of the jet is also discussed. To cite this article: F. Beaubert, S. Viazzo, C. R. Mecanique 330 (2002) 803–810.     

Numerical simulation of turbulent free surface flow with two‐equation k–ε eddy‐viscosity models

This paper presents a finite difference technique for solving incompressible turbulent free surface fluid flow problems. The closure of the time‐averaged Navier–Stokes equations is achieved by using the two‐equation eddy‐viscosity model: the high‐Reynolds k–ε (standard) model, with a time scale proposed by Durbin; and a low‐Reynolds number form of the standard k–ε model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non‐linear terms, a second/third‐order upwinding technique is adopted. The computational method is validated by applying it to the flat plate boundary layer problem and to impinging jet flows. The method is then applied to a turbulent planar jet flow beneath and parallel to a free surface. Computations show that the high‐Reynolds k–ε model yields favourable predictions both of the zero‐pressure‐gradient turbulent boundary layer on a flat plate and jet impingement flows. However, the results using the low‐Reynolds number form of the k–ε model are somewhat unsatisfactory. Copyright © 2004 John Wiley & Sons, Ltd.     

Effects of wall permeability on turbulence

In order to understand the effects of the wall permeability on turbulence near a porous wall, flow field measurements are carried out for turbulent flows in a channel with a porous bottom wall by a two-component particle image velocimetry (PIV) system. The porous media used are three kinds of foamed ceramics which have almost the same porosity (0.8) but different permeability. It is confirmed that the flow becomes more turbulent over the porous wall and tends to be turbulent even at the bulk Reynolds number of Re_b=1300 in the most permeable wall case tested. Corresponding to laminar to turbulent transition, the magnitude of the slip velocity on the porous wall is found to increase drastically in a narrow range of the Reynolds number. To discuss the effects of the wall roughness and the wall permeability, detailed discussions are made of zero-plane displacement and equivalent wall roughness for porous media. The results clearly indicate that the turbulence is induced by not only the wall roughness but the wall permeability. The measurements have also revealed that as Re_b or the wall permeability increases, the wall normal fluctuating velocity near the porous wall is enhanced due to the effects of the wall permeability. This leads to the increase of the turbulent shear stress resulting in higher friction factors of turbulence over porous walls.     

Numerical simulation of turbulent impinging jet on a rotating disk

The calculations of quasi‐three‐dimensional momentum equations were carried out to study the influence of wall rotation on the characteristics of an impinging jet. The pressure coefficient, the mean velocity distributions and the components of Reynolds stress are calculated. The flow is assumed to be steady, incompressible and turbulent. The finite volume scheme is used to solve the continuity equation, momentum equations and k–ε model equations. The flow characteristics were studied by varying rotation speed ω for 0?ω?167.6 rad/s, the distance from nozzle to disk (H/d) was (3, 5, 8 and 10) and the Reynolds number Re base on V_J and d was 1.45 × 10⁴. The results showed that, the radial velocity and turbulence intensity increase by increasing the rotation speed and decrease in the impingement zone as nozzle to disk spacing increases. When the centrifugal force increases, the radial normal stresses and shear stresses increase. The location of maximum radial velocity decreases as the local velocity ratio (α) increases. The pressure coefficient depends on the centrifugal force and it decreases as the distance from nozzle to plate increases. In impingement zone and radial wall jet, the spread of flow increases as the angular velocity decreases The numerical results give good agreement with the experiment data of Minagawa and Obi (Int. J. of Heat and Fluid Flow 2004; 25 :759–766). Copyright © 2006 John Wiley & Sons, Ltd.     

Automatic eduction and statistical analysis of coherent structures in the wall region of a confined plane turbulent impinging jet

This paper describes a vortex detection algorithm used to expose and statistically characterize the coherent flow patterns observable in the velocity vector fields measured by particle image velocimetry in the impingement region of air curtains. The philosophy and the architecture of this algorithm are presented. Its strengths and weaknesses are discussed. The results of a parametrical analysis performed to assess the variability of the response of our algorithm to the three user-specified parameters in our eduction scheme are reviewed. The technique is illustrated in the case of a plane turbulent impinging twin-jet with an opening ratio of 10. The corresponding jet Reynolds number, based on the initial mean flow velocity U ₀ and the jet width e, is 14,000. The results of a statistical analysis of the size, shape, spatial distribution and energetic content of the coherent eddy structures detected in the impingement region of this test flow are provided. Although many questions remain open, new insights into the way these structures might form, organise and evolve are given. Relevant results provide an original picture of the plane turbulent impinging jet.     

Evaluation of erosion in a two‐way coupled fluid–particle system

In this study, the effects of flow turbulence intensity, temperature, particle sizes and impinging velocity on erosion by particle impact are demonstrated numerically. Underlying turbulent flow on an Eulerian frame is described by the compressible Reynolds averaged Navier–Stokes equations with a RNG k–ε turbulence model. The particle trajectories and particle–wall interactions are evaluated by a Eulerian–Lagrangian approach in a two‐way coupling system. An erosion model considering material weight removal from surfaces is used to predict erosive wear. Computational validation against measured data is demonstrated satisfactorily. The analysis of erosion shows that the prevention of erosion is enhanced by increasing the effects of flow temperature and turbulence intensity and reducing particle inertial momentum. Copyright © 2001 John Wiley & Sons, Ltd.