Lattice Boltzmann simulation of the convective heat transfer from a stream-wise oscillating circular cylinder

In this paper, we studied the convective heat transfer from a stream-wise oscillating circular cylinder. Two dimensional numerical simulations are conducted at Re = 100–200, A = 0.1–0.4 and F = f_o/f_s = 0.2–3.0 with the aid of the lattice Boltzmann method. In particular, detailed attentions are paid on the extensive numerical results elucidating the influence of oscillation frequency, oscillation amplitude and Reynolds number on the time-average and RMS value of the Nusselt number. Over the ranges of conditions considered herein, the heat transfer characteristics are observed to be influenced in an intricate manner by the value of the oscillation frequency (F), oscillation amplitude (A) and Reynolds number (Re). Firstly, the heat transfer is enhanced when the cylinder oscillates stream-wise with small amplitude and low frequency, while it will be reduced by large amplitude and high frequency. Secondly, the average Nusselt number (Nu (ave)) decreases against the increasing value of oscillation frequency, while the RMS value of the Nusselt number, Nu (RMS), displays an opposite trend. Third, we obtained a similar frequency effect on the heat transfer over the range of Reynolds numbers investigated in this paper. In addition, detailed analyses on phase portraits, energy spectrum are also made.     

Film cooling effectiveness from trenched shaped and compound holes

This paper presents a comparative-numerical investigation on film cooling from a row of simple and compound-angle holes injected at 35° on a flat plate with four film cooling configurations: (1) cylindrical film hole; (2) 15° forward diffused film hole; (3) trenched cylindrical film hole; (4) trenched 15° forward-diffused film hole. All simulations are at fixed density ratio of 1.6, blowing ratio of 1.25, length-to-diameter L/D = 4 and pitch-to-diameter ratio of 3.0. The effect of length-to-diameter ratio on film cooling has been also investigated using L/D in the range of 1–8. Computational solutions of the steady, Reynolds-averaged Navier–Stokes equations have been obtained using a finite volume method. It has been found that the shape of the hole and the trenched holes can significantly affect the film cooling flow over the protected surface. Further, it has been shown that the film cooling effectiveness by trenched shaped holes is higher than all other configurations both in spanwise and streamwise specially downstream of the injection. Also, a trenched compound angle injection shaped hole produces much higher film cooling protection than the other configurations investigated in the present paper. The length-to-diameter ratio of trenched holes was found to have a significant effect on film cooling effectiveness and the spread of the coolant jets.     

Free convection heat transfer from a horizontal fin attached cylinder between confined nearly adiabatic walls

Steady state two-dimensional free convection heat transfer from a horizontal, isothermal fin attached cylinder, located between nearly two adiabatic walls is studied experimentally using a Mach–Zehnder interferometer. Effects of the walls inclination angel (θ) on heat transfer from the cylinder is investigated for Rayleigh number ranging from 1000 to 15,500. Two cylinders with different diameters of D = 10 and 20 mm are used to cover wide Rayleigh range. Results indicate that, heat transfer phenomena differ for different Rayleigh number. For Rayleigh numbers lower than 5500, heat transfer rate from cylinder surface is lower than the heat transfer from a single cylinder. In this range by the use of walls, heat transfer from the cylinder decreases slightly and walls’ inclination does not change heat transfer rate from the cylinder surface. For Rayleigh number ranging from 5500 to 15,500, amount of heat transfer from the cylinder surface is less than that of a single cylinder. However, by adding nearly adiabatic walls to experimental model heat transfer mechanism differs and chimney effect between fin and walls increases the heat transfer rate from the cylinder surface. By increasing the walls inclination angel from 0° to 20°, the chimney effect between walls and fin diminishes and heat transfer rate from the cylinder surface is approaching to the heat transfer rate of fin attached cylinder without adiabatic walls.     

