The focus of this paper is a numerical simulation study of the flow dynamics in a periodic porous medium to analyse the physics of a symmetry-breaking phenomenon, which causes a deviation in the direction of the macroscale flow from that of the applied pressure gradient. The phenomenon is prominent in the range of porosity from 0.43 to 0.72 for circular solid obstacles. It is the result of the flow instabilities formed when the surface forces on the solid obstacles compete with the inertial force of the fluid flow in the turbulent regime. We report the origin and mechanism of the symmetry-breaking phenomenon in periodic porous media. Large-eddy simulation (LES) is used to simulate turbulent flow in a homogeneous porous medium consisting of a periodic, square lattice arrangement of cylindrical solid obstacles. Direct numerical simulation is used to simulate the transient stages during symmetry breakdown and also to validate the LES method. Quantitative and qualitative observations are made from the following approaches: (1) macroscale momentum budget and (2) two- and three-dimensional flow visualization. The phenomenon draws its roots from the amplification of a flow instability that emerges from the vortex shedding process. The symmetry-breaking phenomenon is a pitchfork bifurcation that can exhibit multiple modes depending on the local vortex shedding process. The phenomenon is observed to be sensitive to the porosity, solid obstacle shape and Reynolds number. It is a source of macroscale turbulence anisotropy in porous media for symmetric solid-obstacle geometries. In the macroscale, the principal axis of the Reynolds stress tensor is not aligned with any of the geometric axes of symmetry, nor with the direction of flow. Thus, symmetry breaking in porous media involves unresolved flow physics that should be taken into consideration while modelling flow inhomogeneity in the macroscale.

Two aspects of homogeneous rotating turbulence are quantified through forced direct numerical simulations in an elongated domain, which, in the direction of rotation, is approximately 340 times larger than the typical initial eddy size. First, by following the time evolution of the integral length scale along the axis of rotation $\ell _{\Vert }$, the growth rate of the columnar eddies and its dependence on the Rossby number $Ro_{\unicode[STIX]{x1D700}}$ is determined as $\unicode[STIX]{x1D6FE}=3.90\exp (-16.72\,Ro_{\unicode[STIX]{x1D700}})$ for $0.06\leqslant Ro_{\unicode[STIX]{x1D700}}\leqslant 0.31$, where $\unicode[STIX]{x1D6FE}$ is the non-dimensional growth rate. Second, a scaling law for the energy dissipation rate $\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D708}}$ is sought. Comparison with current available scaling laws shows that the relation proposed by Baqui & Davidson (Phys. Fluids, vol. 27(2), 2015, 025107), i.e. $\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D708}}\sim {u^{\prime }}^{3}/\ell _{\Vert }$, where $u^{\prime }$ is the root-mean-square velocity, approximates well part of our data, more specifically the range $0.39\leqslant Ro_{\unicode[STIX]{x1D700}}\leqslant 1.54$. However, relations proposed in the literature fail to model the data for the second and most interesting range, i.e. $0.06\leqslant Ro_{\unicode[STIX]{x1D700}}\leqslant 0.31$, which is marked by the formation of columnar eddies. To find a similarity relation for the latter, we exploit the concept of a spectral transfer time introduced by Kraichnan (Phys. Fluids, vol. 8(7), 1965, p. 1385). Within this framework, the energy dissipation rate is considered to depend on both the nonlinear time scale and the relaxation time scale. Thus, by analysing our data, expressions for these different time scales are obtained that result in $\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D708}}\sim (u^{\prime 4}Ro_{\unicode[STIX]{x1D700}}^{0.62}\unicode[STIX]{x1D70F}_{nl}^{iso})/\ell _{\bot }^{2}$, where $\ell _{\bot }$ is the integral length scale in the direction normal to the axis of rotation and $\unicode[STIX]{x1D70F}_{nl}^{iso}$ is the nonlinear time scale of the initial homogeneous isotropic field.

Turbulent channel flow simulation of a shear-thinning fluid is considered – see Arosemena et al. (J. Fluid Mech., vol. 908, 2021, p. A43) – and compared with a Newtonian base case to reveal the effects of the shear-dependent rheology on the near-wall structures. Analyses of different flow statistics revealed that, for the shear-thinning fluid case, the streamwise vortices appear to grow in size, depart from the wall and present a lessening in their intensity. Information regarding variations in the quasi-longitudinal vortices is also obtained from three-dimensional structures identified through a normalized $Q$-criterion. With shear-thinning rheology, it is shown that the structures are comprised of wall-attached and -detached families which are taller than for a Newtonian fluid. Also, for a given height, the structures appear to be longer, with approximately the same width and overall larger volume for the shear-thinning fluid case; albeit their fractal dimension remains the same when compared to the Newtonian base case. Moreover, it is observed that the number density of vortical structures decreases with shear-thinning fluid behaviour. These observations, in conjunction with the known changes to the longitudinal velocity structures which appear to be less streaky, more spanwise separated and thickened with shear-thinning rheology, strongly suggest that the near-wall self-sustaining process has been disrupted. As we move slightly away from the wall and with shear-thinning behaviour, the local increase in viscosity seems to lead to less energetic vortices whereas the streaks are provided with an additional source of energy due to fluctuations in viscosity.

