The transition route and bifurcations of the buoyant flows developing on a heated horizontal circular surface are elaborated using direct numerical simulations and direct stability analysis. A series of bifurcations, as a function of Rayleigh numbers ($Ra$) ranging from $10^6$ to $6.0\times 10^7$, are found on the route to chaos of the flows at $Pr=7$. When $Ra<1.0\times 10^3$, the buoyant flows above the heated horizontal surface are dominated by conduction, because of which the distinct thermal boundary layer and plume are not present. At $Ra=1.1\times 10^6$, a Hopf bifurcation occurs, resulting in the flow transition from a steady state to a periodic puffing state. As $Ra$ increases further, the flow enters a periodic rotating state at $Ra=1.9\times 10^6$, which is a unique state that was rarely discussed in the literature. These critical transitions, leaving from a steady state and subsequently entering a series of periodic states (puffing, rotating, flapping and period-doubling) and finally leading to chaos, are diagnosed using two-dimensional Fourier transforms. Moreover, direct stability analysis is conducted by introducing random numerical perturbations into the boundary condition of the surface heating. We find that when the state of a flow is in the vicinity of critical values (e.g. $Ra=2.0\times 10^6$), the flow is conditionally unstable to perturbations, and it can bifurcate from the rotating state to the flapping state in advance. However, for relatively stable flow states, such as at $Ra=1.5\times 10^6$, the flow remains in its periodic puffing state even though it is being perturbed.

We explore predictions of two models of one-dimensional capillary rise in rigid and partially saturated porous media. One is an existing one from the literature and the second is a free-boundary model based on Richards’ equation with two moving boundaries of the evolving partially saturated region. Both models involve the specification of saturation-dependent functions for local capillary pressure and permeability and connect to classical models for saturated porous media. Existing capillary-rise experiments show two notable regimes: (i) an early-time regime typically well-described by classical capillary-rise theory in a fully saturated porous media, and (ii) a long-time regime that has anomalous dynamics in which the capillary-rise height may scale with a non-classical power law in time or have more complicated dynamics. We demonstrate that the predictions of both models compare well with experimental capillary-rise data over early- and long-time regimes gathered from three independent studies in the literature. The model predictions also shed light on recent scaling laws that relate the capillary pressure and permeability of the partially saturated media to the capillary-rise height. We use these models to probe computationally observed permeability relationships to capillary-rise height. We demonstrate that a recently proposed permeability scaling for the anomalous capillary-rise regime is indeed realized and is particularly apparent in the lower portion of the partially saturated media. For our free-boundary model we also compute capillary pressure measures and show that these reveal the linear relation between the capillary pressure and capillary-rise height expected for a capillarity–gravity balance in the upper portion of the partially saturated porous media.

Utilizing the discrete element method and the pore network model, we numerically investigate the impact of compaction on the longitudinal dispersion coefficient of porous media. Notably, the dispersion coefficient exhibits a non-monotonic dependence on the degree of compaction, which is distinguished by the presence of three distinct regimes in the variation of dispersion coefficient. The non-monotonic variation of dispersion coefficient is attributed to the disparate effect of compaction on dispersion mechanisms. Specifically, the porous medium tightens with an increasing pressure load, reducing the effect of molecular diffusion that primarily governs at small Péclet numbers. On the other hand, heightened pressure loads enhance the heterogeneity of pore structures, resulting in increased disorder and a higher proportion of stagnant zones within porous media flow. These enhancements further strengthen mechanical dispersion and hold-up dispersion, respectively, both acting at higher Péclet numbers. It is crucial to highlight that hold-up dispersion is induced by the low-velocity regions in porous media flow, which differ fundamentally from zero-velocity regions (such as dead-ends or the interior of permeable grains) as described by the classical theory of dispersion. The competition between weakened molecular diffusion and enhanced hold-up dispersion and mechanical dispersion, together with the shift in the dominance of dispersion mechanisms across various Péclet numbers, results in multiple regimes in the variation of dispersion coefficients. Our study provides unique insights into structural design and modulation of the dispersion coefficient of porous materials.

