Unsteady mhd flow and heat transfer on a rotating disk in an ambient fluid

An unsteady flow and heat transfer of a viscous incompressible electrically conducting fluid over a rotating infinite disk in an otherwise ambient fluid are studied. The unsteadiness in the flow field is caused by the angular velocity of the disk which varies with time. The magnetic field is applied normal to the disk surface. The new self-similar solution of the Navier-Stokes and energy equations is obtained numerically. The solution obtained here is not only the solution of the Navier-Stokes equations, but also of the boundary layer equations. Also, for a simple scaling factor, it represents the solution of the flow and heat transfer in the forward stagnation-point region of a rotating sphere or over a rotating cone. The asymptotic behaviour of the solution for a large magnetic field or for a large independent variable is also examined. The surface shear stresses in the radial and tangential directions and the surface heat transfer increase as the acceleration parameter increases. Also the surface shear stress in the radial direction and the surface heat transfer decrease with increasing magnetic field, but the surface shear stress in the tangential direction increases. (C) 2002 Editions scientifiques et medicales Elsevier SAS. All rights reserved.

Effect of Turbine Inlet Temperature on Rotor Blade Tip Leakage Flow and Heat Transfer

One of the most critical gas turbine engine components, the rotor blade tip and casing, is exposed to high thermal load. It becomes a significant design challenge to protect the turbine materials from this severe situation. The purpose of this paper is to study numerically the effect of turbine inlet temperature on the tip leakage flow structure and heat transfer. In this paper, the effect of turbine inlet temperature on the tip leakage flow structure and heat transfer has been studied numerically. Uniform low (LTIT: 444 K) and high (HTIT: 800 K) turbine inlet temperature, as well as non-uniform inlet temperature have been considered. The results showed the higher turbine inlet temperature yields the higher velocity and temperature variations in the leakage flow aerodynamics and heat transfer. For a given turbine geometry and on-design operating conditions, the turbine power output can be increased by 1.33 times, when the turbine inlet temperature increases 1.80 times. Whereas the averaged heat fluxes on the casing and the blade tip become 2.71 and 2.82 times larger, respectively. Therefore, about 2.8 times larger cooling capacity is required to keep the same turbine material temperature. Furthermore, the maximum heat flux on the blade tip of high turbine inlet temperature case reaches up to 3.348 times larger than that of LTIT case. The effect of the interaction of stator and rotor on heat transfer features is also explored using unsteady simulations. The non-uniform turbine inlet temperature enhances the heat flux fluctuation on the blade tip and casing.

Unsteady heat transfer analysis of a film cooling flow

Unsteady heat transfer in a turbine blade film cooling flow is studied using detached eddy simulation (DES). Detailed computation of a single row of 35 degree round holes on a flat plate has been obtained for a blowing ratio of 1.0 and a density ratio of 2.0. The instantaneous flow fields and heat transfer distributions are found to be highly unsteady and oscillatory in nature. The fluctuation of the adiabatic effectiveness and heat transfer coefficient, for example, can be as high as 15 and 50 percent of the time-averaged value, respectively. The correlation between the coherent vortical structures and the unsteady heat transfer is carefully examined. It is shown that the fluctuations in the adiabatic effectiveness and heat transfer coefficient are mainly caused by the spanwise fluctuation of the coolant jet and the thermal turbulent boundary layer accompanying the unsteady flow structures.

Effect of turbine inlet temperature on blade tip leakage flow and heat transfer

