Reducing drag under high turbulence is a critical but challenging issue that has engendered great concern. This study utilizes hydrophilic tips in superhydrophobic (SHP) grooves to enhance the stability of plastron, which results in a considerable drag reduction ($DR$) up to 62 %, at Reynolds number ($Re$) reaching $2.79 \times 10^{4}$. The effect of the spacing width $w$ of the microgrooves on both $DR$ and flow structures is investigated. Experimental results demonstrate that $DR$ increases as either microgroove spacing $w$ or $Re$ increases. The velocity fields obtained using particle image velocimetry indicate that the air-filled SHP grooves induce a considerable wall slip. This slip significantly weakens the intensity of Taylor rolls, reduces local momentum transport, and consequently lowers drag. This phenomenon becomes more pronounced with increasing $w$. Furthermore, to quantify the multiscale relationship between global response and geometrical as well as driving parameters, $DR\sim (w, \phi _s, Re)$, a theoretical model is established based on angular momentum defect theory and magnitude estimate. It is demonstrated that a decrease in the surface solid fraction can reduce wall shear, and an increase in the groove width can weaken turbulence kinetic energy production, rendering enhanced slip and drag reduction. This research has implications for designing and optimizing turbulent-drag-reducing surfaces in various engineering applications, such as transportation and marine engineering.

Geothermal gradients and heterogeneous permeability are commonly observed in natural geological formations for underground CO$_2$ sequestration. In this study, we conduct three-dimensional direct numerical simulations on the double-diffusive convection with both unstable temperature and concentration gradients in homogeneous and heterogeneous porous media. For homogeneous porous media, the root-mean-squared velocity increases linearly with density ratio defined as the buoyancy ratio by temperature and concentration differences. The flow structures show no remarkable changes when temperature Rayleigh number ${Ra}_T$ is less than its critical value, but alter from sheet-like to cellular structures as ${Ra}_T$ surpasses this threshold. The concentration wavenumber scales approximately as $k_{rS}\sim {Ra}_e^{0.47}$ with a defined effective Rayleigh number ${Ra}_e$. By using a scale analysis, the concentration flux exhibits a consistent linear relation with the total driving forces for all simulations. For heterogeneous porous media, where the Dykstra–Parsons coefficient $V_{DP}$ and correlation length $l_{r}$ determine the spatial distribution of the permeability field, the flow is strengthened in places with higher permeability. The velocity and concentration flux are less affected by $l_{r}$ than that by $V_{DP}$. For small correlation length, the flow structures coarsen and their characteristic width generally increases with increasing heterogeneity. For large correlation length, small structures emerge in the regions with large permeability, which can be attributed to the intensified local Rayleigh number triggering more vigorous convection there. The variations of concentration flux with $l_{r}$ and $V_{DP}$ can be explained by the portion of area covered by high concentration with large vertical velocity near the boundaries.

This study proposes a new mechanism that can lead to layering or convection from the finite amplitude perturbation acting on the double diffusive convection with uniform background shear. We focus on the double diffusive convection in the diffusive regime with the cold fresh water laying above the warm salty water. We demonstrate that, although the unperturbed system is linearly stable, the finite amplitude perturbation can trigger the initial flow motions which subsequently obtain energy from the gravitational potential energy and from the uniform background shear, and evolve to layering or convection. By using the linear stability analysis for the initial growth stage and the energy analysis for the following transitional stage, the critical Richardson number can be predicted theoretically. Here the Richardson number measures the relative strength of stratification to the background shear. The dominant wavenumbers and the growth rates of the corresponding modes given by linear theory agree well with the two-dimensional direct numerical simulations, and so does the critical Richardson number predicted by the theoretical model. The layering state is dominated by the double diffusion process, while the convection state at smaller Richardson number exhibits stronger influences from shear and generates smaller heat and salinity fluxes. The theoretical model is further applied to the parameter range which is relevant to the real oceanic environments and reveals that for the typical density ratio observed in the staircase regions in the Arctic Ocean, the current mechanism can lead to layering for relatively weak shear.

