The application of thin-film technology to measure turbine-vane heat transfer and effectiveness in a film-cooled, engine-simulated environment

Thin-film technology has been used to measure the heat transfer coefficient and cooling effectiveness over heavily film cooled nozzle guide vanes (NGVs). The measurements were performed in a transonic annular cascade which has a wide operating range and simulates the flow in the gas turbine jet engine. Engine-representative Mach and Reynolds numbers were employed and the upstream free-stream turbulence intensity was 13%. The aerodynamic and thermodynamic characteristics of the coolant flow (momentum flux and density ratio between the coolant and mainstream) have been modelled to represent engine conditions by using a foreign gas mixture of SF₆ and Argon. Engine-level values of heat transfer coefficient and cooling effectiveness have been obtained by correcting for the different molecular (thermal) properties of the gases used in the engine-simulated experiments to those which exist in the true engine environment. This paper presents the best combined heat transfer coefficient and effectiveness data currently available for a fully cooled, three-dimensional NGVs at engine conditions.     

Experimental and numerical investigation on heat transfer characteristics of supercritical R1234ze(E) cooled in horizontal tubes

The heat transfer characteristics of supercritical fluids in tubes have been considered indispensable for the design and optimization of the heat exchanger and the energy conversion system. Specifically the cooling heat transfer of supercritical R1234ze(E) in horizontal tubes is a promising heat-power conversion technology; however, there is a scarcity of conducted research in available literature. The present work, the first-ever study in this direction, aims to thoroughly investigate the heat transfer characteristics of supercritical R1234ze(E) which is cooled in horizontal tubes. Experimental work was performed to thoroughly explore and inspect the heat transfer characteristics of supercritical R1234ze(E) passing through the tube of 4.12 mm diameter at 4–5 MPa pressure and 240–400 kg/m² s mass flux. Furthermore, the simulation study, supporting the experimental investigation under the same conditions of pressure and mass flux, extended the range of tube diameter up to 9.44 mm. The effects of pressure, mass flux and tube diameter on the heat transfer coefficient were carefully analyzed in the present research work. Based on the simulation results and experimental results, heat transfer correlations were newly developed by separating the region above and below the pseudo-critical temperature. The average absolute deviation between the calculated Nusselt numbers by the numerical correlation and the simulation results was found 2.87%; the average absolute deviation between the calculated Nusselt numbers by the experimental correlation and the experimental results was found 5.3%.     

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.     

The effect of the low Reynolds number k-e turbulence models on simulation of the enhanced and deteriorated convective heat transfer to the supercritical fluid flows

One of the most important challenges in mathematical modeling of convective heat transfer to the turbulent flow of a supercritical fluid is the turbulence modeling. The turbulence modeling, like other aspects of the supercritical fluid flows, seems to be high affected by the large variations of the fluid properties which needs to be further investigated. A two dimensional CFD code has been developed and used in this study to examine a number of the low Reynolds number k-e turbulence models. Both flow conditions corresponding to the heat transfer enhancement and deterioration have been studied. The results appear to be quite sensitive to the choice of the turbulence model, especially in the deteriorated regime of heat transfer. The turbulence model assisting the two-dimensional numerical model of the present study to best fit the experiments has been determined for both cases of the enhanced and deteriorated heat transfer. That is while the jump in the wall temperature occurring in the deteriorated regime of heat transfer is over- predicted by the present numerical code regardless of the turbulence model used.     

Numerical Study of Mixed Convection in Vertical Tubes under Conditions of Unstable Stratification

The present numerical study was performed using a three-parameter turbulence model supplemented by the thermogravitational terms and the transport equation for the transverse turbulent heat flux. The results obtained are compared with the few available experimental data for downward water and air flows through vertical heated tubes and with calculations for upward flow.     

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.     

A DNS study on the effects of convex streamwise curvature on coherent structures in a temporally-developing turbulent boundary layer with supercritical water

Direct numerical simulation (DNS) results are used to establish the effect of convex streamwise curvature on the development of turbulent boundary layers, and the effect of such curvature on the forced-convection heat transfer variations observed at certain supercritical thermodynamic states. The results illustrate the stabilizing effects of this flow geometry through modification of the structure and distribution of hairpin-like vortical flow structures in the boundary layer. Furthermore, enhancement of convective heat transfer realized at a particular heat flux-to-mass flux ratio with the working fluid at a supercritical state is observed to be reduced by the stabilizing effect of convex surface curvature.     

Analysis on heat transfer correlations of supercritical CO<Subscript>2</Subscript> cooled in horizontal circular tubes

The objective of the present study is to analyze the heat transfer correlations of supercritical CO₂ cooled in horizontal circular tubes. In the paper, heat transfer correlations are first reviewed and compared with the experimental data at different heat fluxes. The results show that most of the previous correlations agree well with the experimental data under lower heat flux, but fail to predict the heat transfer coefficient well when the heat flux is as high as 33 kW/m². The study of buoyancy effect on convective heat transfer shows that buoyancy effect significantly affects the heat transfer with the increase of heat flux, and both free and forced convections operate in the turbulence flow during supercritical CO₂ cooling process. The influencing factors on heat transfer coefficient are summarized and the new correlation can be developed with the four dimensionless numbers.     

