The elliptic approximation (EA) – rooted in Taylor’s frozen flow hypothesis, Kolmogorov’s theory of small-scale turbulence, and the Kraichnan–Tennekes random sweeping hypothesis – remains a foundational framework for modelling spatiotemporal velocity correlations in incompressible wall-bounded turbulence. This study revisits the model’s theoretical basis, and extends its applicability to velocity and temperature fluctuations in supersonic channel flows. First, we identify non-elliptic distortions in the viscous sublayer, and introduce a shear-induced acceleration that captures the observed deviation from the assumed constant convection velocity at large time separations. Next, we show that the inertial-range scalings underpinning the EA are not valid in regions where the model remains accurate; instead, its validity is supported by extended self-similarity between spatial and temporal structure functions. Finally, we conduct high-fidelity direct numerical simulations of compressible channel flows with fluctuating Mach numbers up to 0.8; our data confirm the robustness of the EA under supersonic conditions, and its effectiveness in characterising both velocity and temperature correlations. Together, these findings provide new theoretical insights into the spatiotemporal structure of wall-bounded turbulence, and broaden the operational envelope of the EA.

We consider the steady heat transfer between a collection of impermeable obstacles immersed in an incompressible two-dimensional (2-D) potential flow, when each obstacle has a prescribed boundary temperature distribution. Inside the fluid, the temperature satisfies a variable-coefficient elliptic partial differential equation (PDE), the solution of which usually requires expensive techniques. To solve this problem efficiently, we construct multiply connected conformal maps under which both the domain and governing equation are greatly simplified. In particular, each obstacle is mapped to a horizontal slit and the governing equation becomes a constant-coefficient elliptic PDE. We then develop a boundary integral approach in the mapped domain to solve for the temperature field when arbitrary Dirichlet temperature data are specified on the obstacles. The inverse conformal map is then used to compute the temperature field in the physical domain. We construct our multiply connected conformal maps by exploiting the flexible and highly accurate AAA-LS algorithm. In multiply connected domains and domains with non-constant boundary temperature data, we note similarities and key differences in the temperature fields and Nusselt number scalings as compared with the isothermal simply connected problem analysed by Choi et al. (J. Fluid Mech., vol. 536, 2005, pp. 155–184). In particular, we derive new asymptotic expressions for the Nusselt number in the case of arbitrary non-constant temperature data in singly connected domains at low Péclet number, and verify these scalings numerically. While our language focuses on the problem of conjugate heat transfer (the transfer of heat between objects in a flow), our methods and findings are equally applicable to the advection–diffusion of any passive scalar in a potential flow.

Compressible wall-bounded turbulent flows exhibit complex mean profiles because of the pronounced compressibility effects and heat transfer. We propose a hybrid transformation framework to collapse compressible mean velocity and temperature profiles onto incompressible forms through scaling each layer by its effective transformation, with the underlying mapping functions discovered via a physics-informed symbolic regression (PISR) method. The hybrid velocity transformation incorporates an intrinsic compressibility correction for the buffer layer and a PISR-derived mapping function for the logarithmic layer. For temperature, we introduce a hybrid transformation that integrates the Mach-invariant-type transformation in the viscous sublayer and a novel PISR-derived scaling in the logarithmic layer. The performance of these transformations is evaluated across compressible turbulent boundary layers with free-stream Mach numbers ranging from 0.5 to 8 and wall-to-recovery-temperature ratios ranging from 0.25 to 1. The hybrid velocity transformation outperforms Griffin–Fu–Moin transformation for the transformed mean velocity profiles, with the mean integrated percent error across the dataset decreasing from 1.67 % to 0.96 %. The hybrid temperature transformation performs better than the Mach-invariant-type and Trettel–Larsson-type transformations for mean temperature profiles. Moreover, the inverse hybrid velocity and temperature transformations can effectively predict the compressible mean velocity and temperature profiles with only wall conditions.

Although stably stratified shear flows, where the base velocity shear is quasi-continuously forced externally, arise in many geophysically and environmentally relevant circumstances, the emergent dynamics of their ensuing statistically steady stratified turbulence is still an open question. We address this phenomenon in a series of three-dimensional direct numerical simulations using spectral element methods. We consider a forced, stably stratified shear flow with an initial bulk Reynolds number $\textit{Re}_{0} = 50$, an initial bulk Richardson number $\textit{Ri}_{0} = 1/80$ (also corresponding to the initial minimum gradient Richardson number $\textit{Ri}_{{g}}$) and a fluid of Prandtl number ${\textit{Pr}} = 1$ in horizontally extended domains. Although the initial configuration is unstable to a primary Kelvin–Helmholtz instability, the ensuing turbulence is sustained by continuously relaxing the resulting flow back towards the initial profiles of streamwise velocity and buoyancy. We study statistical as well as structural aspects of the final statistically steady flows, including the flux coefficient $\varGamma _{\chi }$ and dynamically emergent length scales $\varLambda$ associated with the large-scale dynamics, respectively. Despite the ongoing stirring and mixing, we find that the shear layer half-depth converges to a finite value of $d \approx 8$ (i.e. $\varLambda _{z} \approx 16$) once the horizontal extent of the domain $L_{{h}} \gtrsim 96$. While this implies a final ${{Re}} \approx 400$ and ${Ri} \approx 0.1$, we hypothesise that such forced flows ‘tune’ themselves eventually to a state of a gradient Richardson number $\textit{Ri}_{{g}} \lesssim 0.2$, consistently with several previous studies. Moreover, provided sufficiently extended domains, we observe the emergence of large-scale flow structures with spanwise $\varLambda _{\!y} \approx 50$ and streamwise $\varLambda _{x} \lesssim 115$. Clearly, these observations demonstrate the marked anisotropy of characteristic emergent length scales, even for such ‘weakly stratified’ forced shear flows. We conjecture that the actual emergent streamwise structures are a vestigial ‘imprint’ in the sheared turbulent flow of the primary linear instability of the converged deepened turbulent shear layer.