In the Prandtl model for anabatic slope flows, a uniform positive buoyancy flux at the surface drives an upslope flow against a stable background stratification. In the present study, we conduct linear stability analysis of the anabatic slope flow under this model and contrast it against the katabatic case as presented in Xiao & Senocak (J. Fluid Mech., vol. 865, 2019, R2). We show that the buoyancy component normal to the sloped surface is responsible for the emergence of stationary longitudinal rolls, whereas a generalised Kelvin–Helmholtz (KH) type of mechanism consisting of shear instability modulated by buoyancy results in a streamwise-travelling mode. In the anabatic case, for slope angles larger than $9^{\circ }$ to the horizontal, the travelling KH mode is dominant whereas, at lower inclination angles, the formation of the stationary vortex instability is favoured. The same dynamics holds qualitatively for the katabatic case, but the mode transition appears at slope angles of approximately $62^{\circ }$. For a fixed slope angle and Prandtl number, we demonstrate through asymptotic analysis of linear growth rates that it is possible to devise a classification scheme that demarcates the stability of Prandtl slope flows into distinct regimes based on the dimensionless stratification perturbation number. We verify the existence of the instability modes with the help of direct numerical simulations, and observe close agreements between simulation results and predictions of linear analysis. For slope angle values in the vicinity of the junction point in the instability map, both longitudinal rolls and travelling waves coexist simultaneously and form complex flow structures.

In this paper, we study the role of shear-induced migration and particle-induced normal stresses in the formation and stability of a particle-laden, gravity-driven shallow flow. We first examine the modification of the base-state Nusselt flow due to the underlying microstructure, how shear-induced migration leads to viscosity stratification. We inspect the development of the base state via the boundary layer formation in the ‘shallow’ limit and find a reduction in entrance length with increasing bulk particle concentration and an increase in entrance length with increasing Péclet number ($Pe_p = \dot {\gamma } a^2 / D_0$, where $\dot{\gamma}$ is the average shear rate, a is the particle size and $D_0$ is the single particle diffusivity). A linear stability analysis is then performed on the fully developed state to identify two modes of instability typically found in gravity-driven falling films – the long-wave surface and the short-wave shear modes. We find that when the associated Péclet number is $Pe_p \ll 1$, increasing bulk particle volume fraction delays the onset of instability for both the surface mode and shear mode. However, with $Pe_p = {O}(1)$, we find an enhancement in both modes of instability. We also find that, beyond a critical Péclet number, for a fixed particle volume fraction, the surface mode is unstable even in the absence of fluid inertia. The enhanced destabilisation is attributed to the combined effects of base-state viscosity stratification and momentum forcing via particle concentration perturbations. We also show that the physics behind the enhancement of instability is independent of the choice of the constitutive model used to describe the dynamics of the particle phase, provided the chosen model has elements of shear-induced migration.

Flow past a rotating step cylinder is investigated through three-dimensional numerical simulations for a diameter ratio of 0.5 and a Reynolds number of 150. The step cylinder comprises two cylinders with different diameters arranged coaxially with a step between them. The rotation rate α is defined as the ratio of the rotation speed of the larger cylinder surface to the free-stream velocity. Vortex shedding happens for both cylinders at α = 0, 0.5 and 1, and is suppressed only for the larger cylinder at α = 2 and 3 and fully suppressed for both cylinders at α = 4. The vortex shedding suppression for the larger cylinder or for both cylinders has significant effects on the wake. The S-, N- and L-cells at α = 0 are in good agreement with those reported in previous studies and still exist at α = 1. The N-cell disappears at α = 0.5, and as a result, the L- and S-cells interact with each other directly at the step position. An additional cellular zone is found at α = 0.5 and 1 and this zone has multiple cells with vortex dislocation between them. At α = 2 and 3, there is a strong hub vortex in the streamwise direction behind the step and the vortices in the wake of the smaller cylinder form helical vortices after they roll around this hub vortex. The hub vortex is generated by the difference between the Magnus effect between the smaller and larger cylinders. At α = 4, the hub vortex still exists but the helical vortices disappear because vortex shedding is suppressed for both cylinders. At this rotation rate, the pressure on the cylinder surface oscillates with a frequency much higher than the vortex shedding frequency. The oscillation of the pressure is caused by the combination of periodic generation of ring vortices and their motion along the span of the larger cylinder.

To explore the impact of surface viscosity on coexisting fluid domains in biomembranes we consider two-phase fluid deformable surfaces as model systems for biomembranes. Such surfaces are modelled by incompressible surface Navier–Stokes–Cahn–Hilliard-like equations with bending forces. We derive this model using the Lagrange–d’Alembert principle considering various dissipation mechanisms. The highly nonlinear model is solved numerically to explore the tight interplay between surface evolution, surface phase composition, surface curvature and surface hydrodynamics. It is demonstrated that hydrodynamics can enhance bulging and furrow formation, which both can further develop to pinch-offs. The numerical approach builds on a Taylor–Hood element for the surface Navier–Stokes part, a semi-implicit approach for the Cahn–Hilliard part, higher-order surface parametrizations, appropriate approximations of the geometric quantities, and mesh redistribution. We demonstrate convergence properties that are known to be optimal for simplified subproblems.