Experimental study of turbulent round jet flow impinging on a square cylinder laid on a flat plate

A turbulent axisymmetric air jet impinging on a square cylinder mounted on a flat plate has been studied experimentally. Turbulence statistics and flow’s topology were investigated. When the surface was heated through uniform heat flux, local heat transfer coefficient was measured. The jet from a long round pipe, 75 pipe diameters (D) in length, at Reynolds number of 23,000, impinged vertically on the square cylinder (3D × 3D × 43D). Measurements were performed using particle image velocimetry, flow visualization using fluorescent dye and infrared thermography. The flow’s topology demonstrated a three-dimensional recirculation after separating from the square cylinder and a presence of foci between the bottom corner and the recirculation’s detachment line. The distribution of heat transfer coefficient was explained by the influence of these flow’s structures and the advection of kinetic energy. On the impingement wall of the square cylinder, a secondary peak in heat transfer coefficient was observed. Its origin can be attributed to very pronounced shear production coupled with the external turbulence coming from the free jet.     

Comparisons and applications of numerical simulation methods for predicting aerodynamic heating around complex configurations

Numerical simulation methods of aerodynamic heating were compared by considering the inuence of numerical schemes and turbulence models,and attempting to investigate the applicability of numerical simulation methods on predicting heat flux in engineering applications. For some typical cases provided with detailed experimental data,four spatial schemes and four turbulence models were adopted to calculate surface heat flux. By analyzing and comparing,some inuencing regularities of numerical schemes and turbulence models on calculating heat flux had been acquired. It is clear that AUSM+-up scheme with rapid compressibilitymodified high Reynolds number k-ω model should be appropriate for calculating heat flux. The numerical methods selected as preference above were applied to calculate the heat flux of a 3-D complex geometry in high speed turbulent flows. The results indicated that numerical simulation can capture the complex flow phenomena and reveal the mechanism of aerodynamic heating. Especially,the numerical result of the heat flux at the stagnation point of the wedge was well in agreement with the prediction of Kemp-Riddel formula,and the surface heat flux distribution was consistent with experiment results,which implied that numerical simulation can be introduced to predict heat flux in engineering applications.     

Direct numerical simulation of turbulent heat transfer in a wall-normal rotating channel flow

Investigations into the characteristics of turbulent heat transfer and coherent flow structures in a plane-channel subjected to wall-normal system rotation are conducted using direct numerical simulation (DNS). In order to investigate the influence of system rotation on the temperature field, a wide range of rotation numbers are tested, with the flow pattern transitioning from being fully turbulent to being quasilaminar, and eventually, fully laminar. In response to the Coriolis force, secondary flows appear as large vortical structures, which interact intensely with the wall shear layers and have a significant impact on the distribution of turbulence kinetic energy (TKE), turbulence scalar energy (TSE), temperature statistics, and turbulent heat fluxes. The characteristic length scales of turbulence structures responsible for the transport of TSE are the largest at the quasilaminar state, which demands a very large computational domain in order to capture the two-dimensional spectra of temperature fluctuations. The effects of the Coriolis force on the turbulent transport processes of the temperature variance and turbulent heat fluxes are thoroughly examined in terms of their respective budget balances.     

The relationship between coherent structures and heat transfer processes in the initial region of a round jet

 This paper describes an experimental study of the relationship between coherent vortical structures and the intensity of heat transport in the initial region of a round, free jet. Simultaneous measurements of velocity and temperature were taken with a four-wire combined probe in a jet that was acoustically stimulated with a frequency corresponding to the jet-column mode. The obtained results suggest that the mutual phase relations between oscillatory and random components of velocity and temperature lead to substantial intensification of the radial heat transport. Due to the same reason the longitudinal heat flux does not reveal a significant change in the presence of coherent structures and, as a result, a much wider spread of the temperature field in comparison with velocity may be observed as a characteristic feature of this flow. It was also observed that heat transfer processes are realized in substantial part by random turbulence generated due to the action of coherent motion. Received: 16 January 1997/Accepted: 20 August 1997     

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

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

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Conjugate heat transfer study of incompressible turbulent offset jet flows

In the present case, the conjugate heat transfer involving a turbulent plane offset jet is considered. The bottom wall of the solid block is maintained at an isothermal temperature higher than the jet inlet temperature. The parameters considered are the offset ratio (OR), the conductivity ratio (K), the solid slab thickness (S) and the Prandtl number (Pr). The Reynolds number considered is 15,000 because the flow becomes fully turbulent and then it becomes independent of the Reynolds number. The ranges of parameters considered are: OR = 3, 7 and 11, K = 1–1,000, S = 1–10 and Pr = 0.01–100. High Reynolds number two-equation model (k–ε) has been used for turbulence modeling. Results for the solid–fluid interface temperature, local Nusselt number, local heat flux, average Nusselt number and average heat transfer have been presented and discussed.     

