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.
-
ItemReorganisation of turbulence by large and spanwise-varying riblets(Cambridge University Press, 2022-12-10)We study the flow above non-optimal riblets, specifically large drag-increasing and two-scale trapezoidal riblets. In order to reach large Reynolds numbers and large scale separation while retaining access to flow details, we employ a combination of boundary-layer hot-wire measurements and direct numerical simulation (DNS) in minimal-span channels. Although the outer Reynolds numbers differ, we observe fair agreement between experiments and DNS at matched viscous-friction-scaled riblet spacings in the overlapping physical and spectral regions, providing confidence that both data sets are valid. We find that hot-wire velocity spectra above very large riblets with are depleted of near-wall energy at scales that are (much) greater than. Large-scale energy likely bypasses the turbulence cascade and is transferred directly to secondary flows of size, which we observe to grow in strength with increasing riblet size. Furthermore, the present very large riblets reduce the von Kármán constant of the spanwise uniform mean velocity in a logarithmic layer and, thus, reduce the accuracy of the roughness-function concept, which we link to the near-wall damping of large flow structures. Half-height riblets in the groove, which we use as a model of imperfectly repeated (spanwise-varying) riblets, impede in-groove turbulence. We show how to scale the drag optimum of imperfectly repeated riblets based on representative measurements of the true geometry by solving inexpensive Poisson equations.
-
ItemRiblet-generated flow mechanisms that lead to local breaking of Reynolds analogy(CAMBRIDGE UNIV PRESS, 2022-11-14)We investigate the Reynolds analogy over riblets, namely the analogy between the fractional increase in Stanton number Ch and the fractional increase in the skin-friction coefficient Cf, relative to a smooth surface. We investigate the direct numerical simulation data of Endrikat et al. (Flow Turbul. Combust., vol. 107, 2021, pp. 1–29). The riblet groove shapes are isosceles triangles with tip angles α = 30◦, 60◦, 90◦, a trapezoid, a rectangle and a right triangle. The viscous-scaled riblet spacing varies between s+ ≈ 10 to 60. The global Reynolds analogy is primarily influenced by Kelvin–Helmholtz rollers and secondary flows. Kelvin–Helmholtz rollers locally break the Reynolds analogy favourably, i.e., cause a locally larger fractional increase in Ch than in Cf. These rollers induce negative wall shear stress patches which have no analogue in wall heat fluxes. Secondary flows at the riblets’ crests are associated with local unfavourable breaking of the Reynolds analogy, i.e., locally larger fractional increase in Cf than in Ch. Only the triangular riblets with α = 30◦ trigger strong Kelvin–Helmholtz rollers without appreciable secondary flows. This riblet shape globally preserves the Reynolds analogy from s+ = 21 to 33. However, the other riblet shapes have weak or non-existent Kelvin–Helmholtz rollers, yet persistent secondary flows. These riblet shapes behave similarly to rough surfaces. They unfavourably break the global Reynolds analogy and do so to a greater extent as s+ increases.
-
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.
-
ItemDispersive stresses in turbulent flow over riblets(Cambridge University Press, 2021-06-25)We carry out direct numerical simulations of turbulent flow over riblets, streamwise- aligned grooves that are designed to reduce drag by modifying the near-wall flow. Twenty riblet geometries and sizes are considered, namely symmetric triangular with tip angle, and, asymmetric triangular, blade and trapezoidal. To save on computational cost, simulations are performed using the minimal-channel flow configuration. With this unprecedented breadth of high-fidelity flow data near the wall, we are able to obtain more general insights into the flow physics of riblets. As observed by García-Mayoral & Jiménez (J. Fluid Mech., vol. 678, 2011, pp. 317-347), we confirm that the drag-change curves of all the present groove geometries better collapse when reported with the viscous-scaled square root of the groove area, rather than the riblet spacing. Using a two-dimensional generalization of the Fukagata-Iwamoto-Kasagi identity in difference form we isolate the different drag-change contributions. We show that the drag increase associated with dispersive stresses carried by secondary flows can be as important as the one associated with the turbulent stresses and the pre-eminence of dispersive stresses can be estimated by the groove width at the riblet mean height.
-
ItemInfluence of riblet shapes on the occurrence of Kelvin-Helmholtz rollers(CAMBRIDGE UNIV PRESS, 2021-03-02)We investigate turbulent flow over streamwise-aligned riblets (grooves) of various shapes and sizes. Small riblets with spacings of typically less than 20 viscous units are known to reduce skin-friction drag compared to a smooth wall, but larger riblets allow inertial-flow mechanisms to appear and cause drag reduction to break down. One of these mechanisms is a Kelvin–Helmholtz instability that García-Mayoral & Jiménez (J. Fluid Mech., vol. 678, 2011, pp. 317–347) identified in turbulent flow over blade riblets. In order to evaluate its dependence on riblet shape and thus gain a broader understanding of the underlying physics, we generate an extensive data set comprising 21 cases using direct numerical simulations of fully developed minimal-span channel flow. The data set contains six riblet shapes of varying sizes between maximum drag reduction and significant drag increase. Comparing the flow fields over riblets to that over a smooth wall, we find that in this data set only large sharp-triangular and blade riblets have a drag penalty associated with the Kelvin–Helmholtz instability and that the mechanism appears to be absent for blunt-triangular and trapezoidal riblets of any size. We therefore investigate two indicators for the occurrence of Kelvin–Helmholtz rollers in turbulent flow over riblets. First, we confirm for all six riblet shapes that the groove cross-sectional area in viscous units serves as a proxy for the wall-normal permeability that is necessary for the development of Kelvin–Helmholtz rollers. Additionally, we find that the occurrence of the instability correlates with a high momentum absorption at the riblet tips. The momentum absorption can be qualitatively predicted using Stokes flow.
-
ItemDirect Numerical Simulations of Turbulent Flow Over Various Riblet Shapes in Minimal-Span Channels(Springer Verlag, 2020-11-20)Riblets reduce skin-friction drag until their viscous-scaled size becomes large enough for turbulence to approach the wall, leading to the breakdown of drag-reduction. In order to investigate inertial-flow mechanisms that are responsible for the breakdown, we employ the minimal-span channel concept for cost-efficient direct numerical simulation (DNS) of rough-wall flows (MacDonald et al., J. Fluid Mech., vol. 816, 2017, pp. 5–42). This allows us to investigate six different riblet shapes and various viscous-scaled sizes for a total of 21 configurations. We verify that the small numerical domains capture all relevant physics by varying the box size and by comparing to reference data from full-span channel flow. Specifically, we find that, close to the wall in the spectral region occupied by drag-increasing Kelvin–Helmholtz rollers (García-Mayoral & Jiménez, J. Fluid Mech., vol. 678, 2011, pp. 317–347), the energy-difference relative to smooth-wall flow is not affected by the narrow domain, even though these structures have large spanwise extents. This allows us to evaluate the influence of the Kelvin–Helmholtz instability by comparing fluctuations of wall-normal and streamwise velocity, pressure and a passive scalar over riblets of different shapes and viscous-scaled sizes to those over a smooth wall. We observe that triangular riblets with a tip angle a = 30° and blades appear to support the instability, whereas triangular riblets with a = 60°–90° and trapezoidal riblets with a = 30° show little to no evidence of Kelvin–Helmholtz rollers.
-
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.
-
ItemTurbulent flow over spanwise-varying roughness in a minimal streamwise channel(IOP Publishing, 2020-06-10)All rights reserved. We report direct numerical simulations in a minimal streamwise domain of turbulent channel flow over spanwise-alternating patches of rough and smooth walls. Despite the minimal streamwise domain overpredicting streamwise-velocity fluctuations and inhibiting the meandering of long turbulent structures, it captures the rotational behaviour of mean secondary flows also observed in other studies with spanwise-varying roughness. To extend the study of spanwise-varying roughness, we prescribe a lateral velocity to the wall roughness to mimic flow over oblique patches of roughness. Far from the wall, long-lived turbulent structures are convected in the direction of the moving roughness, but their speeds are only weakly perturbed from a preferential value of around 40% of the friction velocity. The turbulence-driven secondary flows laterally convect at comparable speeds, but depend on the roughness patch width.