Film cooling on a convex surface: influence of external pressure gradient and Mach number on film cooling performance

 The film cooling performance on a convex surface subjected to zero and favourable pressure gradient free-stream flow was investigated. Adiabatic film cooling effectiveness values were obtained for five different injection geometries, three with cylindrical holes and two with shaped holes. Heat transfer coefficients were derived for selected injection configurations. CO₂ was used as coolant to simulate density ratios between coolant and free-stream close to gas turbine engine conditions. The film cooling effectiveness results indicate a strong dependency on the free-stream Mach number level. Results obtained at the higher free-stream Mach number show for cylindrical holes generally and for shaped holes at moderate blowing rates significant higher film cooling effectiveness values compared to the lower free-stream Mach number data. Free-stream acceleration generally reduced adiabatic film cooling effectiveness relative to constant free-stream flow conditions. The different free-stream conditions investigated indicate no significant effects on the corresponding heat transfer increase due to film injection. The determined heat flux ratios or film cooling performance indicated that coolant injection with shaped film cooling holes is much more efficient than with cylindrical holes especially at higher blowing rates. Heat flux penalties can occur at high blowing rates when using cylindrical holes. Received on 29 May 2000     

Effect of film cooling on the aerodynamic performance of an airfoil

The effect of film cooling on the aerodynamic performance of turbine blades is becoming increasingly important as the gas turbine operating temperature is being increased in order to increase the performance. The current paper investigates the effect of blowing ratio on the aerodynamic losses of a symmetric airfoil by pressure measurements and Particle Image Velocimetry (PIV). The test model features 4 rows of holes located on the suction side at 5%, 10%, 15% and 50% of the chord length. The Reynolds number based on the airfoil chord is 1.2 × 10⁵. Experiments are performed by varying the location of air injection, the angle of attack, and the mainstream velocity. The coolant air is injected at ambient temperature and the blowing ratio is varied from 0 to 1.91. It is observed that the losses due to film cooling increase with blowing ratio of 0 to 0.48, and the wake is shifted towards the suction side. Conversely, the aerodynamic losses decrease when the blowing ratio is increased further from 0.64 to 1.91. This trend has been observed for all the experimental configurations. The effect of blowing ratio on flow separation is investigated with the time-averaged velocity fields obtained from PIV measurements. It is observed that low blowing ratios, the separation point shifts upstream and at high blowing ratios the ejected coolant energizes the flow and delays separation. The pressure field around the airfoil is reconstructed from the integration of the Poisson equation based on the PIV velocity fields. The experimental results can be used for validation of numerical models for predicting losses due to film cooling.     

Heat-transfer characteristics of climbing film evaporation in a vertical tube

Heat-transfer characteristics of climbing film evaporation were experimentally investigated on a vertical climbing film evaporator heated by tube-outside hot water. The experimental setup was designed for determining the effect of the height of feed water inside a vertical tube and the range of temperature difference on local heat transfer coefficient inside a vertical tube (h_i). In this setup, the height of feed water was successfully controlled and the polypropylene shell effectively impedes the heat loss to the ground. The results indicated that a reduction in the height of feed water contributed to a significant increase in h_i if no dry patches around the wall of the heated tube appeared inside the tube. The height ratio of feed water R_h = 0.3 was proposed as the optimal one as dry patches destroyed the continuous climbing film when R_h is under 0.3. It was found that the minimum temperature difference driving climbing film evaporation is suggested as 5 °C due to a sharp reduction in h_i for temperature difference below 5 °C. The experiment also showed that h_i increased with an increase in temperature difference, which proved the superiority of climbing film evaporation in utilizing low-grade surplus heating source due to its wide range of driving temperature difference. The experimental results were compared with the previous literature and demonstrated a satisfactory agreement.     

Heat transfer due to impinging co-axial jets and the jets’ fluid flow characteristics