Turbulence characteristics of radially-confined impinging jet flows

Radially confined, axisymmetric impinging jet flows are investigated by using the standard particle image velocimetry experimental technique. The confinement is achieved by placing a confinement block around a jet, co-axially. The inner diameter of the block is successively varied to nine different values. The inlet-based Reynolds number of the jet is kept constant at 5000. The nine diametric values yielded nine different flows of widely different characteristics. Among other usage, an insight into the flow characteristics can be helpful in designing compact impinging jet applications, as such a radially confined flow is equivalent to passing the pre-impingement jet through a hole perforated in a solid wall (i.e. the jet source can be placed behind a wall). The study has revealed that the flows, in general, form two circulation zones, three mixing layers, and two boundary layers. Based on turbulence characteristics of the five shear layers, overall characteristics of the flows are understood systematically. Mean velocity and various turbulence statistics are also presented, and mechanisms underlying behind their variations are explained. Finally, scaling laws are obtained for the mean velocity and for the turbulence statistics, both in the impingement and in the wall jet regions.     

Buoyancy effects on turbulent swirling flows

A three-parameter model of turbulence applicable to free boundary layers has been developed and applied for the prediction of axisymmetric turbulent swirling flows in uniform and stagnant surroundings under the action of buoyancy forces. The turbulent momentum and heat fluxes appearing in the time-averaged equations for the mean motion have been determined from algebraic expressions, derived by neglecting the convection and diffusion terms in the differential transport equations for these quantities, which relate the turbulent fluxes to the kinetic energy of turbulence, k, the dissipation length scale of turbulence, L, and the temperature covariance, T ′². Differential transport equations have been used to determine these latter quantities. The governing equations have been solved using fully implicit finite difference schemes. The turbulence model is capable of reproducing the gross features of pure jet flows, buoyant flows and swirling flows for weak and moderate swirl. The behaviour of a turbulent buoyant swirling jet has been found to depend solely on exit swirl and Froude numbers. The predicted results indicate that the incorporation of buoyancy can cause significant changes in the behaviour of a swirling jet, particularly when the buoyancy strength is high. The jet exhibits similarity behaviour in the initial region for weak swirl and weak buoyancy strengths only, and the asymptotic case of a swirling jet under the action of buoyancy forces is a pure plume in the far field. The predicted results have been found to be in satisfactory agreement with the available experimental data and in good qualitative agreement with other predicted results.     

Computational Modelling of Turbulent Mixing in Confined Swirling Environment Under Constant and Variable Density Conditions

A high-intensity swirling flow in a model combustor subjected to large density variations has been examined computationally. The focus is on the Favre-averaged Navier–Stokes computations of the momentum and scalar transport employing turbulence models based on the differential second-moment closure (SMC) strategy. An updated version of the basic, high-Reynolds number SMC model accounting for a quadratic expansion of both the pressure–strain and dissipation tensors and a near-wall SMC model were used for predicting the mean velocity and turbulence fields. The accompanied mixing between the annular swirling air flow and the central non-swirling helium jet was studied by applying three scalar flux models differing mainly in the model formulation for the pressure-scalar gradient correlation. The computed axial and circumferential velocities agree fairly well with the reference experiment [So et al., NASA Contractor Report 3832, 1984; Ahmed and So, Exp. Fluids 4 (1986) 107], reproducing important features of such a weakly supercritical flow configuration (tendency of the flow core to separate). Although the length at which the mixing was completed was reproduced in reasonable agreement with the experimental results, the mixing activity in terms of the spreading rate of the shear/mixing layer, that is its thickness, was somewhat more intensive. Prior to these investigations, the model applied was validated by computing the transport of the passive scalar in the non-swirling (Johnson and Bennet, Report NASA CR-165574, UTRC Report R81-915540-9, 1981) and swirling (Roback and Johnson, NASA Contractor Report 168252, 1983) flow in a model combustor.     

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.     

Experimental study of heat flux in mixed convective flow over solid waves

In this experimental study, we address transport processes in a mixed convective flow over a heated wavy surface. Therefore, we combine digital particle image velocimetry (DPIV) and two-color planar laser induced fluorescence (PLIF) to simultaneously measure the velocity and temperature field. For this, we propose to use the dye combination Rhodamine B and Rhodamine 110, both excited with the Nd:YAG laser also used for the PIV measurements. We investigate the influence of mixed convection over a wavy surface on the velocity field, turbulence statistics, the temperature field and the heat flux. By computing these quantities we find a correlation between the maximum in the Reynolds stress profiles and the components of the heat flux vector, thus regions of maximum momentum and scalar transport coincide. In addition, we apply a proper orthogonal decomposition (POD) to extract the most dominant flow structures in a measurement plane above the wavy surface. This first POD mode is identified as streamwise-oriented, counter-rotating vortices whose spanwise scaling is also correlated with the maximum of heat flux.