Instabilities in the oscillatory flow of a complex fluid.

The dynamics of a fluid in a vertical tube, subjected to an oscillatory pressure gradient, is studied experimentally for both a Newtonian and a viscoelastic shear-thinning fluid. Particle image velocimetry is used to determine the two-dimensional velocity fields in the vertical plane of the tube axis, in a range of driving amplitudes from 0.8 to 2.5 mm and of driving frequencies from 2.0 to 11.5 Hz. The Newtonian fluid exhibits a laminar flow regime, independent of the axial position, in the whole range of drivings. For the complex fluid, instead, the parallel shear flow regime exhibited at low amplitudes [Torralba, Phys. Rev. E 72, 016308 (2005)] becomes unstable at higher drivings against the formation of symmetric vortices, equally spaced along the tube. At even higher drivings the vortex structure itself becomes unstable, and complex nonsymmetric structures develop. Given that inertial effects remain negligible even at the hardest drivings (Re < 10(-1)), it is the complex rheology of the fluid that is responsible for the instabilities observed. The system studied represents an interesting example of the development of shear-induced instabilities in nonlinear complex fluids in purely parallel shear flow.

Oscillatory sensitivity patterns for global modes in wakes

Globally unstable wakes with co-flow at intermediate Reynolds numbers are studied, to quantify important spatial regions for the development and control of the global instability. One region of high structural sensitivity is found close to the inlet for all wakes, in agreement with previous findings for cylinder wakes. A second, elongated region of high structural sensitivity is seen downstream of the first one for unconfined wakes at Re = 400. When base flow modifications are considered, a spatially oscillating sensitivity pattern is found inside the downstream high structural sensitivity region. This implies that the same change in the base flow can either destabilize or stabilize the flow, depending on the exact position where it is applied. It is shown that the sensitivity pattern remains unchanged for different choices of streamwise boundary conditions and numerical resolution. Actual base flow modifications are applied in selected configurations, and the linear global modes recomputed. It is confirmed that the linear global eigenvalues move according to the predicted sensitivity pattern for small amplitude base flow modifications, for which the theory applies. We also look at the implications of a small control cylinder on the flow. Only the upstream high sensitivity region proves to be robust in terms of control, but one should be careful not to disturb the flow in the downstream high sensitivity region, in order to achieve control. The findings can have direct implications on the numerical resolution requirements for wakes at higher Reynolds numbers. Furthermore, they provide one more possible explanation to why confined wakes have a more narrow frequency spectrum than unconfined wakes.

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.

Nonlinear thermoacoustics of ducted premixed flames: The influence of perturbation convection speed

When a premixed flame is placed within a duct, acoustic waves induce velocity perturbations at the flame's base. These travel down the flame, distorting its surface and modulating its heat release. This can induce self-sustained thermoacoustic oscillations. Although the phase speed of these perturbations is often assumed to equal the mean flow speed, experiments conducted in other studies and Direct Numerical Simulation (DNS) conducted in this study show that it varies with the acoustic frequency. In this paper, we examine how these variations affect the nonlinear thermoacoustic behaviour. We model the heat release with a nonlinear kinematic G-equation, in which the velocity perturbation is modelled on DNS results. The acoustics are governed by linearised momentum and energy equations. We calculate the flame describing function (FDF) using harmonic forcing at several frequencies and amplitudes. Then we calculate thermoacoustic limit cycles and explain their existence and stability by examining the amplitude-dependence of the gain and phase of the FDF. We find that, when the phase speed equals the mean flow speed, the system has only one stable state. When the phase speed does not equal the mean flow speed, however, the system supports multiple limit cycles because the phase of the FDF changes significantly with oscillation amplitude. This shows that the phase speed of velocity perturbations has a strong influence on the nonlinear thermoacoustic behaviour of ducted premixed flames. © 2013 The Combustion Institute.

An electrokinetic sensor for studying immersed surfaces, using focused ultrasound

A novel method for obtaining information on the charge density of an immersed surface is presented. The technique uses focused ultrasound to excite oscillatory fluid motion in the plane of the solid-liquid interface, over a localised area. The displacement current (resulting from the motion of fluid-borne ions in the outer double-layer) is detected by electrodes in the liquid. The method is demonstrated as a means for monitoring protein adsorption, and for monitoring interactions between two different proteins. A second electrokinetic effect at the interface is identified, isolated from the first, and shown to provide additional information on the compressibility and charge density of the double-layer. © 2001 Elsevier Science B.V. All rights reserved.

Transient Natural Ventilation of a Room with a Distributed Heat Source