This investigation had multiple goals. One goal was to obtain definitive information about the heat transfer characteristics of co-axial impinging jets, and this was achieved by measurements of the stagnation-point, surface-distribution and average heat transfer coefficients. These results are parameterized by the Reynolds number Re which ranged from 5000 to 25,000, the dimensionless separation distance between the jet exit and the impingement plate H/D (4–12), and the ratio of the inner diameters of the inner and outer pipes d/D (0–0.55). The d/D = 0 case corresponds to a single circular jet. The other major goal of this work was to quantify the velocity field of co-axial free jets (impingement plate removed). The velocity-field study included both measurements of the mean velocity and the turbulence intensity.It was found that the variation of the stagnation-point heat transfer coefficient with d/D attained a maximum at d/D = 0.55. Furthermore, the variation of the local heat transfer coefficient across the impingement surface was more peaked for d/D = 0 and became flatter with decreasing d/D. This suggests that for cooling a broad expanse of surface, co-axial jets of high d/D are preferable. On the other hand, for localized cooling, the single jet (d/D = 0) performed the best. In general, for a given Reynolds number, a co-axial jet yields higher heat transfer coefficients than a single jet. Off-axis velocity peaks were encountered for the jets with d/D = 0.105. The measurements of turbulence intensity yielded values as high as 18%.     

Flow boiling heat transfer of R134a,R236fa and R245fa in a horizontal 1.030 mm circular channel

This research focuses on acquiring accurate flow boiling heat transfer data and flow pattern visualization for three refrigerants, R134a, R236fa and R245fa in a 1.030 mm channel. We investigate trends in the data, and their possible mechanisms, for mass fluxes from 200 to 1600 kg/m²s, heat fluxes from 2.3 kW/m² to 250 kW/m² at T_sat = 31 °C and ΔT_sub from 2 to 9 K. The local saturated flow boiling heat transfer coefficients display a heat flux and a mass flux dependency but no residual subcooling influence. The changes in heat transfer trends correspond well with flow regime transitions. These were segregated into the isolated bubble (IB) regime, the coalescing bubble (CB) regime, and the annular (A) regime for the three fluids. The importance of nucleate boiling and forced convection in these small channels is still relatively unclear and requires further research.     

Conjugate turbulent heat transfer in a square cavity with a solar control coating deposited to a vertical semitransparent wall

This paper presents a numerical study of the conjugate heat transfer (natural convection, surface thermal radiation and conduction) in a square cavity with turbulent flow. The cavity has one vertical isothermal wall, two horizontal adiabatic walls and one vertical semitransparent wall with a selective coating applied to the inner side to control the solar radiation transmission. Later on the semitransparent wall is replaced with another one without the selective coating. The mathematical model for the turbulent flow in the cavity was solved using the finite volume method. The system had the following conditions: the uniform temperature in the isothermal wall was 21 °C, the external ambient temperature was fixed at 35 °C and on the semitransparent wall the direct normal solar irradiation of 750 W/m² was considered constant. The Rayleigh number was varied in the range of 10⁹ ? Ra ? 10¹² by changing the lengths of the cavity from 0.70 m to 6.98 m, respectively. The results show that, even though the air temperature of the cavity with the solar control film coating semitransparent wall (case A) is higher compared with the one without solar film coating (case B), the total amount of heat going through the cavity is lower compared to the one going through the cavity without solar control film. The total amount of energy transferred to the air in cavity for the case A was 41.98% less than for the case B. A set of correlations for the Nusselt number was obtained for both cases considering the conjugate heat transfer.     

Influence of single rectangular groove on the flow past a circular cylinder

In the present study, flow control mechanism of single groove on a circular cylinder surface is presented experimentally using Particle image velocimetry (PIV). A square shaped groove is patterned longitudinally on the surface of the cylinder with a diameter of 50 mm. The flow characteristics are studied as a function of angular position of the groove from the forward stagnation point of the cylinder within 0° ≤ θ ≤ 150°. In the current work, instantaneous and time-averaged flow data such as vorticity, ω streamline, Ψ streamwise, u/U_o and transverse, v/U_o velocity components, turbulent kinetic energy, TKE and RMS of streamwise, u_rms and transverse, v_rms velocity components are utilized in order to present the results of quantitative analyses. Furthermore, Strouhal numbers are calculated using Karman vortex shedding frequency, f_k obtained from single point spectral analysis. It is concluded that a critical angular position of the groove, θ = 80° is observed. The flow separation is controlled within 0° ≤ θ < 80°. At θ = 80°, the flow separation starts to occur in the upstream direction. The instability within the shear layer is also induced on grooved side of the cylinder with frequencies different than Karman vortex shedding frequency, f_k.     

