We consider Euler flows on two-dimensional (2-D) periodic domain and are interested in the stability, both linear and nonlinear, of a simple equilibrium given by the 2-D Taylor–Green vortex. As the first main result, numerical evidence is provided for the fact that such flows possess unstable eigenvalues embedded in the band of the essential spectrum of the linearized operator. However, the unstable eigenfunction is discontinuous at the hyperbolic stagnation points of the base flow and its regularity is consistent with the prediction of Lin (Intl Math. Res. Not., vol. 2004, issue 41, 2004, pp. 2147–2178). This eigenfunction gives rise to an exponential transient growth with the rate given by the real part of the eigenvalue followed by passage to a nonlinear instability. As the second main result, we illustrate a fundamentally different, non-modal, growth mechanism involving a continuous family of uncorrelated functions, instead of an eigenfunction of the linearized operator. Constructed by solving a suitable partial differential equation (PDE) optimization problem, the resulting flows saturate the known estimates on the growth of the semigroup related to the essential spectrum of the linearized Euler operator as the numerical resolution is refined. These findings are contrasted with the results of earlier studies of a similar problem conducted in a slightly viscous setting where only the modal growth of instabilities was observed. This highlights the special stability properties of equilibria in inviscid flows.

Three-dimensional (3-D) clustering characteristics of large-Stokes-number sprays interacting with turbulent swirling co-flows are investigated experimentally. The astigmatic interferometric particle imaging (AIPI) technique is used for simultaneous measurement of the spray droplets position in 3-D space and their corresponding diameter. The Stokes number estimated based on the Kolmogorov time scale varies from 34 to 142. The results show that the degree of droplet clustering plateaus at about 0.4 and at large Stokes numbers. It is obtained that the mean length scale of the clusters normalized by the Kolmogorov length scale follows a power-law relation, and the mean void length scale normalized by the integral length scale plateaus at about 1.5 and at large Stokes numbers. It is shown that the ratio of the number density of the droplets residing within the clusters to the global number density increases with increasing Stokes number and is about 8 for the largest Stokes number examined in this study. The joint characteristics of cluster's normalized volume and the mean diameter of droplets residing within the clusters show that small-volume clusters accommodate droplets with a relatively broad range of diameters. However, large clusters carry droplets with the most probable diameter. The developed AIPI technique in the present study and the corresponding spray characteristics are of importance for engineering applications that aim to understand the 3-D clustering characteristics of large-Stokes-number droplets sprayed into turbulent swirling co-flows.

Liquid droplet dynamics are widely used in biological and engineering applications, which contain complex interfacial instabilities and pattern formation such as droplet merging, splitting and transport. This paper studies a class of mean field control formulations for these droplet dynamics, which can be used to control and manipulate droplets in applications. We first formulate the droplet dynamics as gradient flows of free energies in modified optimal transport metrics with nonlinear mobilities. We then design an optimal control problem for these gradient flows. As an example, a lubrication equation for a thin volatile liquid film laden with an active suspension is developed, with control achieved through its activity field. Lastly, we apply the primal–dual hybrid gradient algorithm with high-order finite-element methods to simulate the proposed mean field control problems. Numerical examples, including droplet formation, bead-up/spreading, transport, and merging/splitting on a two-dimensional spatial domain, demonstrate the effectiveness of the proposed mean field control mechanism.

The depth-integrated horizontal momentum equations and continuity equation are employed to develop a new model. The vertical velocity and pressure can be expressed exactly in terms of horizontal velocities and free-surface elevation, which are the only unknowns in the model. Dividing the water column into elements and approximating horizontal velocities using linear shape function in each element, a set of model equations for horizontal velocities at element nodes is derived by adopting the weighted residual method. These model equations can be applied for transient or steady free-surface flows by prescribing appropriate lateral boundary conditions and initial conditions. Here, only the wave–current–bathymetry interaction problems are investigated. Theoretical analyses are conducted to examine various linear wave properties of the new models, which outperform the Green–Naghdi-type models for the range of water depth to wavelength ratios and the Boussinesq-type models as they are capable of simulating vertically sheared currents. One-dimensional horizontal numerical models, using a finite-difference method, are applied to a wide range of wave–current–bathymetry problems. Numerical validations are performed for nonlinear Stokes wave and bichromatic wave group propagation in deep water, sideband instability, regular wave transformation over a submerged shoal and focusing wave group interacting with linearly sheared currents in deep water. Very good agreements are observed between numerical results and laboratory data. Lastly, numerical experiments of wave shoaling from deep to shallow water are conducted to further demonstrate the capability of the new model.

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.

