The dimensional transition in turbulent jets of a shear-thinning fluid is studied via direct numerical simulations. Our findings reveal that under vertical confinement, the flow exhibits a unique mixed-dimensional (or 2.5-dimensional) state, where large-scale two-dimensional and small-scale three-dimensional structures coexist. This transition from three-dimensional turbulence near the inlet to two-dimensional dynamics downstream is dictated by the level of confinement: weak confinement guarantees turbulence to remain three-dimensional, whereas strong confinement forces the transition to two dimensions; the mixed-dimensional state is observed for moderate confinement and it emerges as soon as flow scales are larger than the vertical length. In this scenario, we observed that the mixed-dimensional state is an overall more energetic state, and it shows a multi-cascade process, where the direct cascade of energy at small scales and the direct cascade of enstrophy at large scales coexist. The results provide insights into the complex dynamics of confined turbulent flows, relevant in both natural and industrial settings.

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan (J. Fluid Mech., vol. 974, 2023, A35) and Jacobs & Dalziel (J. Fluid Mech., vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

The research on elasto-inertial turbulence (EIT), a new type of turbulent flow, has reached the stage of identifying the minimal flow unit (MFU). On this issue, direct numerical simulations of FENE-P fluid flow in two-dimensional channels with variable sizes are conducted in this study. We demonstrate with the increase of channel length that the simulated flow experiences several different flow patterns, and there exists an MFU for EIT to be self-sustained. At Weissenberg number ($Wi$) higher than the one required to excite EIT, when the channel length is relatively small, a steady arrowhead regime (SAR) flow structure and a laminar-like friction coefficient is achieved. However, as the channel length increases, the flow can fully develop into EIT characterized with high flow drag. Close to the size of the MFU, the simulated flow behaves intermittently between the SAR state with low drag and EIT state with high drag. The flow falling back to ‘laminar flow’ is caused by the insufficient channel size below the MFU. Furthermore, we give the relationship between the value of the MFU and the effective $Wi$, and explain its physical reasons. Moreover, the intermittent flow regime obtained based on the MFU gives us an opportunity to look into the origin and exciting process of EIT. Through capturing the onset process of EIT, we observed that EIT originates from the sheet-like extension structure located near the wall, which is maybe related to the wall mode rather than the centre mode. The fracture and regeneration of this sheet-like structure is the key mechanism for the self-sustaining of EIT.

We report direct numerical simulations results of the rough-wall channel, focusing on roughness with high $k_{rms}/k_a$ statistics but small to negative $Sk$ statistics, and we study the implications of this new dataset on rough-wall modelling. Here, $k_{rms}$ is the root mean square, $k_a$ is the first-order moment of roughness height, and $Sk$ is the skewness. The effects of packing density, skewness and arrangement of roughness elements on mean streamwise velocity, equivalent roughness height ($z_0$) and Reynolds and dispersive stresses have been studied. We demonstrate that two-point correlation lengths of roughness height statistics play an important role in characterizing rough surfaces with identical moments of roughness height but different arrangements of roughness elements. Analysis of the present as well as historical data suggests that the task of rough-wall modelling is to identify geometric parameters that distinguish the rough surfaces within the calibration dataset. We demonstrate a novel feature selection procedure to determine these parameters. Further, since there is no finite set of roughness statistics that distinguish between all rough surfaces, we argue that obtaining a universal rough-wall model for making equivalent sand-grain roughness ($k_s$) predictions would be challenging, and that each rough-wall model would have its applicable range. This motivates the development of group-based rough-wall models. The applicability of multi-variate polynomial regression and feedforward neural networks for building such group-based rough-wall models using the selected features has been shown.

The present study offers a twofold contribution on counter-gradient transport (CGT) of turbulent scalar flux. First, by examining turbulent scalar mixing through synchronized particle image velocimetry and planar laser-induced fluorescence on an inclined jet in cross-flow, we clarify the previously unexplained phenomenon of CGT, revealing key flow structures, their spatial distribution and modelling implications. Statistical analysis identifies two distinct CGT regions: local cross-gradient transport in the windward shear layer and non-local effects near the wall after injection. These behaviours are driven by specific flow structures, namely Kelvin–Helmholtz vortices (local) and wake vortices (non-local), suggesting that scalar flux can be decomposed into a gradient-type term for gradient diffusion and a term for large-eddy stirring. Second, we propose a new approach for reconstruction of turbulent mean flow and scalar fields using continuous adjoint data assimilation (DA). By rectifying model-form errors through anisotropic correction under observational constraints, our DA model minimizes discrepancies between experimental measurements and numerical predictions. As expected, the introduced forcing term effectively identifies regions where traditional models fall short, particularly in the jet centreline and near-wall regions, thereby enhancing the accuracy of the mean scalar field. These enhancements occur not only within the observation region but also in unseen regions, underscoring present DA approach's reliability and practicality for reproducing mean flow behaviours from limited data. These findings lay a solid foundation for adjoint-based model-consistent data-driven methods, offering promising potential for accurately predicting complex flow scenarios like film cooling.

Direct numerical simulations of spanwise-rotating turbulent channel flow with a parabolic bump on the bottom wall are employed to investigate the effects of rotation on flow separation. Four rotation rates, $Ro_b := 2\varOmega H/U_b = \pm 0.42$, $\pm$1.0, are compared with the non-rotating scenario. The mild adverse pressure gradient induced by the lee side of the bump allows for a variable pressure-induced separation. The separation region is reduced (increased) when the bump is on the anti-cyclonic (cyclonic) side of the channel, compared with the non-rotating separation. The total drag is reduced in all rotating cases. Through several mechanisms, rotation alters the onset of separation, reattachment and wake recovery. The mean momentum deficit is found to be the key. A physical interpretation of the ratio between the system rotation and mean shear vorticity, $S:=\varOmega /\varOmega _s$, provides the mechanisms regarding stability thresholds $S=-0.5$ and $-$1. The rotation effects are explained accordingly, with reference to the dynamics of several flow structures. For anti-cyclonic separation, particularly, the interaction between the Taylor–Görtler vortices and hairpin vortices of wall-bounded turbulence is proven to be responsible for the breakdown of the separating shear layer. A generalized argument is made regarding the essential role of near-wall deceleration and resultant ejection of enhanced hairpin vortices in destabilizing an anti-cyclonic flow. This mechanism is anticipated to have broad impacts on other applications in analogy to rotating shear flows, such as thermal convection and boundary layers over concave walls.

A permeable disk serves as a simplified model for the conversion of wind energy by a horizontal axis wind turbine. In this study, we investigate how inflow turbulence intensity (TI), $I_\infty$, and inflow turbulence integral length scale, $L_\infty$, influence the flow recovery in the wake, the capability of a permeable disk in extracting turbulence kinetic energy (TKE) of the incoming flow, and the statistics of wake-added turbulence using large-eddy simulation. The simulated inflows include various TIs (i.e. $I_\infty =2.5\,\%$–$25\,\%$) and integral length scales (i.e. $L_\infty / D =0.5$–$2.0$) for two thrust coefficients. Simulation results show that both inflow TI and integral length scale influence flow recovery via enhanced ejections and sweeps across the wake boundary, with the former strongly affecting the position where the wake starts to recover and the latter mainly on the recovery rate. Moreover, it is shown that increasing $I_\infty$ and $L_\infty$ increases the TKE extraction by the disk, occurring mainly at scales ($s$) greater than $0.5D$ and frequencies depending on the inflow integral length scale. As for the wake-added TKE, the inflow TI mainly affects its intensity, while the inflow integral length scale affects both its intensity and the sensitive frequencies, with the spectral distributions in scale space ($s$) being similar and the peak located around $s/D=1.0$ for the considered inflows.

Turbulent flows over rough surfaces can be encountered in a wide range of engineering applications. Despite the progress made after several decades of studies, the prediction of drag and roughness function from the surface geometrical parameters remains an open question. Several methods have shown encouraging results. However, they lack generality and present some scatter in the data. In this paper we propose a new parameter, the effective distribution ($ED$), which lays foundation on the effective slope with some changes to take into account the sheltering effect of large roughness elements and the drag induced by pinnacles higher than the average roughness elements. To develop this new correlation between geometrical features of the wall and the drag, we performed a set of simulations of the turbulent flow over a rough surface made of triangular elements varying their height and spatial distribution. The $ED$ correlates quite well both with the drag and the roughness function for a wide range of cases having different mean roughness height, skewness and kurtosis. To further validate the $ED$, and assessing how it can be generalized to real rough wall, an irregular wall made from the superposition of random sinusoidal function was considered. Results were consistent with the correlation here presented.

We present direct numerical simulation (DNS) and modelling of incompressible, turbulent, generalized Couette–Poiseuille flow. A particular example is specified by spherical coordinates $(Re,\theta,\phi )$, where $Re = 6000$ is a global Reynolds number, $\phi$ denotes the angle between the moving plate, velocity-difference vector and the volume-flow vector and $\tan \theta$ specifies the ratio of the mean volume-flow speed to the plate speed. The limits $\phi \to 0^\circ$ and $\phi \to 90^\circ$ give alignment and orthogonality, respectively, while $\theta \to 0^\circ,\ \theta \to 90^\circ$ correspond respectively to pure Couette flow in the $x$ direction and pure Poiseuille flow at angle $\phi$ to the $x$ axis. Competition between the Couette-flow shear and the forced volume flow produces a mean-velocity profile with directional twist between the confining walls. Resultant mean-speed profiles relative to each wall generally show a log-like region. An empirical flow model is constructed based on component log and log-wake velocity profiles relative to the two walls. This gives predictions of four independent components of shear stress and also mean-velocity profiles as functions of $(Re,\theta,\phi )$. The model captures DNS results including the mean-flow twist. Premultiplied energy spectra are obtained for symmetric flows with $\phi =90^\circ$. With increasing $\theta$, the energy peak gradually moves in the direction of increasing $k_x$ and decreasing $k_z$. Rotation of the energy spectrum produced by the faster moving velocity near the wall is also observed. Rapid weakening of a spike maxima in the Couette-type flow regime indicates attenuation of large-scale roll structures, which is also shown in the $Q$-criterion visualization of a three-dimensional time-averaged flow.

High-fidelity simulations are conducted to investigate the turbulent boundary layers around a finite-span NACA0012 wing with a rounded wing-tip geometry at a chord-based Reynolds number of $Re_c=200\,000$ and at various angles of attack up to $10^\circ$. The study aims to discern the differences between the boundary layers on the finite-span wing and those on infinite-span wings at equivalent angles of attack. The finite-span boundary layers exhibit: (i) an altered streamwise and a non-zero spanwise pressure gradient as a result of the variable downwash induced by the wing-tip vortices (an inviscid effect typical of finite-span wings); (ii) differences in the flow history at different wall-normal distances, caused by the variable flow angle in the wall-normal direction (due to constant pressure gradients and variable momentum normal to the wall); (iii) laminar flow entrainment into the turbulent boundary layers near the wing tip (due to a laminar–turbulent interface); and (iv) variations in boundary layer thickness across the span, attributed to the variable wall-normal velocity in that direction (a primarily inviscid effect). These physical effects are then used to explain the differences in the Reynolds stress profiles and other boundary layer quantities, including the reduced near-wall peak of the streamwise Reynolds stress and the elevated Reynolds stress levels near the boundary layer edge, both observed in the finite-span wings. Other aspects of the flow, such as the downstream development of wing-tip vortices and their interactions with the surrounding flow, are reserved for future investigations.

The first two parts of this study showed by experiment and by large-eddy simulation that, in mixing layers formed of processions of spanwise-oriented vortex structures, there is a fundamental change in the dynamics of the large-structure growth at what has been called the ‘mixing’ transition. This third part examines the relationship of the post-transition ‘organised’ flow to the mixing layer of classical turbulence theory. Further large-eddy simulations are presented which, like some of the experiments reported in Part 1, have captured the mixing layer in both its classical and organised turbulent states, thereby allowing them to be characterised separately. The simulation results are then used to calibrate model-free integral analyses of the dynamics of the mixing layer's spatial growth. It is shown that the organised and classical flows are alternative self-preserving turbulent flow states involving fundamentally different exchanges of mass, momentum and energy, one or other of which emerges naturally in particular realisations depending on the initial conditions.

Two-dimensional (2-D) quadrant analysis is generally used for investigating flow and sediment dynamics around a rigid structure in open channel flows. Given that particle distribution around rigid obstacles is not spatially uniform and changes in time, while vortices evolve to become three-dimensional (3-D) structures, 2-D quadrant analysis might be unsuitable to completely determine the sediment transport. Hence, 3-D quadrant and 3-D octant analyses should be considered, using the 3-D instantaneous velocity data and relative 3-D bursting process to define sediment transport surrounding the submerged square and circular cylinders. The turbulent kinetic energy (TKE), transition probabilities, occurrence probabilities, stress fraction and angles of inclination of 3-D bursting events are considered to quantify the coherent structures surrounding the cylinders and their interaction with bed particles. Experiments were conducted at the Hydraulics and Water Resources Engineering Laboratory, School of Infrastructure, Indian Institute of Technology, Bhubaneswar, Odisha, India, and velocity data were recorded at different cross-sections around the submerged cylinders using an acoustic Doppler velocimeter. Results show that the TKE is greater for internal ejection, external ejection, and internal sweep, external sweep in the upstream of the circular and square cylindrical structures. On the other part, the TKE is significant for internal ejection, external ejection, and internal sweep, external sweep in the downstream of the aligned square cylindrical structures, which justifies the highest scour depth that occurred upstream of the circular and square cylindrical structures and downstream of the aligned square cylindrical structure. The transition probability of the bursting events was determined using the Markov process from the measured velocity data to investigate the consecutive occurrence of bursting events. Further, the importance of sweeps and ejections on sediment erosion surrounding the cylinders within the scour hole at various stages of its development was investigated via 3-D quadrant analysis of the bursting occurrences. The outcomes show that external sweep and internal ejection events are active mechanisms for bed particle transport surrounding the cylinder. The maximum transition probability values are found around aligned cylindrical structures in comparison with the circular and square cylindrical structures in the transverse direction. This depicts the formation of a trailing vortex on both sides of the aligned square cylindrical object. The results reveal that the effect of inclination angles with respect to the water flow is greater for internal ejection and external sweep from upstream to downstream within the scour hole surrounding the cylindrical structures at various phases of development as horseshoe vortices and downflow develop upstream of the cylindrical structures while trailing vortices and wake vortices form at the top and downstream of the cylindrical structures. Internal and external ejection have a higher stress fraction than an internal and external sweep for square cylinders with alignment angles of 0°, 20° and circular cylinders over underdeveloped and developed scoured beds, respectively. With the higher percentage of fractional contributions for internal sweeps, the external sweep is predicted close to the cylindrical objects in comparison with the internal ejection and external ejection events because of the formation and warping of the horseshoe vortex close to the cylindrical objects, suggesting a significant probability of 3-D bursting occurrences with sediment movement near the cylindrical structures.

The Reynolds-averaged Navier–Stokes (RANS) models depend on empirical constants to close the Reynolds stress terms. The empirical constants were obtained using experiments conducted at low Reynolds numbers several decades ago. In this paper, we revisit the turbulent viscosity parameter $C_\mu$, based on the stress–intensity ratio $c^2 = {|\overline {uw}|}/{k}$. Here, $\overline {|uw|}$ and $k$ are the absolute values of the Reynolds stress and turbulent kinetic energy, respectively. Through a priori comparisons, we find that the currently accepted value of $C_\mu = 0.09$ does not agree with the latest direct numerical simulation (DNS) and experimental datasets of wall-bounded turbulent planar flows. Therefore, a new value is suggested by averaging $c^2$ in the equilibrium region, where the production ($\mathcal {P}$) of $k$ is within 10 % of the dissipation rate ($\epsilon$), and consequently, $c^4 \approx C_\mu$. We evaluate flows up to friction Reynolds number $Re_\tau \approx 10\,000$ and find that with increasing $Re_\tau$, $C_\mu$ approaches a value of 0.06, which is almost 50 % lower than the prevalent value of 0.09. Finally, we perform an a priori test with the new (proposed) value of $C_\mu = 0.06$ to show that the estimated turbulent viscosity $\nu _T$ for wall-bounded flows is in much closer agreement with the exact (DNS) values than when $\nu _T$ is estimated using $C_\mu = 0.09$.

Large-scale coherent structures in incompressible turbulent pipe flow are studied for a wide range of Reynolds numbers ($Re_\tau =180, 550, 1000, 2000$ and $5200$). Employing the Karhunen–Loève decomposition and a novel approach based on the Voronoi diagram, we identify and classify statistically coherent structures based on their location, dimensions and $Re_{\tau }$. With increasing $Re_{\tau }$, two distinct classes of structures become more energetic, namely wall-attached and detached eddies. The Voronoi methodology is shown to delineate these two classes without the need for specific criteria or thresholds. At the highest $Re_{\tau }$, the attached eddies scale linearly with the wall-normal distance with a slope of approximately $l_y\sim 1.2y/R$, while the detached eddies remain constant at the size of $l_y \approx 0.26R$, with a progressive shift towards the pipe centre. We extract these two classes of structures and describe their spatial characteristics, including radial size, helix angle and azimuthal self-similarity. The spatial distribution could help explain the differences in mean velocity between pipe and channel flows, as well as in modelling large and very-large-scale motions (LSM and VLSM). In addition, a comprehensive description is provided for both wall-attached and detached structures in terms of LSM and VLSM.

We use three-dimensional direct numerical simulations of homogeneous isotropic turbulence in a cubic domain to investigate the dynamics of heavy, chiral, finite-size inertial particles and their effects on the flow. Using an immersed-boundary method and a complex collision model, four-way coupled simulations have been performed, and the effects of particle-to-fluid density ratio, turbulence strength and particle volume fraction have been analysed. We find that freely falling particles on the one hand add energy to the turbulent flow but, on the other hand, they also enhance the flow dissipation: depending on the combination of flow parameters, the former or the latter mechanism prevails, thus yielding enhanced or weakened turbulence. Furthermore, particle chirality entails a preferential angular velocity which induces a net vorticity in the fluid phase. As turbulence strengthens, the energy introduced by the falling particles becomes less relevant and stronger velocity fluctuations alter the solid phase dynamics, making the effect of chirality irrelevant for the large-scale features of the flow. Moreover, comparing the time history of collision events for chiral particles and spheres (at the same volume fraction) suggests that the former tend to entangle, in contrast to the latter which rebound impulsively.

A high-resolution numerical simulation of an air–water turbulent upward bubbly flow in a pipe is performed to investigate the turbulence characteristics and bubble interaction with the wall. We consider three bubble equivalent diameters and three total bubble volume fractions. The bulk and bubble Reynolds numbers are $Re_{bulk}= u_{bulk} D/\nu _w = 5300$ and $Re_{bub}= (\langle u_{bub}\rangle - u_{bulk}) d_{eq}/\nu _w = 533\unicode{x2013}1000$, respectively, where $u_{bulk}$ is the water bulk velocity, $\langle u_{bub}\rangle$ is the overall bubble mean velocity, $D$ is the pipe diameter and $\nu _w$ is the water kinematic viscosity. The mean water velocity near the wall significantly increases due to bubble interaction with the wall, and the root-mean-square water velocity fluctuations are proportional to $\bar {\psi }(r)^{0.4}$, where $\bar {\psi } (r)$ is the mean bubble volume fraction. For the cases considered, the bubble-induced turbulence suppresses the shear-induced turbulence and becomes the dominant flow characteristic at all radial locations including near the wall. Rising bubbles near the wall mostly bounce against the wall rather than slide along the wall or hang around the wall without collision. Low-speed streaks observed in the near-wall region in the absence of bubbles nearly disappear due to the bouncing bubbles. These bouncing bubbles generate counter-rotating vortices in their wake, and increase the skin friction by sweeping high-speed water towards the wall. We also suggest an algebraic Reynolds-averaged Navier–Stokes model considering the interaction between shear-induced and bubble-induced turbulence. This model provides accurate predictions for a wide range of liquid bulk Reynolds numbers.

The influence of outer large-scale motions (LSMs) on near-wall structures in compressible turbulent channel flows is investigated. To separate the compressibility effects, velocity fluctuations are decomposed into solenoidal and dilatational components using the Helmholtz decomposition method. Solenoidal velocity fluctuations manifest as near-wall streaks and outer large-scale structures. The spanwise drifting of near-wall solenoidal streaks is found to be driven by the outer LSMs, while LSMs have a trivial influence on the spanwise density of solenoidal streaks, consistent with the outer LSM impacts found in incompressible flows (Zhou et al., J. Fluid Mech., vol. 940, 2022, p. A23). Dilatational motions are characterized by the near-wall small-scale travelling-wave packets and the large-scale parts in the outer region. The streamwise advection velocity of the near-wall structures remains at $16 \sim 18u_{\tau }$, hardly influenced by Mach numbers, Reynolds numbers and wall temperatures. The spanwise drifting of near-wall dilatational structures, quantified by the particle image velocimetry method, follows a mechanism distinct from solenoidal streaks. This drifting velocity is notably larger than those of the solenoidal streaks, and the influence of outer LSMs is not the primary trigger for this drifting.

The influence of symmetry-breaking effects of ridge-type roughness on secondary currents in turbulent channel flow is investigated using direct numerical simulations. The ridges have triangular cross-section, which is systematically varied from isosceles to right-angled triangle, introducing an imbalance to the slopes of the ridges’ lateral surfaces while the streamwise homogeneity of the surfaces is maintained. In all cases, secondary current vortices are produced, but asymmetric ridge cross-sections break the symmetry of these vortices. As a result of the asymmetry-induced misalignment and imbalance in the secondary current vortices, net spanwise flow emerges. The magnitude of the spanwise flow increases with the slope ratio of the ridge lateral surfaces and significantly modifies the mean flow topology, leading to the merging of critical points in the case of the right-angled triangular ridge shape. Within the cavities, the net spanwise flow is accompanied by a non-zero mean spanwise pressure gradient, while from the perspective of the outer flow, the scalene ridge surfaces have a similar effect as a wall that is slowly moving in the spanwise direction. Overall, the present results suggest the existence of a special type of Prandtl's secondary currents of the second kind, namely those that result in net spanwise flow.

Comprehensive coherent structures around a surface-mounted low aspect ratio square cylinder in uniform flow with an oblique angle of $45^{\circ }$ were investigated for cylinder-width-based Reynolds numbers of 3000 and 10 000 by direct numerical simulation based on a topology-confined mesh refinement framework. High-resolution simulations and the critical-point concept were scrutinized to reveal for the first time the reasonable and compatible topologies of flow separation and complete near-wall structures, due to their extensive impact on various engineering applications. Large-scale horseshoe vortices are observed at two notable foci in the viscous sublayer. Within this layer, a wall-parallel jet is formed by downflow intruding into the bottom surface at the half-saddle point, then deflecting in the upstream direction and finally penetrating the bottom surface until the half-saddle point. A pair of conical vortices on the cylinder's top surface switch themselves on two sides of the square cylinder, where the switching frequency is identical with that of the sway of the side shear layer. The undulation of the Kelvin–Helmholtz instability is identified in the instantaneous development of a conical vortex and side shear layer, where the scaling of the ratio of the Kelvin–Helmholtz and von Kármán frequencies follows the power-law relation obtained by Lander et al. (J. Fluid Mech., vol. 849, 2018, pp. 1096–1119). Large-scale arch-shaped vortex is often detected in the intermediate wake region of a square cylinder, involving two interconnected portions, such as the leg portion separated from leeward surfaces and head portion rolled up from the top surface. The leg portion of the arch-shaped vortex was rooted by two foci near the bottom-surface plane.