The present work aims at clarifying the effects of a solid boundary on the salt fingers in the wall-bounded double diffusive convection turbulence driven by the salinity and temperature differences between the top and bottom plates. The fluid properties are the same as the seawater, and two-dimensional direct numerical simulations are conducted over a wide range of the thermal and salinity Rayleigh numbers which measure the strength of driving salinity difference and stabilising temperature difference. We find that the bulk density ratio $\varLambda _b$, defined by the mean temperature and salinity gradients at the bulk, controls the flow morphology. As $\varLambda _b$ exceeds unity, the bulk flow shifts from wide convection rolls to slender salt fingers. Two different regimes are further identified for the cases of salt-finger type. One is the confined salt-finger regime where the characteristic height of salt fingers is comparable to the bulk height and the influences of the solid boundary are noticeable. The other is the free salt-finger regime where the salt fingers are much shorter than the bulk height. In this latter regime, the transport properties versus $\varLambda _b$ are in quantitative agreement with those obtained in the fully periodic domain (e.g. Traxler et al., J. Fluid Mech., vol. 677, 2011, pp. 530–553). For a limited range of density ratio at the highest salinity Rayleigh number considered here, multiple states can be obtained from different initial conditions. The large-scale secondary instability and spontaneous formation of staircase from finger layers are not observed in the current study.

Convective dissolution is an important mechanism for long-term CO$_2$ sequestration in deep saline aquifers. The presence of an unstable geothermal gradient can affect the process of dissolution. In this paper, we present direct numerical simulations in a three-dimensional porous medium at three different concentration Rayleigh numbers $Ra_S$ with a set of thermal Rayleigh numbers $Ra_T$. Simulations reveal that the flow structures alter when ${\textit {Ra}}_T$ increases for a fixed ${\textit {Ra}}_S$. Strong thermal gradient can yield large-scale convection rolls which change the horizontal distribution and motions of concentration fingers. The time evolution of fluxes also has different responses to different ${\textit {Ra}}_T$. A theoretical model is developed and successfully describes the evolution of concentration flux and volume averaged concentration during the final shutdown stage. We further calculate the dissolved CO$_2$ into the interior over time, which shows non-monotonic variations as ${\textit {Ra}}_T$ increases. At the end of our simulations, the maximum increment of dissolved CO$_2$ occurs when density ratio is around unity for all three concentration Rayleigh numbers we have explored. We apply our results to a typical geological reservoir and discuss their implications.

In this work we investigate the thermal convection driven by heat-releasing point particles. Three-dimensional direct numerical simulations are conducted for $1\times 10^7\le {\textit {Ra}} \le 1\times 10^{10}$ and $0.01\le {\textit {St}} \le 10$, where the Rayleigh number ${\textit {Ra}}$ and Stokes number ${\textit {St}}$ measure the strengths of the heat releasing rate and the Stokes drag, respectively. A regime at intermediate Stokes numbers is identified with most particles accumulating into the top boundary layer region, while for other cases particles are constantly advected over the entire domain. For the latter state, the flow motions are stronger at the upper part of the domain. The thicknesses of both momentum and thermal boundary layers at the top plate follow the same scaling law with ${\textit {Ra}}$ and show minor dependences on ${\textit {St}}$. The volume-averaged temperature and convective flux exhibit non-monotonic variations as ${\textit {St}}$ increases and reach their minimums at intermediate ${\textit {St}}$. The fraction coefficient of heat flux, i.e. the ratio between the heat flux through the bottom plate and the total flux through both plates, shares the similar dependence on ${\textit {St}}$ as the convective flux. The relation between these scaling laws can be explained by using the global balance between the dissipation and convective flux. The scaling laws for the transition between different flow regimes are also proposed and agree with the numerical results. The preferential concentration of particles is observed for all cases and is strongest at intermediate Stokes numbers, for which multiscale clustering emerges with small clusters forming larger ones.