The dynamics of a shock-induced separation unit generated by a 20$^\circ$ sharp fin placed on a cylindrical surface in a Mach 2.5 flow was investigated. Specifically, the present work investigated the mechanisms that govern the mid-frequency range of separation shock unsteadiness in the fin shock wave–boundary layer interaction (SBLI) unit. Two-dimensional pressure fields were obtained over the cylinder surface spanning the entire fin SBLI unit using high-bandwidth pressure-sensitive paint at 40 kHz imaging rate that allowed probing the low- through mid-frequency ranges of the separation shock unsteadiness. The mean pressure field showed a progressive weakening of the separation shock with downstream distance, which is an artifact of the three-dimensional relief offered by the curved mounting surface. The root-mean-square (r.m.s.) pressure field exhibited a banded structure with elevated $p_{r.m.s.}$ levels beneath the intermittent region, separation vortex and adjacent to the fin root. The power spectral density (PSD) of the surface pressure fluctuations obtained beneath the intermittent region revealed that the separation shock oscillations exhibited the mid-frequency content over the majority of its length. Interestingly, neither the PSD nor the length of the intermittent region varied noticeably with downstream distance, revealing a constant separation shock foot velocity along the entire SBLI. The pressure fluctuation PSD beneath the separation vortex also exhibited the broadband peak at the mid-frequency range of the separation shock motions over the majority of its length within the measurement domain. By contrast, the region adjacent to the fin root exhibited pressure oscillations at a substantially lower frequency compared with the separation shock and the separation vortex. Two-point coherence and cross-correlation analysis provided unique insights into the critical sources and mechanisms that drive the separation shock unsteadiness. The separation vortex and separation shock dynamics were found to be driven by a combination of convecting perturbations that originated from the vicinity of the fin leading edge and the local interactions of the separated flow with the incoming boundary layer. The boundary layer locally strengthened or weakened the convecting pressure perturbations depending on the local momentum fluctuations within the boundary layer.

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.

We study the sedimentation of U-shaped circular disks in the Stokes limit of vanishing inertia. We simulate the flow past such disks using a finite-element-based solution of the three-dimensional Stokes equations, accounting for the integrable singularities that develop along their edges. We show that the purely vertical sedimentation of such disks in their upright (upside-down) U orientation is unstable to perturbations about their pitching (rolling) axes. The instability is found to depend only weakly on the size of the container in which the disks sediment, allowing us to analyse their behaviour based on the resistance matrix which governs the evolution of the disk's six rigid-body degrees of freedom in an unbounded fluid. We show that the governing equations can be reduced to two ordinary differential equations which describe the disk's inclination against the direction of gravity. A phase-plane analysis, the results of which are in good agreement with experiments, reveals that the two instabilities generally cause the disk to sediment along complex spiral trajectories while it alternates between pitching- and rolling-dominated motions. The chirality of the trajectories is set by the initial conditions rather than the (non-chiral) shape of the disk. For certain initial orientations, the disk retains its inclination and sediments along a perfectly helical path. The observed behaviour is fundamentally different from that displayed by flat circular disks which sediment without any reorientation. We therefore study the effect of variations in the disk's curvature to show how in the limit of vanishing curvature the behaviour of a flat disk is recovered.

Zonal flows are mean flows in the east–west direction, which are ubiquitous on planets, and can be formed through ‘zonostrophic instability’: within turbulence or random waves, a weak large-scale zonal flow can grow exponentially to become prominent. In this paper, we study the statistical behaviour of the zonostrophic instability and the effect of magnetic fields. We use a stochastic white noise forcing to drive random waves, and study the growth of a mean flow in this random system. The dispersion relation for the growth rate of the expectation of the mean flow is derived, and properties of the instability are discussed. In the limits of weak and strong magnetic diffusivity, the dispersion relation reduces to manageable expressions, which provide clear insights into the effect of the magnetic field and scaling laws for the threshold of instability. The magnetic field mainly plays a stabilising role and thus impedes the formation of the zonal flow, but under certain conditions it can also have destabilising effects. Numerical simulation of the stochastic flow is performed to confirm the theory. Results indicate that the magnetic field can significantly increase the randomness of the zonal flow. It is found that the zonal flow of an individual realisation may behave very differently from the expectation. For weak magnetic diffusivity and moderate magnetic field strengths, this leads to considerable variation of the outcome, that is whether zonostrophic instability takes place or not in individual realisations.

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.

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.