The interaction of a hollow droplet impacting a solid surface occurs in several applications, including controllable biomedicine and thermal spray coating. Understanding the physics of the hollow droplet spreading is the key to maintaining the mass transfer process in all relevant applications. In this work, a comprehensive experimental, numerical and theoretical study is performed on water hollow droplets impacting a rigid surface to better understand the flattening process of a hollow droplet. In the numerical part, compressible Navier–Stokes equations are solved using the volume of fluid (VOF) method in a two-dimensional (2-D)-axisymmetric model. The comparison of simulation results with the experimental photographs shows that the numerical solution can correctly predict the hollow droplet shape evolution. The results show that the spreading diameter and height of the counter-jet formed after the hollow droplet impact grow with impact velocity. Investigating the size and location of the entrapped bubble shows an optimum bubble size that facilitates the hollow droplet flattening. It is also shown that the ripples on splats produced by the hollow droplets with a larger bubble size are higher than those of small bubbles. In the end, a theoretical model is developed to analyse the maximum spreading diameter of the hollow droplet impact analytically. Its prediction is in good agreement with the experimental and numerical results.

The motion of a slightly rarefied gas in a long straight two-dimensional channel caused by a discontinuous surface temperature is investigated on the basis of kinetic theory with a special interest in the fluid-dynamic description. More precisely, the channel is longitudinally divided into two parts and each part is kept at a uniform temperature different from each other, so that the surface temperature of the whole channel has a jump discontinuity at the junction. Under the assumption that the amount of jump in the surface temperature is small, the steady behaviour of the gas induced in the channel is studied on the basis of the linearized Boltzmann equation and the diffuse reflection boundary condition in the case where the Knudsen number, defined by the ratio of the molecular mean free path and the width of the channel, is small. Using a matched asymptotic expansion method combined with Sone’s asymptotics, a Stokes system describing the overall macroscopic behaviour of the gas inside the channel is derived, with a new feature of the ‘slip boundary condition’ for the flow velocity due to the jump discontinuity in the surface temperature of the channel. This condition takes the form of a diverging singularity with source and sink located at the point of discontinuity, with a multiplicative factor determined through the analysis of a spatially two-dimensional Knudsen-layer (or a Knudsen-zone) problem. Some numerical demonstrations based on the Bhatnagar–Gross–Krook equation are also presented.

Identifying the governing parameters of self-sustained oscillation is crucial for the diagnosis, prediction and control of thermoacoustic instabilities. In this paper, we propose and validate a novel method for computing the parameters of thermoacoustic oscillation in a stochastic environment, which exploits a physics-informed neural network (PINN). Specifically, we introduce a negative log-likelihood loss function that integrates the stochastic samples and the solution of the Fokker–Planck equation. The proposed framework is validated using the numerically generated signal and the experimental data obtained from an annular combustor, both before and after the supercritical Hopf bifurcation. The results of PINN-based system identification show good agreement with the actual system parameters and the original stochastic signal, with improved accuracy compared to established methods. To the best of our knowledge, this study constitutes the first demonstration of the PINN-inverse approach that uses the noise-induced dynamics of thermoacoustic systems, opening up new pathways for diagnosing and predicting the thermoacoustic behaviour of various combustion systems.

We report direct numerical simulations of the flow around a spanwise-flexible wing in forward flight. The simulations were performed at $Re=1000$ for wings of aspect ratio 2 and 4 undergoing a heaving and pitching motion at Strouhal number $St_c\approx 0.5$. We have varied the effective stiffness of the wing $\varPi _1$ while keeping the effective inertia constant, $\varPi _0=0.1$. It has been found that there is an optimal aerodynamic performance of the wing linked to a damped resonance phenomenon, that occurs when the imposed frequency of oscillation approaches the first natural frequency of the structure in the fluid, $\omega _{n,f}/\omega \approx 1$. In that situation, the time-averaged thrust is maximum, increasing by factor 2 with respect to the rigid case with an increase in propulsive efficiency of approximately 15 %. This enhanced aerodynamic performance results from the combination of larger effective angles of attack of the outboard wing sections and a delayed development of the leading edge vortex. With increasing flexibility beyond the resonant frequency, the aerodynamic performance drops significantly, in terms of both thrust production and propulsive efficiency. The cause of this drop lies in the increasing phase lag between the deflection of the wing and the heaving/pitching motion, which results in weaker leading edge vortices, negative effective angles of attack in the outboard sections of the wing, and drag generation in the first half of the stroke. Our results also show that flexible wings with the same $\omega _{n,f}/\omega$ but different aspect ratio have the same aerodynamic performance, emphasizing the importance of the structural properties of the wing for its aerodynamic performance.