Mechanical Engineering - Research Publications
Permanent URI for this collection
Search Results
1 - 10 of 10
-
ItemNo Preview AvailableModelling the downstream development of a turbulent boundary layer following a step change of roughness(CAMBRIDGE UNIV PRESS, 2022-09-23)In this study, we develop an analytical model to predict the turbulent boundary layer downstream of a step-change in the surface roughness where upstream flow conditions are given. We first revisit the classical model of Elliott (Trans. Am. Geophys. Union, vol. 39, 1958, pp. 1048–1054), who modelled the velocity distribution within and above the internal layer with a simple piecewise logarithmic profile, and evolved the velocity profile using the streamwise momentum equation. Elliott's model was originally developed for an atmospheric surface layer, and to make the model applicable to a spatially developing turbulent boundary layer with finite thickness, we propose a number of more physical refinements, including adding a wake function to the velocity profile, considering the growth of the entire boundary layer in the streamwise direction, and using a more realistic shear stress profile in the momentum equation. In particular, we implement the blending model (Li et al., J. Fluid Mech., vol. 923, 2021, p. A18) to account for the deviation of the mean flow within the internal layer from a canonical velocity profile based on the local wall condition. These refinements lead to improved agreement between the prediction and the measurement, especially in the vicinity of the rough-to-smooth change.
-
ItemNo Preview AvailableNavier-Stokes-based linear model for unstably stratified turbulent channel flows(AMER PHYSICAL SOC, 2022-04-06)
-
ItemAn energy-efficient pathway to turbulent drag reduction(NATURE PORTFOLIO, 2021-10-04)Simulations and experiments at low Reynolds numbers have suggested that skin-friction drag generated by turbulent fluid flow over a surface can be decreased by oscillatory motion in the surface, with the amount of drag reduction predicted to decline with increasing Reynolds number. Here, we report direct measurements of substantial drag reduction achieved by using spanwise surface oscillations at high friction Reynolds numbers ([Formula: see text]) up to 12,800. The drag reduction occurs via two distinct physical pathways. The first pathway, as studied previously, involves actuating the surface at frequencies comparable to those of the small-scale eddies that dominate turbulence near the surface. We show that this strategy leads to drag reduction levels up to 25% at [Formula: see text] = 6,000, but with a power cost that exceeds any drag-reduction savings. The second pathway is new, and it involves actuation at frequencies comparable to those of the large-scale eddies farther from the surface. This alternate pathway produces drag reduction of 13% at [Formula: see text] = 12,800. It requires significantly less power and the drag reduction grows with Reynolds number, thereby opening up potential new avenues for reducing fuel consumption by transport vehicles and increasing power generation by wind turbines.
-
ItemHeat Transfer Coefficient Estimation for Turbulent Boundary Layers(The University of Queensland, 2020-12-11)Convective heat transfer in rough wall-bounded turbulent flows is prevalent in many engineering applications, such as in gas turbines and heat exchangers. At present, engineers lack the design tools to accurately predict the convective heat transfer in the presence of non-smooth boundaries. Accordingly, a new turbulent boundary layer facility has been commissioned, where the temperature of an interchangeable test surface can be precisely controlled, and conductive heat losses are minimized. Using this facility, we can estimate the heat transfer coefficient (Stanton number, St), through measurement of the power supplied to the electrical heaters and also from measurements of the thermal and momentum boundary layers evolving over this surface. These methods have been initially investigated over a shorter smooth prototype heated surface and compared with existing St prediction models. Preliminary results suggest that we can accurately estimate St in this facility.
-
ItemAn investigation of cold-wire spatial resolution using a DNS database(The University of Queensland, 2020-12-11)The effect of spatial resolution of cold-wire anemometry on both the variance and energy spectrum of temperature fluctuations is analyzed through the use of a numerical database. Temperature fluctuation snapshots from a direct numerical simulation (DNS) of a heated smooth-wall turbulent channel flow are spatially averaged in the spanwise direction to simulate the wire filtering. The results show that the wire length does not affect the mean temperature while it significantly attenuates the variance of temperature fluctuations, particularly in the vicinity of the wall. As the filter length grows, the peaks of the one- and two-dimensional energy spectrograms are further attenuated. Limited attenuation is seen when the filter length is smaller than 30 wall units in the vicinity of the wall, whereas a complete suppression of the near-wall energetic peak is observed when the filter length exceeds 100 wall units.
-
Item
-
ItemCoriolis effect on centrifugal buoyancy-driven convection in a thin cylindrical shell(Cambridge University Press (CUP), 2021-03-10)We study the effect of the Coriolis force on centrifugal buoyancy-driven convection in a rotating cylindrical shell with inner cold wall and outer hot wall. This is done by performing direct numerical simulations for increasing inverse Rossby number Ro−1 from zero (no Coriolis force) to 20 (very large Coriolis force) and for Rayleigh number Ra from 107 to 1010 and Prandtl number Pr=0.7, corresponding to air. We invoke the thin-shell limit, which neglects the curvature and radial variations of the centripetal acceleration. As Ro−1 increases from zero, the system forms an azimuthal bidirectional wind that reaches its maximum momentum at an optimal Ro−1opt, associated with a maximal skin-friction coefficient Cf and a minimal Nusselt number Nu. Just beyond Ro−1opt, the wind weakens and an axial, quasi-two-dimensional cyclone, corotating with the system, begins to form. A local ‘turbulence’ inverse Rossby number (non-dimensionalised by the eddy turnover time) determines the onset of cyclone formation for all Ra, when its value reaches approximately 4. At Ro−1≫Ro−1opt, the system falls into the geostrophic regime with a sudden drop in Nu. The bidirectional wind for Ro−1≤Ro−1opt is a feature of this system, as it hastens the boundary layer transition from laminar to turbulent, towards the ultimate regime. We see the onset of this transition at Ra=1010 and Ro−1≃Ro−1opt, although the mean flow profile has not yet fully collapsed on the Prandtl–von Kármán (logarithmic) law.
-
ItemThe effect of spanwise wavelength of surface heterogeneity on turbulent secondary flows(Cambridge University Press (CUP), 2020-07-10)We examine the behaviour of turbulent boundary layers over surfaces composed of spanwise-alternating smooth and rough strips, where the width of the strips varies such that, where is the boundary-layer thickness averaged over one spanwise wavelength of the heterogeneity. The experiments are configured to examine the influences of spanwise variation in wall shear stress over a large range. Hot-wire anemometry and particle image velocimetry (PIV) reveal that the half-wavelength governs the diameter and strength of the resulting mean secondary flows and hence the observed isovels of the mean streamwise velocity. Three possible cases are observed: limiting cases (either or), where the secondary flows are confined near the wall or near the roughness change, and intermediate cases (), where the secondary flows are space filling and at their strongest. These secondary flows, however, exhibit a time-dependent behaviour which might be masked by time averaging. Further analysis of the energy spectrogram and fluctuating flow fields obtained from PIV show that the secondary flows meander in a similar manner to that of large-scale structures occurring naturally in turbulence over smooth walls. The meandering of the secondary flows is a function of and is most prominent when.
-
ItemGlobal and local aspects of entrainment in temporal plumes(CAMBRIDGE UNIV PRESS, 2017-02-10)To date, the understanding of the role buoyancy plays in the entrainment process in unstable configurations such as turbulent plumes remains incomplete. Towards addressing this question, we set up a flow in which a plume evolves in time instead of space. We demonstrate that the temporal problem is equivalent to a spatial plume in a strong coflow and address in detail how the temporal plume can be realized via direct numerical simulation. Using numerical data of plume simulations up to $Re_{\unicode[STIX]{x1D706}}\approx 100$, we show that the entrainment coefficient can be determined consistently using a global entrainment analysis in an integral framework as well as via a local approach. The latter is based on a study of the local propagation of the turbulent/non-turbulent interface relative to the fluid. Locally, this process is dominated by small-scale diffusion which is amplified by interface convolutions such that the total entrained flux is independent of viscosity. Further, we identify a direct buoyancy contribution to entrainment by baroclinic torque, which accounts for 8 %–12 % of the entrained flux locally, comparable to the 15 % buoyancy contribution at the integral level. It appears that the baroclinic torque is a mechanism that might explain higher values of the entrainment coefficient in spatial plumes compared with jets.
-
ItemCharacteristics of the entrainment velocity in a developing wake(International Symposium on Turbulence and Shear Flow Phenomena, 2015)