Visualization of turbulent wedges under favourable pressure gradients using both shear-sensitive and temperature-sensitive liquid crystals

Turbulent wedges induced by a 3D surface roughness placed in a laminar boundary layer over a flat plate were visualised for the first time using both shear-sensitive and temperature-sensitive liquid crystals. The experiments were carried out at three different levels of favourable pressure gradients. The purpose of this investigation was to examine the spreading angles of the turbulent wedges indicated by their associated surface shear stresses and heat transfer characteristics and hence obtain further insight about the difference in the behaviour of transitional momentum and thermal boundary layers when a streamwise pressure gradient exists. It was shown that under a zero pressure gradient the spreading angles indicated by the two types of liquid crystals are the same, but the difference increases as the level of favourable pressure gradient increases. The result from the present study could have an important implication to the transition modelling of thermal boundary layers over gas turbine blades.

The development of aerodynamic performance measurements in a transient test facility

Transient test facilities offer the potential for the simultaneous study of turbine aerodynamic performance, unsteady flow phenomena and the heat transfer characteristics of a turbine stage. This paper describes the development of aerodynamic performance measurement techniques in the Oxford Rotor Facility (ORF). The solutions to the technological issues involved with transient testing presented in this paper are expected to achieve levels of precision uncertainty comparable with traditional steady flow test rigs. The theoretical background to the measurement of aerodynamic performance is presented together with a comprehensive pre-test uncertainty analysis. The instrumentation scheme for the measurement of stage mass flow rate is discussed in detail, the measurements of shaft power, total inlet enthalpy, and stage pressure ratio are also outlined. The current working section features a 62% scale, 1-1/2 stage, high-pressure shroudless transonic turbine. The required inlet flow conditions are provided by an Isentropic Light Piston Tunnel (ILPT) with a quasi-steady state run time of approximately 70ms. The testing is conducted at engine representative specific speed, pressure ratio, gas-to-wall temperature ratio, Mach number and Reynolds number.

Deposition of small particles from turbulent flows

This paper deals with particle deposition onto solid walls from turbulent flows. The aim of the study is to model particle deposition in industrial flows, such as the one in gas turbines. The numerical study has been carried out with a two fluid approach. The possible contribution to the deposition from Brownian diffusion, turbulent diffusion and shear-induced lift force are considered in the study. Three types of turbulent two-phase flows have been studied: turbulent channel flow, turbulent flow in a bent duct and turbulent flow in a turbine blade cascade. In the turbulent channel flow case, the numerical results from a two-dimensional code show good agreement with numerical and experimental results from other resources. Deposition problem in a bent duct flow is introduced to study the effect of curvature. Finally, the deposition of small particles on a cascade of turbine blades is simulated. The results show that the current two fluid models are capable of predicting particle deposition rates in complex industrial flows.

Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

Numerical simulation of the growth of ZnO nanostructures in a tube furnace by physical vapour deposition

Zinc oxide is a versatile II-VI naturally n-type semiconductor that exhibits piezoelectric properties. By controlling the growth kinetics during a simple carbothermal reduction process a wide range of 1D nanostructures such as nanowires, nanobelts, and nanotetrapods have been synthesized. The driving force: for the nanostructure growth is the Zn vapour supersaturation and supply rate which, if known, can be used to predict and explain the type of crystal structure that results. A model which attempts to determine the Zn vapour concentration as a function of position in the growth furnace is described. A numerical simulation package, COMSOL, was used to simultaneously model the effects of fluid flow, diffusion and heat transfer in a tube furnace made specifically for ZnO nanostructure growth. Parameters such as the temperature, pressure, and flow rate are used as inputs to the model to show the effect that each one has on the Zn concentration profile. An experimental parametric study of ZnO nanostructure growth was also conducted and compared to the model predictions for the Zn concentration in the tube. © 2008 Materials Research Society.

CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

Physics in fish processing technology

To design, develop and put into operation an equipment either to increase the productivity or to improve the existing technique to obtain a better quality of the product, the fishery engineer/scientist should have a comprehensive knowledge of fundamental principles involved in the process. Many a technique in fish processing technology, whether it applies to freezing, dehydration or canning, involves always a type of heat transfer, which is dependent to a certain extent on the external physical parameters like temperature. humidity, pressure, air flow etc. and also on the thermodynamic properties of fish muscle in the temperature ranges encountered. Similarly informations on other physical values like dielectric constant and dielectric loss in the design of quick trawlers and in quality assessment of frozen/iced fish, refractive index and viscosity in the measurement of the saturation and polymerisation of fish oils and shear strength in the judgement of textural qualities of cooked fish are also equally important.

