ORCID Profile
0000-0002-7157-6440
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.
In Research Link Australia (RLA), "Research Topics" refer to ANZSRC FOR and SEO codes. These topics are either sourced from ANZSRC FOR and SEO codes listed in researchers' related grants or generated by a large language model (LLM) based on their publications.
Interdisciplinary Engineering | Turbulent Flows | Fluidization And Fluid Mechanics | Turbulent Flows | Computational Fluid Dynamics | Fluidisation and Fluid Mechanics | Heat And Mass Transfer Operations | Environmental Engineering | Environmental Engineering Modelling | Computational Heat Transfer | Water And Sanitary Engineering |
Physical sciences | Water Allocation and Quantification | Industrial Energy Conservation and Efficiency | Forest and Woodlands Water Management | Atmospheric Processes and Dynamics | Industry | Land and water management | Residential Energy Conservation and Efficiency | Land and water management | Land and water management | Commercial Energy Conservation and Efficiency | Expanding Knowledge in the Earth Sciences | Coastal and Estuarine Water Management | Farmland, Arable Cropland and Permanent Cropland Water Management
Publisher: Wiley
Date: 15-10-2021
Publisher: AIP Publishing
Date: 03-2022
DOI: 10.1063/5.0083230
Abstract: Direct numerical simulation (DNS) results for turbulent open-channel flow through an idealized sine-generated meander with and without an internal heat source that models radiative heating from above are used to analyze the effect of a very sharp meander configuration and thermal stratification on the turbulence structure in the channel with friction Reynolds number Reτ=200. Spatial distributions of temperature, mean velocities, vorticity, mean-flow kinetic energy, and turbulent kinetic energy (TKE) are presented. In both cases, the cross-sectional motion is characterized by three circulation cells: a center-region cell and two weaker outer bank and inner bank cells. However, there is also a small cell observed near the corner of the channel bed inner bank at the channel outlet and the channel bed outer bank at the channel inlet. The tri-cellular cross-stream motions control the distributions of temperature and kinetic energy. In the stratified case, two separated shear layers (SSLs) are found: the first one is formed before the bend apex, and the second one is observed in the wake after the bend apex. In the neutral case, only the first SSL is observed. Turbulent lification can be seen in both cases however, in the stratified case, the second SSL stretches out to the channel outlet and is introduced back to the channel inlet by an anti-symmetric periodic boundary condition and then follows the outer bank line. The two SSLs converge in the region before the bend apex and lify the turbulence more strongly there than in the neutral case. The turbulence kinetic energy budget terms for the stratified case are analyzed to determine the characteristics of production, dissipation and transport of TKE in thermally stratified meandering flow.
Publisher: Trans Tech Publications, Ltd.
Date: 05-2014
DOI: 10.4028/WWW.SCIENTIFIC.NET/AMM.553.193
Abstract: Modern GPUs (graphical processing units) are a common source of processing power inmany supercomputers. Their performance derives from the highly parallel architecture that is em-ployed and have the benefit of low cost, temperature and power consumption. Two finite differencemodels have been implemented on GPU, a semi-implicit and an explicit algorithm, to numericallymodel a stratified shear layer, that needs fine meshes to be modelled accurately. The GPU modelswere shown to improve performance by factors of around 50x and 20x for the semi-implicit and ex-plicit models respectively.
Publisher: Elsevier BV
Date: 07-2013
Publisher: Elsevier BV
Date: 04-2020
Publisher: Cambridge University Press (CUP)
Date: 08-2019
DOI: 10.1017/JFM.2019.543
Abstract: Evolution of thermally stratified open channel flow after removal of a volumetric heat source is investigated using direct numerical simulation. The heat source models radiative heating from above and varies with height due to progressive absorption. After removal of the heat source the initial stable stratification breaks down and the channel approaches a fully mixed isothermal state. The initial state consists of three distinct regions: a near-wall region where stratification plays only a minor role, a central region where stratification has a significant effect on flow dynamics and a near-surface region where buoyancy effects dominate. We find that a state of local energetic equilibrium observed in the central region of the channel in the initial state persists until the late stages of the destratification process. In this region local turbulence parameters such as eddy diffusivity $k_{h}$ and flux Richardson number $R_{f}$ are found to be functions only of the Prandtl number $Pr$ and a mixed parameter ${\\mathcal{Q}}$ , which is equal to the ratio of the local buoyancy Reynolds number $Re_{b}$ and the friction Reynolds number $Re_{\\unicode[STIX]{x1D70F}}$ . Close to the top and bottom boundaries turbulence is also affected by $Re_{\\unicode[STIX]{x1D70F}}$ and vertical position $z$ . In the initial heated equilibrium state the laminar surface layer is stabilised by the heat source, which acts as a potential energy sink. Removal of the heat source allows Kelvin–Helmholtz-like shear instabilities to form that lead to a rapid transition to turbulence and significantly enhance the mixing process. The destratifying flow is found to be governed by bulk parameters $Re_{\\unicode[STIX]{x1D70F}}$ , $Pr$ and the friction Richardson number $Ri_{\\unicode[STIX]{x1D70F}}$ . The overall destratification rate ${\\mathcal{D}}$ is found to be a function of $Ri_{\\unicode[STIX]{x1D70F}}$ and $Pr$ .
Publisher: Begell House
Date: 2011
Publisher: Cambridge University Press (CUP)
Date: 29-12-2021
DOI: 10.1017/JFM.2020.921
Publisher: Cambridge University Press (CUP)
Date: 25-01-2021
Publisher: Cambridge University Press (CUP)
Date: 28-07-2020
DOI: 10.1017/JFM.2020.447
Publisher: Cambridge University Press (CUP)
Date: 30-08-2013
DOI: 10.1017/JFM.2013.392
Abstract: This study considers the convective-type instability of the near-field flow of a planar, pure thermal plume with a finite area source. Previous studies revealed the existence of an off-axis thermal boundary-layer instability, driving a puffing instability in the central ascending column, and qualitatively showed correlations between instabilities in these two flow regions. This paper extends the analysis to examine the effect of Prandtl number on transitional near-field behaviours and reports on the stability characteristics of a near-field, pure thermal plume based on a direct stability analysis. The variations in flow behaviours in response to symmetric and asymmetric disturbances suggest the existence of coupled instability mechanisms in the off-axis thermal boundary layer and the central ascending column.
Publisher: Springer Science and Business Media LLC
Date: 2003
Publisher: AIP Publishing
Date: 02-2013
DOI: 10.1063/1.4792709
Abstract: Numerical simulations were conducted of turbulent natural convection in a shallow tetrahedron domain representing the sidearm of a lake or water reservoir. The tetrahedron cavity is a more realistic, three-dimensional approximation of a lake or reservoir sidearm than the two-dimensional triangle cavity often seen in the literature. Lateral temperature gradients exist due to the varying depth of the cavity, resulting in lateral circulation. These flows are important in a reservoir as they can carry with them particles and various pollutants, transporting and mixing them with the central section. Therefore, study in this area is important in water quality management. The simulations use a Cartesian grid with an Immersed Boundary Method for the sloped bottom surfaces. Heat input is through a solar radiation model consisting of a heat flux from the sloped bottom boundaries and an internal heating source term in the body of the water. Also studied is the night time model where cooling is through a heat flux at the top boundary. Scaling analysis from the literature is extended to suit the new geometry and numerical simulations are used to validate the results. The numerical simulations include calculating horizontal heat transfer profiles, volumetric flow rates, and analysis of complex flow features. The extension to three dimensions results in significant changes to the flow and introduces some complex features, such as helical flow both towards and away from the tip. Some general parameterisations are proposed for the tetrahedron cavity based on the numerical simulations.
Publisher: Wiley
Date: 19-03-2013
DOI: 10.1002/FLD.3789
Publisher: Cambridge University Press (CUP)
Date: 29-11-2018
DOI: 10.1017/JFM.2017.796
Abstract: The entrainment of fluid across a sheared density interface has been examined experimentally in a purging cavity flow. In this flow, a long straight cavity with sloped entry and exit boundaries is located in the base of a straight open channel. Dense cavity fluid is entrained from the cavity into the turbulent overflow. The cavity geometry has been designed to ensure there is no separation of the overflow in the cavity region, with the goal of avoiding cavity-specific entrainment mechanisms as have been encountered in most previous experiments using similar arrangements. Results are obtained over a bulk Richardson number range $Ri_{b}=g\\unicode[STIX]{x0394}\\unicode[STIX]{x1D70C}D/\\unicode[STIX]{x1D70C}_{0}U_{b}^{2}=1$ to 19, where $D$ and $U_{b}$ are the depth of the mixed layer and bulk velocity in the cavity, respectively. The experiments cover the Reynolds number range $Re=U_{b}D/\\unicode[STIX]{x1D708}=7100$ to 15 100 and interface length to mixed layer depth ratios of 2.4 to 16. Particle image velocimetry and laser induced fluorescence measurements indicate the flow regime over this entire range is one dominated by the Holmboe wave instability. The non-dimensional entrainment rate, $E=u_{e}/U_{b}$ , is shown to scale with the bulk Richardson number. We find that the entrainment scaling $E=CRi_{b}^{-1.38}$ applies over the entire experimental range, with no apparent dependence on interface length. The exponent in the scaling is similar to previous non-cavity-based sheared interface flows, however, the constant $C$ is up to an order of magnitude smaller. Close agreement is, however, obtained by instead correlating entrainment with the local gradient Richardson number centred on the interface, rather than bulk quantities. We obtain $E=0.0021Ri_{g}^{-0.63}$ for data over $10 Ri_{g} $ , where $Ri_{g}=\\langle g\\unicode[STIX]{x2202}\\unicode[STIX]{x1D70C}/\\unicode[STIX]{x1D70C}_{0}\\unicode[STIX]{x2202}z\\rangle /\\langle (\\unicode[STIX]{x2202}U/\\unicode[STIX]{x2202}z)^{2}\\rangle$ . The density interface is much thinner and therefore more stable in the present flow configuration compared with other published results for the same bulk Richardson number. We suggest that our configuration ensures a sharp mixing layer profile at the upstream end of the cavity even at relatively low bulk Richardson numbers of $Ri_{b}=1$ and that the reduced mixing in the Holmboe wave regime allows the interface to retain its sharp character over the cavity length, resulting in weak sensitivity to cavity length.
Publisher: Elsevier BV
Date: 2003
Publisher: American Geophysical Union (AGU)
Date: 12-2012
DOI: 10.1029/2012WR012218
Publisher: Cambridge University Press (CUP)
Date: 31-03-2022
DOI: 10.1017/JFM.2022.152
Abstract: Turbulent fountain flow consists of two distinct stages, the initial ‘negatively buoyant jet’ (NBJ) stage, and the fully developed ‘fountain’ stage. The present study investigates both stages of the flow using particle image velocimetry and planar laser-induced fluorescence, over a range of source Froude numbers, $10\\lesssim Fr_o\\lesssim 30$ , and Reynolds numbers, $5500\\lesssim Re_o\\lesssim 7700$ . While the velocity and buoyancy profiles in NBJs take similar Gaussian shapes over a wide range of axial locations, this was not observed in fountains. The changing profile shape is most evident in the outer flow (OF) region, while there is a degree of similarity in the inner flow (IF). Entrainment between IF and OF is shown to depend on the local Richardson number, $Ri$ . The fountains are found to have a negative entrainment coefficient, $\\alpha $ , for the majority of their height, implying a net radial outflow of fluid from the IF to the OF. An alternative description of entrainment is considered, the ‘decomposed top-hat’ model, where the radial flow is separated into inflow and outflow components that are then estimated using the present experimental data. The inflow component was found to be proportional to the axial IF velocity, which is similar to the classical description of entrainment in pure jets lumes, while the outflow depends on the local $Ri$ . Entrainment in NBJs may also be described by this framework, which, despite not having an OF, is still subject to an $Ri$ -dependent radial outflow.
Publisher: AIP Publishing
Date: 02-2012
DOI: 10.1063/1.3680861
Abstract: Experimental data and large eddy simulation results are analysed to investigate shear driven entrainment of a negatively buoyant fluid from trapezoidal depressions and cavities. This flow is of relevance to a number of environmentally significant applications including purging of saline pools in rivers and pollutant dispersion in cities and towns situated within topographic depressions. New scaling relations for the entrainment rate are developed based on physical arguments. Our scaling relations are shown to agree well with both experiments and numerical simulations of this flow in trapezoidal cavities with aspect ratios ranging between 7 and 17, entry beach angles between 8° and 33°, and an exit beach angle of 33°. For the numerical simulations, a sub-filter scale turbulence model is used that combines the dynamic mixed model of Zang et al. [“A dynamic mixed subgrid-scale model and its application to recirculating flows,” Phys. Fluids A 5, 3186 (1993)]10.1063/1.858675 with the dynamic localization procedure of Piomelli and Liu [“Large eddy simulation of rotating channel flows using a localized dynamic model,” Phys. Fluids 7, 839 (1995)]10.1063/1.868607 and a buoyancy correction similar to that proposed by Brown et al. [“Large-eddy simulation of stable atmospheric boundary-layers with a revised stochastic subgrid model,” Q. J. R. Meteorol. Soc. 120, 1485 (1994)]10.1002/qj.49712052004 based on data measured in the stable atmospheric boundary layer. Simulations run using this model give significantly closer agreement with the experimental data than simulations run using a version of the dynamic Smagorinsky model with similar modifications. Support for the physical arguments upon which the scaling relations are based is obtained through a statistical analysis of the turbulent flow fields generated by the numerical simulations.
Publisher: Elsevier BV
Date: 04-2020
Publisher: Elsevier BV
Date: 03-2017
Publisher: Elsevier BV
Date: 04-2020
Publisher: Elsevier BV
Date: 02-2017
Publisher: Informa UK Limited
Date: 12-09-2019
Publisher: Elsevier BV
Date: 06-2016
Publisher: Cambridge University Press (CUP)
Date: 09-02-2015
DOI: 10.1017/JFM.2014.711
Abstract: Direct numerical simulations (DNS) of turbulent stratified flow in an open channel with an internal heat source following the Beer–Lambert law from the surface are used to investigate the transition from neutral to strongly stable flow. Our buoyancy bulk parameter is defined through the ratio of the domain height ${\\it\\delta}$ to $\\mathscr{L}$ , a bulk Obukhov length scale for the flow. We cover the range ${\\it\\lambda}={\\it\\delta}/\\mathscr{L}=0{-}2.0$ , from neutral conditions to the onset of the stable regime, with the Reynolds number range $Re_{{\\it\\tau}}=200{-}800$ , at a Prandtl number of 0.71. The result is a boundary layer flow where the effects of stratification are weak in the wall region but progressively stronger in the outer layer up to the free surface. At ${\\it\\lambda}\\simeq 1$ the turbulent kinetic energy (TKE) budget is in local equilibrium over a region extending from the near-wall region to a free-surface affected region a distance $l_{{\\it\\nu}}$ from the surface, with $l_{{\\it\\nu}}/{\\it\\delta}\\sim Re^{-1/2}$ . In this equilibrium region the flow can be characterised by the flux Richardson number $R_{f}$ and the local Obukhov length scale ${\\it\\Lambda}$ . At higher ${\\it\\lambda}$ local mixing limit conditions are observed over an extended region. At ${\\it\\lambda}=2$ the flux Richardson number approaches critical limit values of $R_{f,c}\\simeq 0.18$ and gradient Richardson number $Ri_{c}\\simeq 0.2$ . At high ${\\it\\lambda}$ , we obtain a flow field where buoyancy interacts with the smallest scales of motion and the turbulent shear stress and buoyancy flux are suppressed to molecular levels. We find that this regime can be identified in terms of the parameter $Re_{\\mathscr{L},c}=\\mathscr{L}u_{{\\it\\tau}}/{\\it\\nu}\\lesssim 200{-}400$ (where $u_{{\\it\\tau}}$ is the friction velocity and ${\\it\\nu}$ the kinematic viscosity), which is related to the $L_{\\ast }$ parameter of Flores and Riley ( Boundary-Layer Meteorol. , vol. 139 (2), 2011, pp. 241–259) and buoyancy Reynolds number $\\mathscr{R}$ . With energetic equilibrium attained, the local buoyancy Reynolds number, $Re_{{\\it\\Lambda}}={\\it\\Lambda}\\langle u^{\\prime }w^{\\prime }\\rangle ^{1/2}/{\\it\\nu}$ , is directly related to the separation of the Ozmidov ( $l_{O}$ ) and Kolmogorov ( ${\\it\\eta}$ ) length scales in the outer boundary layer by $Re_{{\\it\\Lambda}}\\simeq \\mathscr{R}\\equiv (l_{O}/{\\it\\eta})^{4/3}$ . The inner wall region has the behaviour $\\mathscr{R}\\sim Re_{\\mathscr{L}}Re_{{\\it\\tau}}$ , in contrast to stratified boundary layer flows where the buoyancy flux is non-zero at the wall and $\\mathscr{R}\\sim Re_{\\mathscr{L}}$ .
Publisher: Elsevier BV
Date: 2005
Publisher: IOP Publishing
Date: 06-2010
Publisher: AIP Publishing
Date: 09-2022
DOI: 10.1063/5.0106526
Abstract: Thermal stratification can lead to the d ing of turbulence, which reduces the mixing of solutes in a fluid body. A series of direct numerical simulation (DNS) solutions sweeping through a range of four different meandering channel curvatures, from a sharp to mild curvature range, are obtained to investigate the effect of curvature on stratification in meandering thermally stratified turbulent open channel flow with an internal heat source that models radiative heating from above. Based on the DNS results, the present paper addresses two issues. First, the influence of changing curvature on the complex bi-cellular pattern of the secondary flow is investigated, including the distribution of the temperature field. Second, the effects of changing curvature on the degree of stratification are analyzed. Stratification can be characterized by the friction Richardson number Riτ and the bulk Richardson number Rib. Stratification can also be viewed in terms of the transfer of energy from mean flow kinetic energy to potential energy via buoyancy fluxes. We study the effect of curvature on stratification by investigating its effect on the friction and bulk Richardson numbers. We also study the transfers between the global potential and kinetic energy reservoirs, including the global available Ea, background Eb, and total potential energy Ep, and the domain-averaged mean kinetic and turbulent kinetic energy. It is found that, in meandering channels, with the increase in curvature, Ep increases and Riτ and Rib decrease, indicating that increasing curvature leads to a decrease in the level of stratification. On the other hand, we also find that a low curvature meandering channel has a higher level of stratification than a straight channel.
Publisher: Elsevier BV
Date: 09-2013
DOI: 10.1016/J.WATRES.2013.05.045
Abstract: Large Eddy Simulation (LES) is used to model two lake sidearms subject to heating from solar radiation and cooling from a surface flux. The sidearms are part of Lake Audrey, NJ, USA and Lake Alexandrina, SA, Australia. The simulation domains are created using bathymetry data and the boundary is modelled with an Immersed Boundary Method. We investigate the cooling and heating phases with separate quasi-steady state simulations. Differential heating occurs in the cavity due to the changing depth. The resulting temperature gradients drive lateral flows. These flows are the dominant transport process in the absence of wind. Study in this area is important in water quality management as the lateral circulation can carry particles and various pollutants, transporting them to and mixing them with the main lake body.
Publisher: Trans Tech Publications, Ltd.
Date: 2014
DOI: 10.4028/WWW.SCIENTIFIC.NET/AMM.493.215
Abstract: Parallel performance of a fractional step Navier-Stokes solver is investigated. Parallelisation is performed using Message Passing Interface, with domain partitioning. Block preconditioning is applied to the solution of the pressure Poisson equation, which is often the bottleneck in the computation of the fractional step method. Preconditioners tested are classes of incomplete matrix decompositions and sparse approximate inverses. The computational domain is decomposed into eight parts of about equal size in terms of the number of cells, and solved on eight parallel processors. Several aspects of the parallelisation, such as domain splitting directions, speed-up and scalability of the preconditioners, are discussed.
Publisher: Cambridge University Press (CUP)
Date: 28-11-2012
DOI: 10.1017/JFM.2011.471
Abstract: Numerical evidence is presented for previously unreported flow behaviour in a two-dimensional rectangular side-heated cavity partitioned in the centre by vertical wall with an infinite conductivity. In this flow heat is transferred between both sides of the cavity through the conducting wall with natural convection boundary layers forming on all vertical surfaces. Simulations have been conducted over the range of Rayleigh numbers $\\mathit{Ra}= 0. 6\\text{{\\ndash}} 1. 6\\ensuremath{\\times} 1{0}^{10} $ at Prandtl number $\\mathit{Pr}= 7. 5$ and at aspect ratios of $H/ W= 1\\text{{\\ndash}} 2$ where $H$ and $W$ are the height and width of the cavity. It was found that the thermal coupling of the boundary layers on either side of the conducting partition causes the cavity flow to become absolutely unstable for a Rayleigh number at which otherwise similar non-partitioned cavity flow is steady but convectively unstable. Additionally, unlike the non-partitioned cavity, which eventually bifurcates to a multi-modal oscillatory regime, this bifurcation is manifested as a single mode oscillation with ${f}^{+ } = f{\\nu }^{1/ 3} / \\mathop{ (g\\ensuremath{\\beta} \\mrm{\\Delta} \\theta )}\\nolimits ^{2/ 3} \\approx 0. 0145$ , where $ \\mrm{\\Delta} \\theta $ is the temperature difference between the hot and cold walls, $g$ is the gravitational acceleration, $f$ is the oscillation frequency and $\\nu $ and $\\ensuremath{\\beta} $ are the fluid viscosity and coefficient of thermal expansion respectively. The critical Rayleigh number for this transition occurs between $\\mathit{Ra}= 1. 0\\text{{\\ndash}} 1. 2\\ensuremath{\\times} 1{0}^{10} $ for $H/ W= 2$ and $\\mathit{Ra}= 1. 2\\text{{\\ndash}} 1. 4\\ensuremath{\\times} 1{0}^{10} $ for $H/ W= 1$ , indicating that the instability has an aspect ratio dependence.
Publisher: Elsevier BV
Date: 06-2012
Publisher: Elsevier BV
Date: 04-2013
Publisher: Begell House
Date: 2015
Publisher: Author(s)
Date: 2018
DOI: 10.1063/1.5046230
No related organisations have been discovered for Michael Kirkpatrick.
Start Date: 06-2022
End Date: 06-2025
Amount: $469,000.00
Funder: Australian Research Council
View Funded ActivityStart Date: 02-2006
End Date: 12-2010
Amount: $132,400.00
Funder: Australian Research Council
View Funded ActivityStart Date: 06-2013
End Date: 12-2017
Amount: $425,000.00
Funder: Australian Research Council
View Funded ActivityStart Date: 06-2016
End Date: 12-2020
Amount: $300,000.00
Funder: Australian Research Council
View Funded ActivityStart Date: 04-2015
End Date: 12-2019
Amount: $335,100.00
Funder: Australian Research Council
View Funded ActivityStart Date: 06-2009
End Date: 12-2014
Amount: $300,000.00
Funder: Australian Research Council
View Funded ActivityStart Date: 07-2011
End Date: 06-2016
Amount: $350,000.00
Funder: Australian Research Council
View Funded ActivityStart Date: 2008
End Date: 12-2010
Amount: $450,000.00
Funder: Australian Research Council
View Funded Activity