One of the most critical gas turbine engine components, rotor blade tip and casing, are exposed to high thermal load. It becomes a significant design challenge to protect the turbine materials from this severe situation. As a result of geometric complexity and experimental limitations, Computational Fluid Dynamics (CFD) tools have been used to predict blade tip leakage flow aerodynamics and heat transfer at typical engine operating conditions. In this paper, the effect of turbine inlet temperature on the tip leakage flow structure and heat transfer has been studied numerically. Uniform low (LTIT: 444 K) and high (HTIT: 800 K) turbine inlet temperature have been considered. The results showed the higher turbine inlet temperature yields the higher velocity and temperature variations in the leakage flow aerodynamics and heat transfer. For a given turbine geometry and on-design operating conditions, the turbine power output can be increased by 1.48 times, when the turbine inlet temperature increases 1.80 times. Whereas the averaged heat fluxes on the casing and the blade tip become 2.71 and 2.82 times larger, respectively. Therefore, about 2.8 times larger cooling capacity is required to keep the same turbine material temperature. Furthermore, the maximum heat flux on the blade tip of high turbine inlet temperature case reaches up to 3.348 times larger than that of LTIT case. The effect of the interaction of stator and rotor on heat transfer features is also explored using unsteady simulations.

Unsteady mixed convection over two-dimensional and axisymmetric bodies

The unsteady laminar incompressible mixed convection flow over a two-dimensional body (cylinder) and an axisymmetric body (sphere) has been studied when the buboyancy forces arise from both thermal and mass diffusion and the unsteadiness in the flow field is introduced by the time dependent free stream velocity. The nonlinear partial differential equations with three independent variables governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The results indicate that for the thermally assisting flow the local skin friction, heat transfer and mass diffusion are enhanced when the buoyancy force from mass diffusion assists the thermal buoyancy force. But this trend is opposite for the thermally opposing flow. The point of zero skin friction moves upstream due to unsteadiness. No singularity is observed at the point of zero skin friction for unsteady flow unlike steady flow. The flow reversal is observed after a certain instant of time. The velocity overshoot occurs for assisting flows.

Flow-acoustic Characterisation of a Cavity-based Combustor Configuration

This study concerns the flow-acoustic characterisation of a cavity-based combustor configuration. A well-validated numerical tool has been used to simulate the unsteady, two-dimensional reacting flow. Initially, a conventional flow over a cavity with dimensions and conditions corresponding to a compact cavity combustor was studied. Cavity mass injections in the form of fuel and air injections required for trapped vortex formation were then employed and the resonance features of this configuration were studied. The results indicate that the cavity depth mode resonance mechanism is dominant at the conditions studied in this work and that the oscillation frequencies do not change with cavity air injection. This observation is important since it implies that the only important variable which can alter resonant frequencies is the cavity depth. With combustion, the pressure oscillation amplitude was observed to increases significantly due to periodic entrainment of the cavity air jet and fluctuation of fuel-air mixture composition to produce highly fluctuating heat-release rates. The underlying mechanisms of the unsteady flow in the cavity combustor identified in this study indicate the strong dependence of the acoustics on the cavity injection strategies.

Two-dimensional flow past circular cylinders using finite volume lattice Boltzmann formulation

In this article, an extension to the total variation diminishing finite volume formulation of the lattice Boltzmann equation method on unstructured meshes was presented. The quadratic least squares procedure is used for the estimation of first-order and second-order spatial gradients of the particle distribution functions. The distribution functions were extrapolated quadratically to the virtual upwind node. The time integration was performed using the fourth-order RungeKutta procedure. A grid convergence study was performed in order to demonstrate the order of accuracy of the present scheme. The formulation was validated for the benchmark two-dimensional, laminar, and unsteady flow past a single circular cylinder. These computations were then investigated for the low Mach number simulations. Further validation was performed for flow past two circular cylinders arranged in tandem and side-by-side. Results of these simulations were extensively compared with the previous numerical data. Copyright (C) 2011 John Wiley & Sons, Ltd.

Unsteady separation and vortex shedding from a laminar separation bubble over a bluff body