We conduct direct numerical simulations and study the evolution of a pair of counter-rotating vortices in a stratified and turbulent environment beyond the first vortex linking (which is due to the Crow instability). The initial position of the two vortices is perpendicular to the direction of thermal stratification. We vary the Froude number, the background turbulence intensity and the Reynolds number, covering strong, weak and neutral stratification; low, moderate and high background turbulence; and low and moderate Reynolds numbers. We observe a second vortex linking in a weakly stratified environment for the first time. The second vortex linking is followed by turbulence bursts that lead to the rapid decay of the residual vortices after the first vortex linking, leading to a short-living vortex pair. This challenges the conventional view that residual vortices have a long life span, which is true only in an unstratified environment. In addition to the second vortex linking, we observe the following. First, strong background turbulence quickly breaks the two vortices. Second, the two vortices induce secondary and tertiary vortices and lose energy to gravitational waves in a strongly stratified ($Fr\leqslant 0.7$) environment.

A sequence of two- and three-dimensional simulations are conducted for the double-diffusive convection (DDC) flows in the diffusive regime subjected to an imposed shear. For a wide range of control parameters, and for sufficiently strong perturbation of the conductive initial state, staircase-like structures spontaneously develop, with relatively well-mixed layers separated by sharp interfaces of enhanced scalar gradient. Such staircases appear to be robust even in the presence of strong shear over very long times, with early-time coarsening of the observed layers. For the same set of control parameters, different asymptotic layered states, with markedly different vertical scalar fluxes, can arise for different initial perturbation structures. The imposed shear significantly spatio-temporally modifies the vertical transport of the various scalars. The flux ratio $\gamma ^*$ (i.e. the ratio between the density fluxes due to the total salt flux and the total heat flux) is found, at steady state, to be essentially equal to the square root of the ratio of the salt diffusivity to the thermal diffusivity, consistent with the physical model proposed by Linden & Shirtcliffe (J. Fluid Mech., vol. 87, 1978, pp. 417–432) and the variational arguments presented by Stern (J. Fluid Mech., vol. 114, 1982, pp. 105–121) for unsheared double-diffusive convection.

In this work we study the convection flow driven by a heat-releasing concentration field which itself is stably stratified. The heat-releasing rate is linearly proportional to concentration. Linear stability analysis is conducted to determine the critical heat-releasing rate for given fluid properties and concentration differences. The most unstable mode associated with the critical heat-releasing rate can be oscillatory for a large concentration Rayleigh number, i.e. the non-dimensionalized concentration difference and large Schmidt number, i.e. the ratio of viscosity to diffusivity of the concentration component. Fully developed flows are then investigated by direct numerical simulations. Flow structures near the bottom plate have larger horizontal scales than those near the top plate. The concentration in the bulk is almost constant and takes a similar value for all the explored parameters, which results in the convective flux increasing linearly with height. To explain the dependences of the global transport properties, we extend the unifying theory of the Rayleigh–Bénard convection to the current system and develop scaling laws for the global fluxes. The numerical results can be described accurately by the theoretical model.

In this work we show that horizontal gradients of temperature and salinity with compensating effects on density can drive thermohaline intrusion in the fluid layer below. Specifically, different types of double diffusive convection generate differential vertical fluxes from the top boundary, which then sustain horizontal temperature and salinity gradients within the bulk. Interleaving layers develop in the bulk and slope downward towards the cold fresh side, which are of the diffusive type. New layers emerge near the bottom boundary and shift the existing layers upward due to the density difference induced by the divergence of the vertical fluxes through the top surface. Detailed analyses reveal that the present intrusion is consistent with those in the narrow fronts, and both layer thickness and current velocity follow the corresponding scaling laws. Such intrusion process provides an extra path to transfer heat and salinity horizontally towards the cold and fresh side, but transfer the density anomaly towards the warm and salty side. These findings extend the circumstances where thermohaline intrusions may be observed.