Influence of a building's integrated-photovoltaics on heating and cooling loads

Building integrated photovoltaics (BIPV) has the potential to become a major source of renewable energy in the urban environment. BIPV has significant influence on the heat transfer through the building envelope because of the change of the thermal resistance by adding or replacing the building elements. Four different roofs are used to assess the impacts of BIPV on the building's heating-and-cooling loads; namely ventilated air-gap BIPV, non-ventilated (closed) air-gap BIPV, closeroof mounted BIPV, and the conventional roof with no PV and no air gap. One-dimensional transient models of four cases are derived to evaluate the PV performances and building cooling-and-heating loads across the different roofs in order to select the appropriate PV building integration method in Tianjin, China. The simulation results show that the PV roof with ventilated air-gap is suitable for the application in summer because this integration leads to the low cooling load and high PV conversion efficiency. The PV roof with ventilation air-gap has a high time lag and small decrement factor in comparison with other three roofs and has the same heat gain as the cool roof of absorptance 0.4. In winter, BIPV of non-ventilated air gap is more appropriate due to the combination of the low heating-load through the PV roof and high PV electrical output. © 2005 Elsevier Ltd. All rights reserved.

Electricity Storage Using a Thermal Storage Scheme

The increasing use of renewable energy technologies for electricity generation, many of which have an unpredictably intermittent nature, will inevitably lead to a greater demand for large-scale electricity storage schemes. For example, the expanding fraction of electricity produced by wind turbines will require either backup or storage capacity to cover extended periods of wind lull. This paper describes a recently proposed storage scheme, referred to here as Pumped Thermal Storage (PTS), and which is based on "sensible heat" storage in large thermal reservoirs. During the charging phase, the system effectively operates as a high temperature-ratio heat pump, extracting heat from a cold reservoir and delivering heat to a hot one. In the discharge phase the processes are reversed and it operates as a heat engine. The round- trip efficiency is limited only by process irreversibilities (as opposed to Second Law limitations on the coefficient of performance and the thermal efficiency of the heat pump and heat engine respectively). PTS is currently being developed in both France and England. In both cases, the schemes operate on the Joule-Brayton (gas turbine) cycle, using argon as the working fluid. However, the French scheme proposes the use of turbomachinery for compression and expansion, whereas for that being developed in England reciprocating devices are proposed. The current paper focuses on the impact of the various process irreversibilities on the thermodynamic round-trip efficiency of the scheme. Consideration is given to compression and expansion losses and pressure losses (in pipe-work, valves and thermal reservoirs); heat transfer related irreversibility in the thermal reservoirs is discussed but not included in the analysis. Results are presented demonstrating how the various loss parameters and operating conditions influence the overall performance.

The mechanical integrity of fuel pin cladding in a pulsed-beam accelerator driven subcritical reactor

The Accelerator Driven Subcritical Reactor (ADSR) is one of the reactor designs proposed for future nuclear energy production. Interest in the ADSR arises from its enhanced and intrinsic safety characteristics, as well as its potential ability to utilize the large global reserves of thorium and to burn legacy actinide waste from other reactors and decommissioned nuclear weapons. The ADSR concept is based on the coupling of a particle accelerator and a subcritical core by means of a neutron spallation target interface. One of the candidate accelerator technologies receiving increasing attention, the Fixed Field Alternating Gradient (FFAG) accelerator, generates a pulsed proton beam. This paper investigates the impact of pulsed proton beam operation on the mechanical integrity of the fuel pin cladding. A pulsed beam induces repetitive temperature changes in the reactor core which lead to cyclic thermal stresses in the cladding. To perform the thermal analysis aspects of this study a code that couples the neutron kinetics of a subcritical core to a cylindrical geometry heat transfer model was developed. This code, named PTS-ADS, enables temperature variations in the cladding to be calculated. These results are then used to perform thermal fatigue analysis and to predict the stress-life behaviour of the cladding. © 2011 Elsevier Ltd. All rights reserved.

