The minimal-span channel for rough-wall turbulent flows

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.This work was partly funded through the Multi ow program by the European Research Council. Computational time was granted under the Victoria Life Sciences Computational Initiative, which is supported by the Victorian Government, Australia

Turbulent flows over sparse canopies

Turbulent flows over sparse and dense canopies exerting a similar drag force on the flow are investigated using Direct Numerical Simulations. The dense canopies are modelled using a homogeneous drag force, while for the sparse canopy, the geometry of the canopy elements is represented. It is found that on using the friction velocity based on the local shear at each height, the streamwise velocity fluctuations and the Reynolds stress within the sparse canopy are similar to those from a comparable smooth-wall case. In addition, when scaled with the local friction velocity, the intensity of the off-wall peak in the streamwise vorticity for sparse canopies also recovers a value similar to a smooth-wall. This indicates that the sparse canopy does not significantly disturb the near-wall turbulence cycle, but causes its rescaling to an intensity consistent with a lower friction velocity within the canopy. In comparison, the dense canopy is found to have a higher damping effect on the turbulent fluctuations. For the case of the sparse canopy, a peak in the spectral energy density of the wall-normal velocity, and Reynolds stress is observed, which may indicate the formation of Kelvin-Helmholtz-like instabilities. It is also found that a sparse canopy is better modelled by a homogeneous drag applied on the mean flow alone, and not the turbulent fluctuations

Pressure fluctuations and interfacial robustness in turbulent flows over superhydrophobic surfaces

Superhydrophobic surfaces can entrap gas pockets within their grooves when submerged in water. Such a mixed-phase boundary is shown to result in an effective slip velocity on the surface, and has promising potential for drag reduction and energy-saving in hydrodynamic applications. The target flow regime, in most such applications, is a turbulent flow. Previous analyses of this problem involved direct numerical simulations of turbulence with the superhydrophobic surface modelled as a flat boundary, but with a heterogeneous mix of slip and no-slip boundary conditions corresponding to the surface texture. Analysis of the kinematic data from these simulations has helped to establish the magnitude of drag reduction for various texture topologies. The present work is the first investigation that, alongside a kinematic investigation, addresses the robustness of superhydrophobic surfaces by studying the load fields obtain from data from direct numerical simulations (DNS). The key questions at the focus of this work are: does a superhydrophobic surface induce a different pressure field compared to a flat surface? If so, how does this difference scale with system parameters, and when does it become significant that it can deform the air-water interface and potentially rapture the entrapped gas pockets? To this end, we have performed DNS of turbulent channel flows subject to superhydrophobic surfaces over a wide range of texture sizes spanning values from to when expressed in terms of viscous units. The pressure statistics at the wall are decomposed into two contributions: one coherent, caused by the stagnation of slipping flow hitting solid posts, and one time-dependent, caused by overlying turbulence. The results show that the larger texture size intensifies the contribution of stagnation pressure, while the contribution from turbulence is essentially insensitive to . The two-dimensional stagnation pressure distribution at the wall and the pressure statistics in the wall-normal direction are found to be self-similar for different . The scaling of the induced pressure and the consequent deformations of the air-water interface are analysed. Based on our results, an upper bound on the texture wavelength is quantified that limits the range of robust operation of superhydrophobic surfaces when exposed to high-speed flows. Our results indicate that when the system parameters are expressed in terms of viscous units, the main parameters controlling the problem are and a Weber number based on inner dimensions; We obtain good collapse when all our results are expressed in wall units, independently of the Reynolds number

Manipulation of near-wall turbulence by surface slip and permeability

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