In this paper, curved detonation equations with gradients for the pre-wave and post-wave are constructed followed by analysis, verification and applications. The study focuses on shock induced chemical reaction such as detonation, with the energy effect for the main attention. Equations consider both planar and transverse curvature to accommodate both planar and axisymmetric flow problems. Influence coefficients are derived and used to analyse the effect of energy and curvature on the post-wave gradient. Good agreement with the simulation results demonstrates that the equations presented in this paper can calculate various post-wave gradients accurately. After verification, the equations can be applied to applications, including not only solution and analysis but also in the inverse design. First, the method can be applied with polar analysis to provide a new perspective and higher order parameters for the study of detonation. Second, the equations can be used for the capture of detonation waves, where both planar and axisymmetric examples show better performance. Furthermore, the equations can be used in the inverse design of detonation waves in combination with the method of characteristics, which is one of the unique benefits of the present equations.

Three-dimensional vortex dynamics around two pitching foils arranged in side-by-side (parallel) configurations is numerically examined at a range of separation (gap) distances ($0.5c \leqslant y^* \leqslant 1.5c$). In-phase ($\phi =0$) and out-of-phase ($\phi ={\rm \pi}$) motions are considered for Strouhal numbers of $0.3$ and $0.5$ at a Reynolds number of $8000$. In this work, we show that the foil proximity effect, defined as the influence of one foil on the flow characteristics around the other, induces a spanwise instability in the braids of trailing-edge vortices (TEVs) during their roll-up. This is a newly identified instability that manifests itself in the form of secondary vortical structures with opposite circulation compared with the TEVs formed on the foils, which leads to the formation of double necking on the braids of the TEVs. We provide quantitative evidence linking the formation of these secondary structures to the braid instability. The first neck merges with the TEV, while the second neck detaches from the braid region and moves downstream independently. As the foil proximity effect intensifies (spacing between the foils decreases), secondary vortical structures, as well as the necks, become more prominent, leading to the emergence of three-dimensional wake features. Lastly, the influence of kinematics of the foils on three-dimensionality of the wake is investigated. At higher Strouhal numbers, broader regions of high strain are developed near the trailing edge, associated with the detachment of stronger structures from the braids of TEVs. The characterized instability demonstrates consistent properties for in-phase and out-of-phase motions, albeit with specific differences in dynamics of leading-edge vortices.

Effects of different filter kernels, namely, spectral cutoff ($\mathcal {S}$-filter) and Gaussian ($\mathcal {G}$-filter), on the geometrical properties of the subfilter stress (SFS) tensor and the filtered strain-rate (FSR) tensor are analysed in a forced homogeneous isotropic turbulence. Utilizing the Euler angle–axis methodology, it is observed that despite similar mean behaviour, the eigenframe alignment between SFS and FSR exhibits a non-trivially different statistical distribution for two different filters. Besides the eigenframe alignment, the eigenstructure of these tensors is also investigated. It is found that in contrast to the eigenstructure of the FSR which does not show sensitive dependence on the filter kernel type, the eigenstructure of the SFS tensor is significantly influenced by the filter type. Subsequently, the impact of different filter kernels on the subfilter energy flux (SFEF) is investigated. It is observed that energy transfer in $\mathcal {G}$-filtering is preferably distributed over the forward region, whereas for the $\mathcal {S}$-filter, the SFEF is more evenly distributed over both forward–backward regions, leading to a heavy energy transfer cancellation. Additionally, by decomposing the SFEF into different partial energy fluxes, it is found that the impact of the $\mathcal {S}$-filtering on the eigenstructure of the SFS leads to the amplification of the backward energy transfer. Conversely, the $\mathcal {G}$-filtering amplifies the forward energy transfer by producing a more pronounced alignment between the contractive–extensive eigenvectors.

Experimental studies on the sloshing of fluid layers are usually performed in rectangular tanks with fixed boundaries. In contrast, the present study uses a 4.76-m-long circular channel, a geometry with open periodic boundaries. Surface waves are excited by means of a submerged hill that, together with the tank, performs a harmonic oscillation. Laboratory measurements are made using 18 ultrasonic probes, evenly distributed over the channel to track the wave propagation. It is shown that a two-dimensional long-wave numerical model derived via the Kármán–Pohlhausen approach reproduces the experimental data as long as the forcing is monochromatic. The sloshing experiments imply a highly complex surface wave field. Different wave types such as solitary waves, undular bores and antisolitary waves are observed. For order one $\delta _{hill} = h_{hill}/h_0$, where $h_0$ is the mean water level and $h_{hill}$ the obstacle's height, the resonant reflections of solitary waves by the submerged obstacle give rise to an amplitude spectrum for which the main resonance peaks can be explained by linear theory. For smaller $\delta _{hill}$, wave transmissions lead to major differences with respect to the more common cases of sloshing with closed ducts having fully reflective ends for which wave transmission through the end walls is not possible. This ultimately results in more complex resonance diagrams and a pattern formation that changes rather abruptly with the frequency. The experiments are of interest not only for engineering applications but also for tidal flows over bottom topography.

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs (Phys. Fluids, vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$, where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.