We use $123$ three-dimensional direct numerical simulations to study fingering convection in non-rotating spherical shells. We investigate the scaling behaviour of the flow length scale, the non-dimensional heat and compositional fluxes $Nu$ and $Sh$ and the mean convective velocity over the fingering convection instability domain defined by $1 \leq R_\rho < Le$, $R_\rho$ being the ratio of density perturbations of thermal and compositional origins and $Le$ the Lewis number. We show that the chemical boundary layers are marginally unstable and adhere to the laminar Prandtl–Blasius model, hence explaining the asymmetry between the inner and outer spherical shell boundary layers. We develop scaling laws for two asymptotic regimes close to the two edges of the instability domain, namely $R_\rho \lesssim Le$ and $R_\rho \gtrsim 1$. For the former, we develop novel power laws of a small parameter $\epsilon$ measuring the distance to onset, which differ from theoretical laws published to date in Cartesian geometry. For the latter, we find that the Sherwood number $Sh$ gradually approaches a scaling $Sh\sim Ra_\xi ^{1/3}$ when $Ra_\xi \gg 1$; and that the Péclet number accordingly follows $Pe \sim Ra_\xi ^{2/3} |Ra_T|^{-1/4}$, $Ra_T$ and $Ra_{\xi}$ being the thermal and chemical Rayleigh numbers. When the Reynolds number exceeds a few tens, we report on a secondary instability which takes the form of large-scale toroidal jets which span the entire spherical domain. Jets distort the fingers, resulting in Reynolds stress correlations, which in turn feed the jet growth until saturation. This nonlinear phenomenon can yield relaxation oscillation cycles.

We investigate two-dimensional (2-D) axisymmetric flow in toroidal geometry, with a focus on a transition between 2-D three-component flow and 2-D two-component flow. This latter flow state allows a self-organization of the system to a quiescent dynamics, characterized by long-living coherent structures. When these large-scale structures orient in the azimuthal direction, the radial transport is reduced. Such a transition, if it can be triggered in toroidally confined fusion plasmas, is beneficial for the generation of zonal flows and should consequently result in a flow field beneficial for confinement.

We present experiments on oscillating hydrofoils undergoing combined heaving and pitching motions, paying particular attention to connections between propulsive efficiency and coherent wake features extracted using modal analysis. Time-averaged forces and particle image velocimetry measurements of the flow field downstream of the foil are presented for a Reynolds number of $Re=11\times 10^3$ and Strouhal numbers in the range $St=0.16\unicode{x2013}0.35$. These conditions produce 2S and 2P wake patterns, as well as a near-momentumless wake structure. A triple decomposition using the optimized dynamic mode decomposition method is employed to identify dominant modal components (or coherent structures) in the wake. These structures can be connected to wake instabilities predicted using spatial stability analyses. Examining the modal components of the wake provides insightful explanations into the transition from drag to thrust production, and conditions that lead to peak propulsive efficiency. In particular, we find modes that correspond to the primary vortex development in the wakes. Other modal components capture elements of bluff body shedding at Strouhal numbers below the optimum for peak propulsive efficiency and characteristics of separation for Strouhal numbers higher than the optimum.