Experimental study on turbulent mixing of spray droplets in crossflow

Centrifugal spray injected at various angles in gas crossflow has been studied experimentally using PIV visualization system and image-processing techniques. Experiments were carried out inside a rectangular duct (95 mm × 95 mm in cross-section) at ambient temperature and pressure, with different gas Reynolds numbers (vary from 12,900 to 45,000) and three injection angles (60°, 90° and 120°). The spray angle of the centrifugal nozzle is 80°, with D₃₂ of 80 μm. The instantaneous images of droplets distribution and the values of the line-averaged D₃₂ at different positions on the cross-sections along the flow field for each condition were obtained, and their flow field configurations were achieved. Quantitative assessments of mixing degree between two phases for different injection angles were determined using a spatial unmixedness parameter. It is found that the addition of droplets into the gas crossflow enhanced the turbulence intensity of the gas crossflow and caused different-scale vortices. The flow field structure, to a great extent, is dependent on the injection angle. The entrainment and centrifugal force of large vortex lead to uneven droplet distribution and moreover influence the mixing of droplets and gas crossflow. A better mixing result can be obtained with the injection angle of nozzles of 60°.     

Analysis of upstream,double-row,cylindrical holes on primary and secondary effects of endwall flow and film cooling

Flow features and film cooling performance of five configurations of double-row, cylindrical holes, upstream of an E3 vane, in a linear cascade are numerically investigated. This simulation is completed using a verified turbulence model at four blowing ratios (M = 0.5, 1.0, 1.5, 2.0). The first three configurations have two rows of cylindrical holes, each row with the same compound angle (β=-45°, 0° or 45°), while the other two have two rows with opposite compound angles (β=-45°, 45° and β=45°, -45°), which are also referred to as double-jet film cooling (DJFC) holes. The primary effects on the downstream endwall and the secondary effects on the nearby airfoil of the cooled passage are analyzed and discussed in detail. Results show that at low blowing ratios the movement of the coolant is denominated by the interaction between the jets and vortices resulting in similar film coverage on both the endwall and airfoil. The effect of vortices is reduced at high blowing ratios. It is also shown that the movement of the coolant is determined by the initial velocity direction, as well as the film cooling configuration.     

Augmentation of heat transfer from heat source placed downstream a guide fence: An experimental study

The present study investigated experimentally the heat transfer from a heat source simulating an electronic chip mounted on a printed circuit board placed downstream of a guide fence on the lower wall of the flow passage with two different aspect ratios (H/W = 0.3 and 1). The channel height to the heat source height ratios (H/B) are of 10 and 3. The effect of the guide fence height (b) and the spacing between the guide fence and the heat source (S) were investigated. The guide fence was orientated such that guide fence extension point was varied from the midpoint of the front face of the heat source to the endpoint of the side face at 5000 ? Re_L ? 30,000. The results for the heat source without guide fence displayed noticeable difference when compared with the flow over smooth plate placed on the lower wall of the flow passage. An enhancement in the convective heat transfer coefficient up to 20% is obtained when decreasing the flow passage height to the heat source height ratio from 10 to 3. Also, higher Nusselt number is located at the front face and the vertical sides of the heat source compared with that of the top surface. Nusselt number increases with the increase in both Reynolds number and the guide fence height while the effect of spacing between the guide fence and the heat source depending on the guide fence height. Correlations for the average Nusselt number were obtained utilizing the present measurements within the investigated range of the different parameters.     

Effect of magnetic field on 3D flow and heat transfer during solidification from a melt