We report on an experimental and theoretical study of the transient flows which develop as a naturally ventilated room adjusts from one temperature to another. We focus on a room heated from below by a uniform heat source, with both high- and low-level ventilation openings. Depending on the initial temperature of the room relative to (i) the final equilibrium temperature and (ii) the exterior temperature, three different modes of ventilation may develop. First, if the room temperature lies between the exterior and the equilibrium temperature, the interior remains well-mixed and gradually heats up to the equilibrium temperature. Secondly, if the room is initially warmer than the equilibrium temperature, then a thermal stratification develops in which the upper layer of originally hot air is displaced upwards by a lower layer of relatively cool inflowing air. At the interface, some mixing occurs owing to the effects of penetrative convection. Thirdly, if the room is initially cooler than the exterior, then on opening the vents, the original air is displaced downwards and a layer of ambient air deepens from above. As this lower layer drains, it is eventually heated to the ambient temperature, and is then able to mix into the overlying layer of external air, and the room becomes well-mixed. For each case, we present new laboratory experiments and compare these with some new quantitative models of the transient flows. We conclude by considering the implications of our work for natural ventilation of large auditoria.

Continuous spectrum growth and modal instability in swirling duct flow

We present results on the stability of compressible inviscid swirling flows in an annular duct. Such flows are present in aeroengines, for example in the by-pass duct, and there are also similar flows in many aeroacoustic or aeronautical applications. The linearised Euler equations have a ('critical layer') singularity associated with pure convection of the unsteady disturbance by the mean flow, and we focus our attention on this region of the spectrum. By considering the critical layer singularity, we identify the continuous spectrum of the problem and describe how it contributes to the unsteady field. We find a very generic family of instability modes near to the continuous spectrum, whose eigenvalue wavenumbers form an infinite set and accumulate to a point in the complex plane. We study this accumulation process asymptotically, and find conditions on the flow to support such instabilities. It is also found that the continuous spectrum can cause a new type of instability, leading to algebraic growth with an exponent determined by the mean flow, given in the analysis. The exponent of algebraic growth can be arbitrarily large. Numerical demonstrations of the continuous spectrum instability, and also the modal instabilities are presented.

Thermal transport in high porosity cellular metal foams

The heat dissipation capability of highly porous cellular metal foams with open cells subject to forced air convection is studied using a combined experimental and analytical approach. The cellular morphologies of six FeCrAlY (an iron-based alloy) foams and six copper alloy foams with a range of pore sizes and porosities are quantified with the scanning electronic microscope and image analysis. Experimental measurements on pressure drop and heat transfer for copper foams are carried out. A numerical model for forced convection across open-celled metal foams is subsequently developed, and the predictions are compared with those measured. Reasonably good agreement with test data is obtained, given the complexity of the cellular foam morphology and the associated momentum/energy transport. The results show that cell size has a more significant effect on the overall heat transfer than porosity. An optimal porosity is obtained based on the balance between pressure drop and overall heat transfer, which decreases as the Reynolds number is increased.

Convective heat dissipation with lattice-frame materials

This paper presents experimental results on heat transfer and pressure drop for a compact heat sink made of fully triangulated, lightweight (porosity∼0.938), aluminum lattice-frame materials (LFMs). Due to the inherent structural anisotropy of the LFMs, two mutually perpendicular orientations were selected for the measurements. Constant heat flux was applied to the heat sink under steady state conditions, and dissipated by forced air convection. The experimental data were compared with those predicted from an analytical model based on fin analogy. The experimental results revealed that pressure drop is strongly dependent upon the orientation of the structure, due mainly to the flow blockage effect. For heat transfer measurements, typical local temperature distributions on the substrate under constant heat flux conditions were captured with infrared camera. The thermal behavior of LFMs was found to follow closely that of cylinder banks, with early transition Reynolds number (based on strut diameter) equal to about 300. The Nusselt number prediction from the fin-analogy correlates well with experimental measurements, except at low Reynolds numbers where a slightly underestimation is observed. Comparisons with empty channels and commonly used heat exchanger media show that the present LFM heat sink can remove heat approximately seven times more efficient than an empty channel and as efficient as a bank of cylinders at the same porosity level. The aluminum LFMs are extremely stiff and strong, making them ideal candidates for multifunctional structures requiring both heat dissipation and mechanical load carrying capabilities. © 2003 Elsevier Ltd. All rights reserved.

The temperature dependence of effective thermal conductivity of open-celled steel alloy foams

The effective thermal conductivity of steel alloy FeCrAlY (Fe-20 wt.% Cr-5 wt.% Al-2 wt.% Y-20 wt.%) foams with a range of pore sizes and porosities was measured between 300 and 800 K, under both vacuum and atmospheric conditions. The results show that the effective thermal conductivity increases rapidly as temperature is increased, particularly in the higher temperature range (500-800 K) where the transport of heat is dominated by thermal radiation. The effective conductivity at temperature 800 K can be three times higher than that at room temperature (300 K). Results obtained under vacuum conditions reveal that the effective conductivity increases with increasing pore size or decreasing porosity. The contribution of natural convection to heat conduction was found to be significant, with the effective thermal conductivity at ambient pressure twice the value of vacuum condition. The results also show that natural convection in metal foams is strongly dependent upon porosity. © 2003 Elsevier B.V. All rights reserved.