Increased gas production from hydrates by combining depressurization with heating of the wellbore

To extract gas from hydrate reservoirs, it has to be dissociated in situ. This endothermic dissociation process absorbs heat energy from the formation and pore fluid. The heat transfer governs the dissociation rate, which is proportional to the difference between the actual temperature and the equilibrium temperature. This study compares three potential gas production schemes from hydrate-bearing soil, where the radial heat transfer is governing. Cylindrical samples with 40% pore-filling hydrate saturation were tested. The production tests were carried out over 90 min by dissociating the hydrate from a centered miniature wellbore, by either lowering the pressure to 6, 4, or 6 MPa with simultaneous heating of the wellbore to 288 K. All tests were replicated by a numerical simulation. With additional heating at the same wellbore pressure, the gas production from hydrates could, on average, be increased by 1.8 and 3.6 times in the simulation and experiments, respectively. If the heat influx from the outer boundary is limited, a simulation showed that the specific heat of the formation is rapidly used up when the wellbore is only depressurized and not heated. © 2012 American Chemical Society.

Strategies for modeling turbulent flows in electronics

Hybrid methods based on the Reynolds Averaged Navier Stokes (RANS) equations and the Large Eddy Simulation (LES) formulation are investigated to try and improve the accuracy of heat transfer and surface temperature predictions for electronics systems and components. Two relatively low Reynolds number flows are studied using hybrid RANS-LES, RANS-Implicit-LES (RANS-ILES) and non-linear LES models. Predictions using these methods are in good agreement with each other, even using different grid resolutions. © 2008 IEEE.

Predictive techniques for turbulent oscillatory flows

The prediction of turbulent oscillatory flow at around transitional Reynolds numbers is considered for an idealized electronics system. To assess the accuracy of turbulence models, comparison is made with measurements. A stochastic procedure is used to recover instantaneous velocity time traces from predictions. This procedure enables more direct comparison with turbulence intensity measurements which have not been filtered to remove the oscillatory flow component. Normal wall distances, required in some turbulence models, are evaluated using a modified Poisson equation based technique. A range of zero, one and two equation turbulence models are tested, including zonal and a non-linear eddy viscosity models. The non-linear and zonal models showed potential for accuracy improvements.

Computation of transient and steady turbulent isothermal flows with an impinging rectangular jet

The computation of both transient and steady turbulent incompressible isothermal flows is studied. The flow is very complex, having streamline curvature, large vortex structures and stagnation resulting from an impinging rectangular jet. For transient computations, the standard k-ε model is adopted. For steady flows, the k-ε, high and low Reynolds number k-l and mixing length models are tried. Zonal approaches combining the above turbulence models are also investigated. None of the models are found to give satisfactory agreement with velocity measurements.

Numerical computation of laminar flow in a heated rotating cavity with an axial throughflow of air

A heated rotating cavity with an axial throughflow of cooling air is used as a model for the flow in the cylindrical cavities between adjacent discs of a high-pressure gas-turbine compressor. In an engine the flow is expected to be turbulent, the limitations of this laminar study are fully realised but it is considered an essential step to understand the fundamental nature of the flow. The three-dimensional, time-dependent governing equations are solved using a code based on the finite volume technique and a multigrid algorithm. The computed flow structure shows that flow enters the cavity in one or more radial arms and then forms regions of cyclonic and anticyclonic circulation. This basic flow structure is consistent with existing experimental evidence obtained from flow visualization. The flow structure also undergoes cyclic changes with time. For example, a single radial arm, and pair of recirculation regions can commute to two radial arms and two pairs of recirculation regions and then revert back to one. The flow structure inside the cavity is found to be heavily influenced by the radial distribution of surface temperature imposed on the discs. As the radial location of the maximum disc temperature moves radially outward, this appears to increase the number of radial arms and pairs of recirculation regions (from one to three for the distributions considered here). If the peripheral shroud is also heated there appear to be many radial arms which exchange fluid with a strong cyclonic flow adjacent to the shroud. One surface temperature distribution is studied in detail and profiles of the relative tangential and radial velocities are presented. The disc heat transfer is also found to be influenced by the disc surface temperature distribution. It is also found that the computed Nusselt numbers are in reasonable accord over most of the disc surface with a correlation found from previous experimental measurements. © 1994, MCB UP Limited.