ORCID Profile
0000-0001-6348-6392
Current Organisation
The University of Auckland
Does something not look right? The information on this page has been harvested from data sources that may not be up to date. We continue to work with information providers to improve coverage and quality. To report an issue, use the Feedback Form.
Publisher: IOP Publishing
Date: 04-2016
Publisher: The University of Queensland
Date: 11-12-2020
DOI: 10.14264/5586462
Publisher: American Meteorological Society
Date: 09-2020
Abstract: We present a turbulent kinetic energy (TKE) closure scheme for the stably stratified atmosphere in which the mixing lengths for momentum and heat are not parameterized in the same manner. The key difference is that, while the mixing length for heat tends toward the stability independent mixing length for momentum in neutrally stratified conditions, it tends toward one based on the Brunt–Väisälä time scale and square root of the TKE in the limit of large stability. This enables a unique steady-state solution for TKE to be obtained, which we demonstrate would otherwise be impossible if the mixing lengths were the same. Despite the model’s relative simplicity, it is shown to perform reasonably well with observational data from the 1999 Cooperative Atmosphere–Surface Exchange Study (CASES-99) using commonly employed model constants. Analyzing the scaling behavior of the nondimensional velocity and potential temperature gradients, or of the stability (correction) functions, reveals that for large stability the present model scales in the same manner as the first-order operational scheme of Viterbo et al. Alternatively, it appears as a blend of two cases of the TKE closure scheme of Baas et al. Critically, because a unique steady-state TKE can be obtained, the present model avoids the nonphysical behavior identified in one of the cases of Baas et al.
Publisher: Cambridge University Press (CUP)
Date: 27-04-2015
DOI: 10.1017/JFM.2015.172
Abstract: Direct numerical simulations (DNS) are conducted for turbulent flow through pipes with three-dimensional sinusoidal roughnesses explicitly represented by body-conforming grids. The same viscous-scaled roughness geometry is first simulated at a range of different Reynolds numbers to investigate the effects of low Reynolds numbers and low $R_{0}/h$ , where $R_{0}$ is the pipe radius and $h$ is the roughness height. Results for the present class of surfaces show that the Hama roughness function ${\\rm\\Delta}U^{+}$ is only marginally affected by low Reynolds numbers (or low $R_{0}/h$ ), and observations of outer-layer similarity (or lack thereof) show no signs of sensitivity to Reynolds number. Then, building on this, a systematic approach is taken to isolate the effects of roughness height $h^{+}$ and wavelength ${\\it\\lambda}^{+}$ in a turbulent wall-bounded flow in both transitionally rough and fully rough regimes. Current findings show that while the effective slope $\\mathit{ES}$ (which for the present sinusoidal surfaces is proportional to $h^{+}/{\\it\\lambda}^{+}$ ) is an important roughness parameter, the roughness function ${\\rm\\Delta}U^{+}$ must also depend on some measure of the viscous roughness height. A simplistic linear–log fit clearly illustrates the strong correlation between ${\\rm\\Delta}U^{+}$ and both the roughness average height $k_{a}^{+}$ (which is related to $h^{+}$ ) and $\\mathit{ES}$ for the surfaces simulated here, consistent with published literature. Various definitions of the virtual origin for rough-wall turbulent pipe flow are investigated and, for the surfaces simulated here, the hydraulic radius of the pipe appears to be the most suitable parameter, and indeed is the only virtual origin that can ever lead to collapse in the total stress. First- and second-order statistics are also analysed and collapses in the outer layer are observed for all cases, including those where the largest roughness height is a substantial proportion of the reference radius (low $R_{0}/h$ ). These results provide evidence that turbulent pipe flow over the present sinusoidal surfaces adheres to Townsend’s notion of outer-layer similarity, which pertains to statistics of relative motion.
Publisher: Cambridge University Press (CUP)
Date: 26-05-2015
DOI: 10.1017/JFM.2015.230
Abstract: We describe a fast direct numerical simulation (DNS) method that promises to directly characterise the hydraulic roughness of any given rough surface, from the hydraulically smooth to the fully rough regime. The method circumvents the unfavourable computational cost associated with simulating high-Reynolds-number flows by employing minimal-span channels (Jiménez & Moin, J. Fluid Mech. , vol. 225, 1991, pp. 213–240). Proof-of-concept simulations demonstrate that flows in minimal-span channels are sufficient for capturing the downward velocity shift, that is, the Hama roughness function, predicted by flows in full-span channels. We consider two sets of simulations, first with modelled roughness imposed by body forces, and second with explicit roughness described by roughness-conforming grids. Owing to the minimal cost, we are able to conduct direct numerical simulations with increasing roughness Reynolds numbers while maintaining a fixed blockage ratio, as is typical in full-scale applications. The present method promises a practical, fast and accurate tool for characterising hydraulic resistance directly from profilometry data of rough surfaces.
Publisher: IOP Publishing
Date: 11-04-2017
Publisher: American Physical Society (APS)
Date: 22-07-2019
Publisher: Cambridge University Press (CUP)
Date: 08-09-2016
DOI: 10.1017/JFM.2016.459
Abstract: We investigate rough-wall turbulent flows through direct numerical simulations of flow over three-dimensional transitionally rough sinusoidal surfaces. The roughness Reynolds number is fixed at $k^{+}=10$ , where $k$ is the sinusoidal semi- litude, and the sinusoidal wavelength is varied, resulting in the roughness solidity $\\unicode[STIX]{x1D6EC}$ (frontal area ided by plan area) ranging from 0.05 to 0.54. The high cost of resolving both the flow around the dense roughness elements and the bulk flow is circumvented by the use of the minimal-span channel technique, recently demonstrated by Chung et al. ( J. Fluid Mech. , vol. 773, 2015, pp. 418–431) to accurately determine the Hama roughness function, $\\unicode[STIX]{x0394}U^{+}$ . Good agreement of the second-order statistics in the near-wall roughness-affected region between minimal- and full-span rough-wall channels is observed. In the sparse regime of roughness ( $\\unicode[STIX]{x1D6EC}\\lesssim 0.15$ ) the roughness function increases with increasing solidity, while in the dense regime the roughness function decreases with increasing solidity. It was found that the dense regime begins when $\\unicode[STIX]{x1D6EC}\\gtrsim 0.15{-}0.18$ , in agreement with the literature. A model is proposed for the limit of $\\unicode[STIX]{x1D6EC}\\rightarrow \\infty$ , which is a smooth wall located at the crest of the roughness elements. This model assists with interpreting the asymptotic behaviour of the roughness, and the rough-wall data presented in this paper show that the near-wall flow is tending towards this modelled limit. The peak streamwise turbulence intensity, which is associated with the turbulent near-wall cycle, is seen to move further away from the wall with increasing solidity. In the sparse regime, increasing $\\unicode[STIX]{x1D6EC}$ reduces the streamwise turbulent energy associated with the near-wall cycle, while increasing $\\unicode[STIX]{x1D6EC}$ in the dense regime increases turbulent energy. An analysis of the difference of the integrated mean momentum balance between smooth- and rough-wall flows reveals that the roughness function decreases in the dense regime due to a reduction in the Reynolds shear stress. This is predominantly due to the near-wall cycle being pushed away from the roughness elements, which leads to a reduction in turbulent energy in the region previously occupied by these events.
Publisher: Cambridge University Press (CUP)
Date: 19-12-2019
DOI: 10.1017/JFM.2018.900
Abstract: We conducted direct numerical simulations of turbulent flow over three-dimensional sinusoidal roughness in a channel. A passive scalar is present in the flow with Prandtl number $Pr=0.7$ , to study heat transfer by forced convection over this rough surface. The minimal-span channel is used to circumvent the high cost of simulating high-Reynolds-number flows, which enables a range of rough surfaces to be efficiently simulated. The near-wall temperature profile in the minimal-span channel agrees well with that of the conventional full-span channel, indicating that it can be readily used for heat-transfer studies at a much reduced cost compared to conventional direct numerical simulation. As the roughness Reynolds number, $k^{+}$ , is increased, the Hama roughness function, $\\unicode[STIX]{x0394}U^{+}$ , increases in the transitionally rough regime before tending towards the fully rough asymptote of $\\unicode[STIX]{x1D705}_{m}^{-1}\\log (k^{+})+C$ , where $C$ is a constant that depends on the particular roughness geometry and $\\unicode[STIX]{x1D705}_{m}\\approx 0.4$ is the von Kármán constant. In this fully rough regime, the skin-friction coefficient is constant with bulk Reynolds number, $Re_{b}$ . Meanwhile, the temperature difference between smooth- and rough-wall flows, $\\unicode[STIX]{x0394}\\unicode[STIX]{x1D6E9}^{+}$ , appears to tend towards a constant value, $\\unicode[STIX]{x0394}\\unicode[STIX]{x1D6E9}_{FR}^{+}$ . This corresponds to the Stanton number (the temperature analogue of the skin-friction coefficient) monotonically decreasing with $Re_{b}$ in the fully rough regime. Using shifted logarithmic velocity and temperature profiles, the heat-transfer law as described by the Stanton number in the fully rough regime can be derived once both the equivalent sand-grain roughness $k_{s}/k$ and the temperature difference $\\unicode[STIX]{x0394}\\unicode[STIX]{x1D6E9}_{FR}^{+}$ are known. In meteorology, this corresponds to the ratio of momentum and heat-transfer roughness lengths, $z_{0m}/z_{0h}$ , being linearly proportional to the inner-normalised momentum roughness length, $z_{0m}^{+}$ , where the constant of proportionality is related to $\\unicode[STIX]{x0394}\\unicode[STIX]{x1D6E9}_{FR}^{+}$ . While Reynolds analogy, or similarity between momentum and heat transfer, breaks down for the bulk skin-friction and heat-transfer coefficients, similar distribution patterns between the heat flux and viscous component of the wall shear stress are observed. Instantaneous visualisations of the temperature field show a thin thermal diffusive sublayer following the roughness geometry in the fully rough regime, resembling the viscous sublayer of a contorted smooth wall.
Publisher: IOP Publishing
Date: 09-2020
DOI: 10.1088/1742-6596/1618/3/032008
Abstract: A scaling or nondimensionalisation of atmospheric buoyancy vortices for power generation (based on the Oberbeck-Boussinesq assumption) is proposed that uses a published formulation from the study of Rotating Rayleigh-Benard Convection. This is combined with assumptions that the vortex flows are pseudo-cyclostrophic and that a radial Richardson number can serve as a predictor of the onset of Kelvin-Helmholtz instability leading to a transition to a turbulent plume, in order to locate the cold reservoir of the vortex when viewed as a heat engine. This permits the prediction of the behaviour of large vortices in atmosphere using data from experiments on small vortices.
Publisher: The Royal Society
Date: 08-07-2014
Abstract: In this study, we consider suspensions of swimming microorganisms in situations where we might expect the promotion of two-dimensional flow, such as within thin fluid films. Given that two-dimensional, inertialess flows are notoriously long-ranged (although not afflicted by Stokes paradox in the case of self-motile bodies), this raises interesting questions around the care which must be taken with the semi-dilute assumption in such situations. Adopting the prototype squirmer as a model of a swimming microorganisms of the type previously considered, we find that although the flowfield decays algebraically with the characteristic separation distance between microorganisms, there remains a finite interaction between the squirmers even at asymptotically large distances. This finding is further borne out by asymptotic analysis, which confirms that the limiting form of the far-field interaction depends solely upon the relative orientation between the microorganisms. Those which swim in the same general direction are seen to experience very large lateral displacements (many times the size of the displacements experienced owing to interactions between less well-aligned swimmers). This clearly has potential implications for very dilute suspensions in which squirmers become broadly aligned in their swimming direction (e.g. during chemotaxis). We show that hydrodynamically enhanced cell spreading, previously reported for denser suspensions, can persist even at extreme dilutions. Moreover, we demonstrate that this induced spreading can continue in the presence of potentially decohering Brownian effects.
Publisher: Springer Science and Business Media LLC
Date: 20-11-2020
Publisher: The University of Queensland
Date: 11-12-2020
DOI: 10.14264/278BAA7
Publisher: Springer Science and Business Media LLC
Date: 27-03-2020
Publisher: Cambridge University Press (CUP)
Date: 28-02-2017
DOI: 10.1017/JFM.2017.69
Abstract: 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.
Publisher: Cambridge University Press (CUP)
Date: 31-08-2018
DOI: 10.1017/JFM.2018.570
Abstract: 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.
Publisher: Elsevier BV
Date: 11-2022
Publisher: Cambridge University Press (CUP)
Date: 19-03-2018
DOI: 10.1017/JFM.2018.150
Abstract: 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 $\\unicode[STIX]{x0394}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 $\\unicode[STIX]{x0394}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, ${\\mathcal{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/{\\mathcal{W}}\\gtrsim 3$ , the roughness function appears to depend only on ${\\mathcal{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.
Publisher: IOP Publishing
Date: 04-2018
Location: United States of America
No related grants have been discovered for Michael MacDonald.