Turbulent convection plays a crucial role in many natural environments and engineering applications. One of the most fundamental questions is how the heat flux depends on the thermal driving and fluid property. It has been proposed that when the fluid layer experiences extremely strong thermal driving, the boundary layers will become fully turbulent and the so-called ultimate regime will emerge. In this work we proposed a numerical experiment in which the thermal boundary layer can be spatially decoupled with the viscous one. We demonstrate that, once the thermal boundary layer is fully decoupled from the viscous boundary layer and locates entirely inside the turbulent bulk, the scaling laws corresponding to the ultimate regime can be obtained, namely $Nu \sim Ra^{1/2}Pr^{1/2}$ and $Re \sim Ra^{1/2} Pr^{-1/2}$ with $Nu$ being the Nusselt number, $Re$ the Reynolds number, $Ra$ the Rayleigh number and $Pr$ the Prandtl number, respectively. Therefore, our results support the physical conjecture of the ultimate regime for the turbulent convection.

We conduct direct numerical simulations for turbulent Rayleigh–Bénard (RB) convection, driven simultaneously by two scalar components (say, temperature and concentration) with different molecular diffusivities, and measure the respective fluxes and the Reynolds number. To account for the results, we generalize the Grossmann–Lohse theory for traditional RB convection (Grossmann & Lohse, J. Fluid Mech., vol. 407, 2000, pp. 27–56; Phys. Rev. Lett., vol. 86 (15), 2001, pp. 3316–3319; Stevens et al., J. Fluid Mech., vol. 730, 2013, pp. 295–308) to this two-scalar turbulent convection. Our numerical results suggest that the generalized theory can successfully capture the overall trends for the fluxes of two scalars and the Reynolds number without introducing any new free parameters. In fact, for most of the parameter space explored here, the theory can even predict the absolute values of the fluxes and the Reynolds number with good accuracy. The current study extends the generality of the Grossmann–Lohse theory in the area of buoyancy-driven convection flows.

Direct numerical simulations are conducted for double diffusive convection (DDC) bounded by two parallel plates. The Prandtl numbers, i.e. the ratios between the viscosity and the molecular diffusivities of scalars, are similar to the values of seawater. The DDC flow is driven by an unstable salinity difference (here across the two plates) and stabilized at the same time by a temperature difference. For these conditions the flow can be in the finger regime. We develop scaling laws for three key response parameters of the system: the non-dimensional salinity flux $\mathit{Nu}_{S}$ mainly depends on the salinity Rayleigh number $\mathit{Ra}_{S}$, which measures the strength of the salinity difference and exhibits a very weak dependence on the density ratio $\unicode[STIX]{x1D6EC}$, which is the ratio of the buoyancy forces induced by two scalar differences. The non-dimensional flow velocity $Re$ and the non-dimensional heat flux $\mathit{Nu}_{T}$ are dependent on both $\mathit{Ra}_{S}$ and $\unicode[STIX]{x1D6EC}$. However, the rescaled Reynolds number $Re\unicode[STIX]{x1D6EC}^{\unicode[STIX]{x1D6FC}_{u}^{eff}}$ and the rescaled convective heat flux $(\mathit{Nu}_{T}-1)\unicode[STIX]{x1D6EC}^{\unicode[STIX]{x1D6FC}_{T}^{eff}}$ depend only on $\mathit{Ra}_{S}$. The two exponents are dependent on the fluid properties and are determined from the numerical results as $\unicode[STIX]{x1D6FC}_{u}^{eff}=0.25\pm 0.02$ and $\unicode[STIX]{x1D6FC}_{T}^{eff}=0.75\pm 0.03$. Moreover, the behaviours of $\mathit{Nu}_{S}$ and $Re\unicode[STIX]{x1D6EC}^{\unicode[STIX]{x1D6FC}_{u}^{eff}}$ agree with the predictions of the Grossmann–Lohse theory which was originally developed for the Rayleigh–Bénard flow. The non-dimensional salt-finger width and the thickness of the velocity boundary layers, after being rescaled by $\unicode[STIX]{x1D6EC}^{\unicode[STIX]{x1D6FC}_{u}^{eff}/2}$, collapse and obey a similar power-law scaling relation with $\mathit{Ra}_{S}$. When $\mathit{Ra}_{S}$ is large enough, salt fingers do not extend from one plate to the other and horizontal zonal flows emerge in the bulk region. We then show that the current scaling strategy can be successfully applied to the experimental results of a heat–copper–ion system (Hage & Tilgner, Phys. Fluids, vol. 22, 2010, 076603). The fluid has different properties and the exponent $\unicode[STIX]{x1D6FC}_{u}^{eff}$ takes a different value $0.54\pm 0.10$.