Cowl-induced incident shock wave/boundary layer interactions (ISWBLIs) under the influence of shoulder expansion represent one of the dominant phenomena in supersonic inlets. To provide a more comprehensive understanding of how an expansion corner affects the ISWBLI, a detailed experimental and analytical study is performed in a Mach 2.73 flow in this work. Pressure measurement, schlieren photography and surface oil-flow visualisation are used to record flow features, including the pressure distribution, separation extent and surface-flow topological structures. Our results reveal three types of ISWBLIs influenced by the expansion corner. When the shock intensity is weak, the separation is small scale with the expansion waves emanating from the expansion corner. This is the first type of expansion-corner-affected ISWBLI (EC-ISWBLI). When the incident shock wave is strong, large-scale separation occurs, accompanied by the disappearance of expansion waves, forming the second type of EC-SWBLI. The expansion corner induces a ‘lock-in’ effect in which the separation onset is consistently locked near the expansion corner regardless of the incident shock intensity and impingement position. The third type of EC-ISWBLI occurs when the shock is sufficiently strong and the impingement point is close to the expansion corner. In this interaction, the ‘lock-in’ effect ceases to manifest. Moreover, a shock polar-incorporating inviscid model is employed to elucidate the shock patterns. Two criteria are established by combining free interaction theory with this model. The first criterion provides valuable insights into the evolution of separations with a minimal overall pressure rise and the second criterion determines the threshold for the occurrence of the ‘lock-in’ effect.

The turbulent boundary layer (TBL) is a widely existing flow phenomenon in nature and engineering applications. Its strong mixing effect can achieve more sufficient material mixing, heat transport, etc. The understanding of the entrainment process and mechanism of irrotational fluids entering the turbulent region can be promoted by studying the geometric and dynamic characteristics of turbulent${/}$non-turbulent interfaces (TNTI). In compressible flow, it is unclear whether the properties of TNTI will change and whether the entrainment will show different features due to the influence of compressibility. Based on the direct numerical simulation results of supersonic compressible plate TBLs with Mach number of 2.9, the geometric and dynamic characteristics of TNTI are investigated in this paper. The interface is identified by the enstrophy method, and the height, thickness, fractal dimension, enstrophy transportation and entrainment characteristics of the interface are investigated. It is found that for the enstrophy transportation in a TBL, the contribution of compressibility-related terms accounts for approximately 13.4 % of the total enstrophy transportation, which tends to transfer the enstrophy of turbulence near the interface to both directions vertical to the interface. This promotes the expansion of the turbulent region towards the non-turbulent region, and the mean height, thickness and entrainment velocity are increased by approximately 3.7 %, 7.0 % and 8.5 %, respectively, while the fractal dimension is basically unaffected. Different from the incompressible flow, the contribution of the compressibility-related terms to the entrainment velocity is independent of the local curvature, and the intense entrainment process is more likely to occur on a highly curved concave surface.

In the present study, we performed direct numerical simulations for a hypersonic turbulent boundary layer over the windward side of a lifting body, the HyTRV model, at Mach number $6$ and attack angle 2$^{\circ }$ to investigate the global and local turbulent features, and evaluate its difference from canonical turbulent boundary layers. By scrutinizing the instantaneous and averaged flow fields, we found that the transverse curvature on the windward side of the HyTRV model induces the transverse opposing pressure gradients that push the flow on both sides towards the windward symmetry plane, yielding significant effects of the azimuthal inhomogeneity and large-scale cross-stream circulations, moderate and azimuthal independent influences of adverse pressure gradient, and negligible impact of the mean flow three-dimensionality. Further inspecting the local turbulent statistics, we identified that the mean and fluctuating velocity become increasingly similar to the highly decelerated turbulent boundary layers over flat plates in that the mean velocity deficit is enhanced, and the outer layer Reynolds stresses are amplified as it approaches the windward symmetry plane, and prove to be self-similar under the scaling of Wei & Knopp (J. Fluid Mech., vol. 958, 2023, A9) for adverse-pressure-gradient turbulent boundary layers. Conditionally averaged Reynolds stresses based on strong sweeping and ejection events demonstrated that the Kelvin–Helmholtz instability of the strong embedded shear layer induced by the large-scale cross-stream circulations is responsible for the turbulence amplification in the outer layer. The strong Reynolds analogy that relates the mean velocity and temperature was refined to incorporate the non-canonical effects, showing considerable improvements in the accuracy of such a formula. On the other hand, the temperature fluctuations are still transported passively, as indicated by their resemblance to the velocity. The conclusions obtained in the present study provide potentially profitable information for turbulent modelling modification for the accurate predictions of skin friction and wall heat transfer.