Boundary layers are subject to favorable and adverse pressure gradients because of both the temporal and spatial components of the pressure gradient. The adverse pressure gradient may cause the flow to separate. In a closed loop unsteady tunnel we have studied the initiation of separation in unsteady flow past a constriction (bluff body) in a channel. We have proposed two important scalings for the time when boundary layer separates. One is based on the local pressure gradient and the other is a convective time scale based on boundary layer parameters. The flow visualization using a dye injection technique shows the flow structure past the body. Nondimensional shedding frequency (Strouhal number) is calculated based on boundary layer and momentum thicknesses. Strouhal number based on the momentum thickness shows a close agreement with that for flat plate and circular cylinder.

Spectral response of a droplet in pulsating external flow field

A droplet introduced in an external convective flow field exhibits significant multimodal shape oscillations depending upon the intensity of the aerodynamic forcing. In this paper, a theoretical model describing the temporal evolution of normal modes of the droplet shape is developed. The fluid is assumed to be weakly viscous and Newtonian. The convective flow velocity, which is assumed to be incompressible and inviscid, is incorporated in the model through the normal stress condition at the droplet surface and the equation of motion governing the dynamics of each mode is derived. The coupling between the external flow and the droplet is approximated to be a one-way process, i.e., the external flow perturbations effect the droplet shape oscillations and the droplet oscillation itself does not influence the external flow characteristics. The shape oscillations of the droplet with different fluid properties under different unsteady flow fields were simulated. For a pulsatile external flow, the frequency spectra of the normal modes of the droplet revealed a dominant response at the resonant frequency, in addition to the driving frequency and the corresponding harmonics. At driving frequencies sufficiently different from the resonant frequency of the prolate-oblate oscillation mode of the droplet, the oscillations are stable. But at resonance the oscillation amplitude grows in time leading to breakup depending upon the fluid viscosity. A line vortex advecting past the droplet, simulated as an isotropic jump in the far field velocity, leads to the resonant excitation of the droplet shape modes if and only if the time taken by the vortex to cross the droplet is less than the resonant period of the P-2 mode of the droplet. A train of two vortices interacting with the droplet is also analysed. It shows clearly that the time instant of introduction of the second vortex with respect to the droplet shape oscillation cycle is crucial in determining the amplitude of oscillation. (C) 2014 AIP Publishing LLC.

Penetration Characteristics of Air, Carbon Dioxide and Helium Transverse Sonic Jets in Mach 5 Cross Flow

An experimental investigation of sonic air, CO2 and Helium transverse jets in Mach 5 cross flow was carried out over a flat plate. The jet to freestream momentum flux ratio, J, was kept the same for all gases. The unsteady flow topology was examined using high speed schlieren visualisation and PIV. Schlieren visualisation provided information regarding oscillating jet shear layer structures and bow shock, Mach disc and barrel shocks. Two-component PIV measurements at the centreline, provided information regarding jet penetration trajectories. Barrel shocks and Mach disc forming the jet boundary were visualised/quantified also jet penetration boundaries were determined. Even though J is kept the same for all gases, the penetration patterns were found to be remarkably different both at the nearfield and the farfield. Air and CO2 jet resulted similar nearfield and farfield penetration pattern whereas Helium jet spread minimal in the nearfield.

Receptivity To Free-Stream Disturbance Waves For Hypersonic Flow Over A Blunt Cone

A high-order shock-fitting finite difference scheme is studied and used to do direction numerical simulation (DNS) of hypersonic unsteady flow over a blunt cone with fast acoustic waves in the free stream, and the receptivity problem in the blunt cone hypersonic boundary layers is studied. The results show that the acoustic waves are the strongest disturbance in the blunt cone hypersonic boundary layers. The wave modes of disturbance in the blunt cone boundary layers are first, second, and third modes which are generated and propagated downstream along the wall. The results also show that as the frequency decreases, the amplitudes of wave modes of disturbance increase, but there is a critical value. When frequency is over the critial value, the amplitudes decrease. Because of the discontinuity of curvature along the blunt cone body, the maximum amplitudes as a function of frequencies are not monotone.

Interpretation of measured profile losses in unsteady wake-turbine blade interaction studies

