Mechanical Engineering - Research Publications
Permanent URI for this collection
Search Results
1 - 10 of 18
-
Item
-
ItemRoughness and Reynolds Number Effects on the Flow Past a Rough-to-Smooth Step Change(Springer International Publishing, 2019)We report direct numerical simulations (DNSs) of open-channel flow with a step change from three-dimensional sinusoidal rough surface to smooth surface. We investigate the persistence of non-equilibrium behaviour beyond this step change (i.e. departures from the equilibrium smooth open-channel flow) and how this depends on (1) roughness virtual origin ϵ/hϵ/h? (scaled by the channel height h), (2) roughness size k / h?, (3) roughness shape? and (4) Reynolds number ReτReτ? To study (1), the roughness origin was placed aligned with, below (step-up) and above (step-down) the smooth patch. To study (2), the equivalent sand-grain roughness of the aligned case was decreased from k+s≃ks+≃ 160 to k+s≃106ks+≃106. To study (3) and (4) the step-down case at Reτ≃395Reτ≃395 was compared with a backward-facing step case at Reτ≃527Reτ≃527, and DNS of square rib rough-to-smooth case at Reτ≃1160Reτ≃1160 (Ismail et al., J. Fluid Mech., vol. 843, 2018, pp. 419–449). Results showed that ϵ/hϵ/h affects the departure from equilibrium by a large extent, while k / h, roughness shape and ReτReτ have a marginal influence. The departure from equilibrium was found to be related to the near-wall amplification of Reynolds shear stress, which in turn depends on ϵ/hϵ/h, i.e. higher ϵ/hϵ/h leads to higher amplification.
-
ItemSecondary motion in turbulent pipe flow with three-dimensional roughness(Cambridge University Press (CUP), 2018-08-31)The occurrence of secondary flows is investigated for three-dimensional sinusoidal roughness where the wavelength and height of the roughness elements are systematically altered. The flow spanned from the transitionally rough regime up to the fully rough regime and the solidity of the roughness ranged from a wavy, sparse roughness to a dense roughness. Analysing the time-averaged velocity, secondary flows are observed in all of the cases, reflected in the coherent stress profile which is dominant in the vicinity of the roughness elements. The roughness sublayer, defined as the region where the coherent stress is non-zero, scales with the roughness wavelength when the roughness is geometrically scaled (proportional increase in both roughness height and wavelength) and when the wavelength increases at fixed roughness height. Premultiplied energy spectra of the streamwise velocity turbulent fluctuations show that energy is reorganised from the largest streamwise wavelengths to the shorter streamwise wavelengths. The peaks in the premultiplied spectra at the streamwise and spanwise wavelengths are correlated with the roughness wavelength in the fully rough regime. Current simulations show that the spanwise scale of roughness determines the occurrence of large-scale secondary flows.
-
ItemDirect numerical simulation of high aspect ratio spanwise-aligned bars(Cambridge University Press (CUP), 2018-03-19)We conduct minimal-channel direct numerical simulations of turbulent flow over two-dimensional rectangular bars aligned in the spanwise direction. This roughness has often been described as d -type, as the roughness function ΔU+ is thought to depend only on the outer-layer length scale (pipe diameter, channel half-height or boundary layer thickness). This is in contrast to conventional engineering rough surfaces, named k -type, for which ΔU+ depends on the roughness height, k. The minimal-span rough-wall channel is used to circumvent the high cost of simulating high Reynolds number flows, enabling a range of bars with varying aspect ratios to be investigated. The present results show that increasing the trough-to-crest height, k, of the roughness while keeping the width between roughness bars, W, fixed in viscous units, results in non- k -type behaviour although this does not necessarily indicate d -type behaviour. Instead, for deep surfaces with k/W≳3, the roughness function appears to depend only on W in viscous units. In these situations, the flow no longer has any information about how deep the roughness is and instead can only ‘see’ the width of the fluid gap between the bars.
-
ItemDirect numerical simulation of high aspect ratio spanwise-aligned bars(Cambridge University Press (CUP), 2017-01-01)We conduct minimal-channel direct numerical simulations of turbulent flow over two-dimensional rectangular bars aligned in the spanwise direction. This roughness has been often described as dtype, as the roughness function ΔU+ is thought to depend only on the outer-layer length scale (pipe diameter, channel half height or boundary layer thickness). This is in contrast to conventional engineering rough surfaces, named k-type, for which ΔU+ depends on the roughness height, k. The minimal-span rough-wall channel is used to circumvent the high cost of simulating high Reynolds number flows, enabling a range of bars with varying aspect ratios to be investigated. The present results show that increasing the trough-to-crest height (k) of the roughness while keeping the width between roughness bars, W, fixed in wall units, results in non-k-type behaviour. The roughness function appears to scale with W, suggesting that this is the only relevant parameter for very deep rough surfaces with k/W≥ 3. In these situations, the flow no longer has any information about how deep the roughness is and instead can only 'see' the width of the fluid gap between the bars.
-
ItemNo Preview AvailableAn experimental investigation into the breakdown of riblet drag reduction at post-optimal conditions(Australasian Fluid Mechanics Society, 2018-01-01)A long-standing question in riblet research is why drag reduction only occurs within a small, non-dimensionally scaled envelope, outside of which drag is significantly increased. For riblets with viscous-scaled spacings that are much larger than those required for drag reduction, one hypothesis is that the riblets exhibit k-type, ‘fully rough’ behaviour. However, this seems counter-intuitive since fully rough behaviour is typically associated with a dominance of pressure drag over viscous drag, and yet riblets can sustain no pressure drag. This study aims to investigate this issue by conducting single normal hot-wire traverses above a trapezoidal riblet surface, over a range of drag-increasing viscous-scaled riblet spacings. Novelty was added by also measuring within the riblet valleys, providing a unique look at the turbulent behaviour within them. Previously proposed mechanisms for the breakdown in drag reduction have included lodgement of turbulence within the riblet valleys, and the development of a Kelvin–Helmholtz instability, but neither mechanism appears active in our results. They instead show a reduction in turbulent energy as riblet spacing increases, despite a significant increase in drag, which does seem to be approaching a k-type roughness asymptote as hypothesised. This may be caused by the generation of time-invariant secondary flows above the riblet tips and corners of the riblet valleys, although this will require further investigation.
-
ItemForm-induced stress in turbulent flow over riblets(Australasian Fluid Mechanics Society, 2018-01-01)We carry out direct numerical simulation of minimal openchannel flow over riblets. Several riblet geometries are simulated, namely symmetric triangular, asymmetric triangular, blade and trapezoidal, and with this unprecedented high-fidelity dataset, we are able to obtain broad insights into the flow physics of riblets. We find that the roughness sublayer thickness, above which the flow is statistically homogeneous, is proportional to the square root of the riblet groove cross-sectional area ℓ+ g in both the drag-reducing and the drag-increasing regime, consistent with the ability of this parameter to collapse the roughness function corresponding to different groove geometries. Large grooves are associated with mean secondary velocities and they carry additional stress that contributes up to 40% of the total shear stress at the crest, comparable to the contribution from the turbulent fluctuations.
-
ItemKelvin–Helmholtz rollers in turbulent flow over riblets(Australasian Fluid Mechanics Society, 2018-01-01)Structures resulting from a Kelvin–Helmholtz instability have been shown to contribute to skin-friction drag in turbulent flow over blade-shaped riblets [4]. Using Direct Numerical Simulation (DNS) data, the present survey of several different riblet shapes reveals that the contribution to wall-shear stress due to the Kelvin–Helmholtz instability depends on riblet shape, in addition to a previously known dependence on riblet size. For a given drag change, sharp triangular and blade riblets promote development of the instability whilst blunt triangular and trapezoidal riblets appear to suppress it.
-
ItemManipulation of near-wall turbulence by surface slip and permeability(IOP PUBLISHING LTD, 2018-01-01)We study the effect on near-wall turbulence of tangential slip and wall-normal transpiration, typically produced by textured surfaces and other surface manipulations. For this, we conduct direct numerical simulations (DNSs) with different virtual origins for the different velocity components. The different origins result in a relative wall-normal displacement of the near-wall, quasi-streamwise vortices with respect to the mean flow, which in turn produces a change in drag. The objective of this work is to extend the existing understanding on how these virtual origins affect the flow. In the literature, the virtual origins for the tangential velocities are typically characterised by slip boundary conditions, while the wall-normal velocity is assumed to be zero at the boundary plane. Here we explore different techniques to define and implement the three virtual origins, with special emphasis on the wall-normal one. We investigate impedance conditions relating the wall-normal velocity to the pressure, and linear relations between the velocity components and their wall-normal gradients, as is typically done to impose slip conditions. These models are first tested to represent a smooth wall below the boundary plane, with all virtual origins equal, and later for different tangential and wall-normal origins. Our results confirm that the change in drag is determined by the offset between the origins perceived by mean flow and the quasi-streamwise vortices or, more generally, the near-wall turbulent cycle. The origin for the latter, however, is not set by the spanwise virtual origin alone, as previously proposed, but by a combination of the spanwise and wall-normal origins, and mainly determined by the shallowest of the two. These observations allow us to extend the existing expression to predict the change in drag, accounting for the wall-normal effect when the transpiration is not negligible.
-
ItemThe minimal-span channel for rough-wall turbulent flows(Cambridge University Press (CUP), 2017-04-10)Roughness predominantly alters the near-wall region of turbulent flow while the outer layer remains similar with respect to the wall shear stress. This makes it a prime candidate for the minimal-span channel, which only captures the near-wall flow by restricting the spanwise channel width to be of the order of a few hundred viscous units. Recently, Chung et al. (J. Fluid Mech., vol. 773, 2015, pp. 418–431) showed that a minimal-span channel can accurately characterise the hydraulic behaviour of roughness. Following this, we aim to investigate the fundamental dynamics of the minimal-span channel framework with an eye towards further improving performance. The streamwise domain length of the channel is investigated with the minimum length found to be three times the spanwise width or 1000 viscous units, whichever is longer. The outer layer of the minimal channel is inherently unphysical and as such alterations to it can be performed so long as the near-wall flow, which is the same as in a full-span channel, remains unchanged. Firstly, a half-height (open) channel with slip wall is shown to reproduce the near-wall behaviour seen in a standard channel, but with half the number of grid points. Next, a forcing model is introduced into the outer layer of a half-height channel. This reduces the high streamwise velocity associated with the minimal channel and allows for a larger computational time step. Finally, an investigation is conducted to see if varying the roughness Reynolds number with time is a feasible method for obtaining the full hydraulic behaviour of a rough surface. Currently, multiple steady simulations at fixed roughness Reynolds numbers are needed to obtain this behaviour. The results indicate that the non-dimensional pressure gradient parameter must be kept below 0.03–0.07 to ensure that pressure gradient effects do not lead to an inaccurate roughness function. An empirical costing argument is developed to determine the cost in terms of CPU hours of minimal-span channel simulations a priori. This argument involves counting the number of eddy lifespans in the channel, which is then related to the statistical uncertainty of the streamwise velocity. For a given statistical uncertainty in the roughness function, this can then be used to determine the simulation run time. Following this, a finite-volume code with a body-fitted grid is used to determine the roughness function for square-based pyramids using the above insights. Comparisons to experimental studies for the same roughness geometry are made and good agreement is observed.