We present the effect of a magnetic field on three-dimensional fluid flow and heat transfer during solidification from a melt in a cubic enclosure. The walls of the enclosure are considered perfectly electrically conducting and the magnetic field is applied separately in three directions. The finite-volume method with enthalpy formulation is used to solve the mathematical model in the solid and liquid phases. The results obtained by our computer code are compared with the numerical and experimental data found in the literature. For Gr = 5 × 10⁵ and Ha = 0, 25, 50, 75, and 100 (where Gr and Ha are the Grashof and Hartmann numbers, respectively), the effects of magnetic field on flow and thermal fields, and on solid/liquid interface shape are presented and discussed. The interface is localized with and without magnetic field. The results show a strong dependence between the interface shape and the intensity and orientation of magnetic field. When the magnetic field is applied along the X-direction, the magnetic stability diagrams (V_max–Ha) and (Nu_avg–Ha) show the strongest stabilization of the flow field and heat transfer.     

Effect of side ratio on fluid flow and heat transfer from rectangular cylinders using the PANS method

A computational study of heat transfer from rectangular cylinders is carried out. Rectangular cylinders are distinguished based on the ratio of the length of streamwise face to the height of the cross-stream face (side ratio, R). The simulations were performed to understand the heat transfer in a flow field comprising separation, reattachment, vortex shedding and stagnation. The Partially-Averaged Navier–Stokes (PANS) modeling approach is used to solve the turbulent flow physics associated and the wall resolve approach is used for the near wall treatment because of the flow separation involved. The simulations were performed using a finite volume based opensource software, OpenFOAM, at Reynolds number (Re) = 22,000 for rectangular cylinder at constant temperature kept in an air stream. Two critical side ratios were obtained, R = 0.62 and 3.0. At R = 0.62, the maximum value of the drag coefficient (C_d) = 2.681 was observed which gradually reduced by 54% at R = 4.0. The base pressure coefficient and global Nusselt number also attained the maximum value at R = 0.62 and from R = 2.5 to 3.0 a sharp discontinuous increase by 140% in the Strouhal number was observed. At R = 0.62, it was observed that the separated flow reattaches at the trailing edge after rolling over the side face and therefore increases the overall Nusselt number. The phase averaging was also performed to analyze the unsteady behavior of heat transfer.     

Chebyshev finite difference approach to modeling the thermoviscosity effect in a power-law liquid film on an unsteady stretching surface

The effect of thermoviscosity (temperature-dependent viscosity) on the heat transfer in a power-law liquid film over an unsteady stretching sheet is investigated. Similarity analysis is used to transform the governing equations for mass, momentum and energy into a system of ordinary differential equations, which contain a thermoviscosity parameter θ_r, unsteadiness parameter S, generalized Prandtl number Pr and power-law index n. The film thickness, the temperature distributions, the local heat transfer rate, and the local skin-friction coefficient were obtained using the Chebyshev finite difference method (ChFD). The results show that thermoviscosity significantly increases the film thickness and the local heat transfer rate while decreasing the local skin-friction coefficient as θ_r → 1. It is found that liquids with a higher power-law index will have a larger film thickness and a higher free-surface temperature, which indicate a lower local heat transfer rate, ?θ′(0).     

Preparation of 3D micro/nanostructured CeO2: Influence of organic/inorganic acids

CeO₂ is an important porous material with a wide range of applications in the abatement of volatile organic compounds (VOCs). In this paper, we prepared a series of novel three-dimensional (3D) micro/nanostructured CeO₂ materials via a solvothermal method. Organic acid-assisted synthesis and inorganic acid post-treatment were used to adjust the CeO₂ microstructures. The size of the 3D micro/nanostructures could be controlled in the range from 180 nm to 1.5 μm and the surface morphology changed from rough to smooth with the use of different organic acids. The CeO₂ synthesized with acetic acid featured a hierarchical porosity and showed good performance for toluene catalytic combustion: a T₅₀ of 187 °C and a T₉₀ of 195 °C. Moreover, the crystallite size, textural properties, and surface chemical states could be tuned by inorganic acid modification. After treatment with HNO₃, the modified CeO₂ materials exhibited improved catalytic activity, with a T₅₀ of ∼175 °C and a T₉₀ of ∼187 °C. We concluded that the toluene combustion activity is related to the porosity and the amount of surface active oxygen of the CeO₂. Both these features can be tuned by the co-work of organic and inorganic acids.     