The interaction of wakes shed by a moving bladerow with a downstream bladerow causes unsteady flow. The meaning of the freestream stagnation pressure and stagnation enthalpy in these circumstances has been examined using simple analyses, measurements and CFD. The unsteady flow in question arises from the behaviour of the wakes as so-called negative-jets. The interactions of the negative-jets with the downstream blades lead to fluctuations in static pressure which in turn generate fluctuations in the stagnation pressure and stagnation enthalpy. It is shown that the fluctuations of the stagnation quantities created by unsteady effects within the bladerow are far greater than those within the incoming wake. The time-mean exit profiles of the stagnation pressure and stagnation enthalpy are affected by these large fluctuations. This phenomenon of energy separation is much more significant than the distortion of the time-mean exit profiles that is caused directly by the cross-passage transport associated with the negative-jet, as described by Kerrebrock and Mikolajczak. Finally, it is shown that if only time-averaged values of loss are required across a bladerow, it is nevertheless sufficient to determine the time-mean exit stagnation pressure.

Performance of choked unsteady ejector-nozzles for use in pressure-gain combustors

If the conventional steady flow combustor of a gas turbine is replaced with a device which achieves a pressure gain during the combustion process then the thermal efficiency of the cycle is raised. All such 'Pressure Gain Combustors' (e.g. PDEs, pulse combustors or wave rotors) are inherently unsteady flow devices. For such a device to be practically installed in a gas turbine it is necessary to design a downstream row of turbine vanes which will both accept the combustors unsteady exit flow and deliver a flow which the turbine rotor can accept. The design requirements of such a vane are that its exit flow both retains the maximum time-mean stagnation pressure gain (the pressure gain produced by the combustor is not lost) and minimises the amplitude of unsteadiness (reduces unsteadiness entering the downstream rotor). In this paper the exit of the pressure gain combustor is simulated with a cold unsteady jet. The first stage vane is simulated by a one-dimensional choked ejector nozzle with no turning. The time-mean and rms stagnation pressure at nozzle exit is measured. A number of geometric configurations are investigated and it is shown that the optimal geometry both maximizes time mean stagnation pressure gain (75% of that in the exit of the unsteady jet) and minimizes the amplitude of unsteadiness (1/3 of that in the primary jet). The structure of the unsteady flow within the ejector nozzle is determined computationally. Copyright © 2009 by J Heffer and R Miller.

Airfoil flow and noise computation using monotonically integrated large eddy simulation and acoustic analogy: Effect of the grid resolution.

A new scalable Monotonically Integrated Large Eddy Simulation (MILES) method based on the Compact Accurately Boundary-Adjusting high-REsolution Technique (CABARET) has been applied for the simulation of unsteady flow around NACA0012 airfoil at Re = 400,000 and M = 0.058. The flow solution is coupled with the Ffowcs Williams-Hawkings formulation for far-field noise prediction. The computational modeling results are presented for several computational grid resolutions: 8, 16, and 32 million grid cells and compared with the experimental data available.

Numerical simulation of flow-induced cavity noise in self-sustained oscillations

The generation and near-field radiation of aerodynamic sound from a low-speed unsteady flow over a two-dimensional automobile door cavity is simulated by using a source-extraction-based coupling method. In the coupling procedure, the unsteady cavity flow field is first computed solving the Reynolds averaged Navier–Stokes (RANS) equations. The radiated sound is then calculated by using a set of acoustic perturbation equations with acoustic source terms which are extracted from the time-dependent solutions of the unsteady flow. The aerodynamic and its resulting acoustic field are computed for the Reynolds number of 53,266 based on the base length of the cavity. The free stream flow velocity is taken to be 50.9m/s. As first stage of the numerical investigation of flow-induced cavity noise, laminar flow is assumed. The CFD solver is based on a cell-centered finite volume method. A dispersion-relation-preserving (DRP), optimized, fourth-order finite difference scheme with fully staggered-grid implementation is used in the acoustic solver