We investigate how compressibility affects the turbulent statistics from a Lagrangian point of view, particularly in the parameter range where the flow transits from the incompressible type to a state dominated by shocklets. A series of three-dimensional simulations were conducted for different types of driving and several Mach numbers. For purely solenoidal driving, as the Mach number increases a new self-similar region first emerges in the Lagrangian structure functions at sub-Kolmogorov time scale and gradually extends to larger time scale. In this region the relative scaling exponent saturates and the saturated value decreases as the compressibility becomes stronger, which can be attributed to the shocklets. The scaling exponent for the inertial range is still very close to that of incompressible turbulence for small Mach number, and discrepancy becomes visible when the Mach number is high enough. When the driving force is dominated by the compressive component the shocklet-induced self-similar region occupies a much wider range of time scales than that in the purely solenoidal driving case. Regardless of the type of driving force, the probability density functions of the velocity increment collapse onto one another for the time scales in the new self-similar region after proper normalization.

The double diffusive convection between two parallel plates is numerically studied for a series of parameters. The flow is driven by the salinity difference and stabilised by the thermal field. Our simulations are directly compared with experiments by Hage & Tilgner (Phys. Fluids, vol. 22, 2010, 076603) for several sets of parameters and reasonable agreement is found. This, in particular, holds for the salinity flux and its dependence on the salinity Rayleigh number. Salt fingers are present in all simulations and extend through the entire height. The thermal Rayleigh number seems to have a minor influence on the salinity flux but affects the Reynolds number and the morphology of the flow. In addition to the numerical calculation, we apply the Grossmann–Lohse theory for Rayleigh–Bénard flow to the present problem without introducing any new coefficients. The theory successfully predicts the salinity flux both with respect to the scaling and even with respect to the absolute value for the numerical and experimental results.

In the traditional hybrid RANS/LES approaches for the simulation of wall-bounded fluid turbulence, such as detached-eddy simulation (DES), the whole flow domain is divided into an inner layer and an outer layer. Typically the Reynolds-averaged Navier–Stokes (RANS) equations are used for the inner layer, while large-eddy simulation (LES) is used for the outer layer. The transition from the inner-layer solution to the outer-layer solution is often problematic due to the lack of small-scale dynamics in the RANS region. In this paper, we propose to simulate the whole flow domain by large-eddy simulation while enforcing a Reynolds-stress constraint on the subgrid-scale (SGS) stress model in the inner layer. Both the algebraic eddy-viscosity model and the one-equation Spalart–Allmaras (SA) model have been used to constrain the Reynolds stress in the inner layer. In this way, we improve the LES methodology by allowing the mean flow of the inner layer to satisfy the RANS solution while small-scale dynamics is included. We validate the Reynolds-stress-constrained large-eddy simulation (RSC-LES) model by simulating three-dimensional turbulent channel flow and flow past a circular cylinder. Our model is able to predict mean velocity, turbulent stress and skin-friction coefficients more accurately in turbulent channel flow and to estimate the pressure coefficient after separation more precisely in flow past a circular cylinder compared with the pure dynamic Smagorinsky model (DSM) and DES using the same grid resolution. Furthermore, the computational cost of the RSC-LES is almost the same as that of DES.