In this study, we investigate the dynamic behaviour of reconfigurable circular plates under acceleration as a model problem to understand the interplay between kinematics and shape deformation in biological propulsion. A high-resolution force transducer and time-resolved particle image velocimetry were employed to simultaneously capture both hydrodynamic forces and vortex dynamics. The results reveal that, unlike rigid plates that exhibit Reynolds number independence, the force evolution of reconfigurable plates is governed by the dimensionless bending stiffness ${\textit{EI}}^*$. A distinct load-shifting phenomenon is observed – characterized by a reduction in peak force amplitude and an elevation of the postpeak force trough, contrasting with the ‘peak-valley’ behaviour typical of rigid plates. Based on ${\textit{EI}}^*$, reconfigurable plates are classified into three regimes: extra-flexible (${\textit{EI}}^* \lt 2.28 \times 10^{-3}$), flexible ($2.28 \times 10^{-3} \leqslant {\textit{EI}}^* \leqslant 0.143$) and rigid (${\textit{EI}}^* \gt 0.143$). Notably, only plates within the flexible regime exhibit the load-shifting phenomenon. Flow visualizations show that the flexible plates, due to their shape reconfiguration, produce flow fields with two distinct features: initially, the formation of three-dimensional, non-axisymmetric vortex rings; subsequently, vortex breakdown occurs due to instability. By applying the vorticity moment theorem, force generation is accurately estimated from the flow field. Using a vortex-based low-order force model, the radial distribution of vorticity is identified as the key mechanism underlying the load-shifting phenomenon. This finding suggests that biological morphing structures in real propulsion scenarios can reduce force fluctuations without compromising average thrust by ‘load-shifting’, offering insights into efficient propulsion strategies.

We report experimental evidence of an Eulerian-mean flow, $\overline {u}(z)$, created by the interaction of surface waves and tailored ambient sub-surface turbulence, which partly cancels the Stokes drift, $u_s(z)$, and present supporting theory. Water-side turbulent velocity fields and Eulerian-mean flows were measured with particle image velocimetry before vs after the passage of a wave group, and with vs without the presence of regular waves. We compare different wavelengths, steepnesses and turbulent intensities. In all cases, a significant change in the Eulerian-mean current is observed, strongly focused near the surface, where it opposes the Stokes drift. The observations support the picture that, when waves encounter ambient sub-surface turbulence, the flow undergoes a transition during which Eulerian-mean momentum is redistributed vertically (without changing the depth-integrated mass transport) until a new equilibrium state is reached, wherein the near-surface ratio between $|{\rm d}\overline {u}/{\rm d}z|$ and $|{\rm d}u_s/{\rm d} z|$ approximately equals the ratio between the streamwise and vertical Reynolds normal stresses. This accords with a simple statistical theory derived here and holds regardless of the absolute turbulence level, whereas stronger turbulence means faster growth of the Eulerian-mean current. We present a model based on Rapid distortion theory which describes the generation of the Eulerian-mean flow as a consequence of the action of the Stokes drift on the background turbulence. Predictions are in qualitative, and reasonable quantitative, agreement with experiments on wave groups, where equilibrium has not yet been reached. Our results could have substantial consequences for predicting the transport of water-borne material in the oceans.

In this study, direct numerical simulation of a turbulent flame–wall interaction (FWI) has been done for premixed H$_2/$air and NH$_3/$H$_2/$air flames in a fully developed channel flow at Re$_\tau$ $\approx$ 300. Both isothermal and adiabatic walls are considered. The results contribute to further clarification of the underlying mechanisms of FWIs. First, the underlying mechanism for the rapid increase of chemical flame thickness near the wall is found to be the zero-flux boundary condition for diffusion. Effects of wall heat loss and wall turbulence are minor. Then, a ridge-based flame surface identification method is proposed to track the flame front, which is found to be more accurate than an isosurface of $C$ (the progress variable), especially during FWIs. Using this technique, the near-wall flame geometry and orientation are correctly captured. It is found that the flames are laminarised near the wall and almost parallel to the isothermal wall shortly before quenching. Flame–vortex interactions lead to entrained flame pockets for H$_2$ as a fuel and to a distributed reaction zone for the case of NH$_3/$H$_2$. Finally, the turbulent combustion regime is investigated by checking wall-distance-dependent Reynolds number and Karlovitz number. It is found that the flames enter the laminar flame regime shortly before wall quenching, instead of the broken reaction regime suggested in previous studies. To support the analysis, the turbulent flame dynamics, including turbulent burning rate, turbulent flame surface area, flame stretch factor, local displacement speed, flame dilatation, flame strain rate (both tangential and normal) and flame alignment with the principal strain rate are quantified, providing a full picture of near-wall turbulent flames for the considered conditions.