Influence of temperature on a carbon–fibre epoxy composite subjected to static and fatigue loading under mode-I delamination

In this paper, interlaminar crack initiation and propagation under mode-I with static and fatigue loading of a composite material are experimentally assessed for different test temperatures. The material under study is made of a 3501-6 epoxy matrix reinforced with AS4 unidirectional carbon fibres, with a symmetric laminate configuration [0°]_16/S. In the experimental programme, DCB specimens were tested under static and fatigue loading. Based on the results obtained from static tests, fatigue tests were programmed to analyse the mode-I fatigue behaviour, so the necessary number of cycles was calculated for initiation and propagation of the crack at the different temperatures. G–N curves were determined under fatigue loading, N being the number of cycles at which delamination begins for a given energy release rate. G_ICmax–a, a–N and da/dN–a curves were also determined for different G_cr rates (90%, 85%, 75%, etc.) and different test temperatures: 90 °C, 50 °C, 20 °C, 0 °C, ?30 °C and ?60 °C.     

Wavelet decomposition method decoupled boiling/evaporation oscillation mechanisms over two to three timescales: A study for a microchannel with pin fin structure

Boiling/evaporation heat transfer in a microchannel with pin fin structure was performed with water as the working fluid. Simultaneous measurements of various parameters were performed. The chip wall temperatures were measured by a high spatial-time resolution IR image system, having a sensitivity of 0.02 °C. The flow pattern variations synchronously changed wall temperatures due to ultra-small Bi number. The wavelet decomposition method successfully identified the noise signal and decoupled various temperature oscillations with different amplitudes and frequencies. Three types of temperature oscillations were identified according to heat flux q and mass flux G. The first type of oscillation occurred at q/G < 0.62 kJ/kg. The approximation coefficient of wavelet decomposition decided the dominant cycle period which was ∼3 times of the fluid residence time in the microchannel, behaving the density wave oscillation characteristic. The detail coefficients of wavelet decomposition decided the dominant cycle period, which matched the flow pattern transition determined value well, representing the flow pattern transition induced oscillation. For the second type of oscillation, the wavelet decomposition decoupled the three oscillation mechanisms. The pressure drop oscillation caused the temperature oscillation amplitudes of 5–10 °C and cycle periods of 10–15 s. The density wave oscillation and flow pattern transition induced oscillation are embedded with both the pressure rise and decrease stages of the pressure drop oscillation. The third type of oscillation happened at q/G > 1.13 kJ/kg, having the density wave oscillation coupled with the varied liquid film evaporation induced oscillation. The liquid island, retention bubble induced nucleation sites and cone-shape two-phase developing region are unique features of microchannel boiling with pin fin structure. This study illustrated that pressure drop oscillation and density wave oscillation, usually happened in large size channels, also take place in microchannels. The flow pattern transition and varied liquid film evaporation induced oscillations are specific to microchannel boiling/evaporation flow.     

Large-scale synthesis of hierarchical MnO2 for benzene catalytic oxidation

Hierarchical sea-urchin-shaped manganese oxide microspheres were synthesized via a facile method based on the reaction between KMnO₄ and MnSO₄ in HNO₃ solution at 50 °C. The average diameter of the microspheres is ∼850 nm. The microspheres consist of a core of diameter of ∼800 nm and nanorods of width ∼50 nm. The nanorods exist at the edge of the core. The Brunauer–Emmett–Teller surface area of the sea-urchin-shaped microspheres is 259.4 m²/g. A possible formation mechanism of the hierarchical sea-urchin-shaped microspheres is proposed. The temperature for 90% conversion of benzene (T_90%) on the hierarchical urchin-shaped MnO₂ microspheres is about 218 °C.