The nonlinear waves in a sheared liquid film on a horizontal plate at small Reynolds numbers are examined by theoretical and numerical approaches. The analysis employs the long-wave approximation along with finite difference schemes. The results show that the surface tension can suppress disturbances and prevent the occurrence of singularities. While the film flow is driven by the shear stress on the interface, its instability highly depends on the magnitude and direction of gravity. Specifically, when the direction of gravity is opposite to the wall-normal direction, perturbations are stabilized by gravity. In contrast, when these two directions are the same, the gravitational force is destabilizing, and stationary travelling waves can exist if a balance is reached between the effects of gravity and surface tension. For the steady solitary waves, there are quasi-periodic oscillations occurring between two stationary points, indicating the presence of heteroclinic trajectories. For periodic waves, the evolutions are sensitive to several parameters and initial disturbances, while one steady-state wave exhibits a sine function-like behaviour.

In transonic flow conditions, buffeting associated with finite-amplitude lift fluctuations can limit the operational envelope of an aircraft. For both airfoils and wings, these oscillations have been linked to global flow instabilities that arise from a Hopf bifurcation. We employ a combination of numerical simulations and global stability analysis to investigate the near-critical behaviour of the oscillatory buffet-onset instability on airfoils. The flow is governed by the unsteady Reynolds-averaged Navier–Stokes equations, with a basic state provided by a steady-state solution. In the weakly nonlinear formulation, the disturbance amplitude is described by the Landau equation. The linear growth rate can be determined from either the simulations or the stability analysis, and the Landau constant is derived from simulations resulting in finite-amplitude equilibrium states. The results show that the Landau constant is nearly independent of Mach number and angle of attack for a given airfoil. Using the Landau constant derived from a small number of simulations, the stability analysis can be employed to efficiently capture the essential finite-amplitude behaviour needed to estimate the buffet-onset boundary. The stability analysis is shown to capture the envelope of lift oscillations during a continuous pitch of an airfoil, from pre-buffet through post-buffet lift levels.

Understanding the generation of large-scale magnetic fields and flows in magnetohydro-dynamical (MHD) turbulence remains one of the most challenging problems in astrophysical fluid dynamics. Although much work has been done on the kinematic generation of large-scale magnetic fields by turbulence, relatively little attention has been paid to the much more difficult problem in which fields and flows interact on an equal footing. The aim is to find conditions for long-wavelength instabilities of stationary MHD states. Here, we first revisit the formal exposition of the long-wavelength linear instability theory, showing how long-wavelength perturbations are governed by four mean field tensors; we then show how these tensors may be calculated explicitly under the ‘short-sudden’ approximation for the turbulence. For MHD states with relatively little disorder, the linear theory works well: average quantities can be readily calculated, and stability to long-wavelength perturbations determined. However, for disordered basic states, linear perturbations can grow without bound and the purely linear theory, as formulated, cannot be applied. We then address the question of whether there is a linear response for sufficiently weak mean fields and flows in a dynamical (nonlinear) evolution, where perturbations are guaranteed to be bounded. As a preliminary study, we first address the nature of the response in a series of one-dimensional maps. For the full MHD problem, we show that in certain circumstances, there is a clear linear response; however, in others, mean quantities – and hence the nature of the response – can be difficult to calculate.

Self-sustained thermoacoustic oscillations as observed in low-emission combustion- involved gas turbines and aero-engines involve complicated thermal fluid–acoustics interaction and rich nonlinear dynamics. Such pulsating oscillations are known as thermoacoustic instability. When it occurs, large-amplitude limit cycle oscillations (LCOs) of thermodynamic parameters are frequently observed. These LCOs could cause overheating, flame flashback, and even engine failures. Thus it is critical to understand and predict the generation mechanisms and nonlinear dynamics behaviours, and then develop corresponding control approaches to prevent or control the onset of such instabilities. In this work, we develop and extend the classical van der Pol oscillators by integrating a physics-informed neural networks (PINNs) algorithm with a modelled nonlinear Rijke-type thermoacoustic combustor. The theoretical Rijke tube system (with Galerkin expansion and modified King's law implemented) and a CFD simulation model are applied to provide ‘training/calibration data’ for the extended van der Pol (EVDP)-PINNs model. The optimized EVDP oscillators are confirmed to be capable of capturing the key nonlinear characteristics by comparing the transient growth behaviours of thermodynamic perturbations and LCO amplitude and frequency. Further investigations are conducted to obtain Hopf bifurcation and amplitude death (AD) characteristics. Comparison is then made to the coupled EVDP systems. Quite similar Hopf bifurcation features, but differences in regions of AD, are observed. In general, we demonstrate an applicable approach to intelligently ‘learn’ a nonlinear thermoacoustic system and to create reliable EVDP oscillator systems, which have great potential to contribute to the development and testing of control approaches, such as the coupling described in this work, which may replace costly experimental tests.

The nonlinear evolution of free-stream vortical disturbances entrained in the entrance region of a circular pipe is investigated using asymptotic and numerical methods. Attention is focused on the low-frequency disturbances that induce streamwise elongated structures. A pair of vortical modes with opposite azimuthal wavenumbers is used to model the free-stream disturbances. Their amplitude is assumed to be intense enough for nonlinear interactions to occur inside the pipe. The formation and evolution of the perturbation flow are described by the nonlinear unsteady boundary-region equations in the cylindrical coordinate system, derived and solved herein for the first time. Matched asymptotic expansions are employed to construct appropriate initial conditions and the initial–boundary value problem is solved numerically by a marching procedure in the streamwise direction. Numerical results show the stabilising effect of nonlinearity on the intense algebraic growth of the disturbances and an increase of the wall-shear stress due to the nonlinear interactions. A parametric study is carried out to evince the effect of the Reynolds number, the streamwise and azimuthal wavelengths, and the radial length scale of the inlet disturbance on the nonlinear flow evolution. Elongated pipe-entrance nonlinear structures (EPENS) occupying the whole pipe cross-section are discovered. EPENS with $h$-fold rotational symmetry comprise $h$ high-speed streaks positioned near the wall, and $h$ low-speed streaks centred around the pipe core. These distinct structures display a striking resemblance to nonlinear travelling waves found numerically and observed experimentally in fully developed pipe flow. Good agreement of our mean-flow and root mean square data with experimental measurements is obtained.

The free-stream turbulence (FST) induced transition in perfect and non-ideal gas zero-pressure-gradient flat-plate boundary layers is investigated by means of large-eddy simulations. The study focuses on the influence of large incoming disturbances over the laminar-to-turbulent transition, by comparing two different integral length scales $L_f$, which differ by a factor of seven, at different FST intensities $T_u$. High-subsonic dense-gas boundary layers of the organic vapour Novec649, representative of organic Rankine cycle applications, are compared with air flows at Mach numbers 0.1 and 0.9. Compressibility and non-ideal gas effects are shown to be of minor importance in comparison to the influence of the FST integral length scale $L_f$. An increase of the inlet turbulent intensity always promotes transition, whereas an increase of $L_f$ has a double effect on the transition onset. At $T_u=2.5\,\%$, increasing $L_f$ promotes the transition, while it tends to delay transition for an FST intensity of 4 %. Larger FST integral scales tend to increase the spanwise distance between laminar streaks generated in the boundary layer. Two competing transition scenarios are observed. When the incoming turbulence intensity and length scale are moderate, the classical bypass route consists in the linear non-modal growth of streaks, which then experience secondary instabilities (sinuous or varicose) and lead to the generation of turbulent spots. The second scenario is characterized by the appearance of $\Lambda$-shaped structures near the inlet, which are further stretched to hairpin vortices before breaking down to turbulence. Spot inceptions can therefore occur at earlier locations than the streak growth. We are then faced with a competition between the classical bypass transition and nonlinear response mechanisms that ‘bypass’ this route. The present case at high $L_f$ and low $T_u$ is an example of a competing scenario, but even for the higher $T_u$ and $L_f$ conditions, only approximately one-third of turbulent spots are due to the $\Lambda$-shaped events. The nonlinear alternative route has strong similarities with scenarios described previously in the literature in the presence of leading edge effects or due to passing wakes. Such a path is governed by the turbulence intensity, but also by the integral length scale, with both parameters playing a critical role in the generation of the $\Lambda$-shaped structures near the inlet. This alternative mechanism is found to be robust under varying flow and thermodynamic conditions.

It is well known that buoyancy suppresses, and can even laminarise, turbulence in upward heated pipe flow. Heat transfer seriously deteriorates in this case. A new direct numerical simulation model is established to simulate flow-dependent heat transfer in an upward heated pipe. The model shows good agreement with experimental results. Three flow states are simulated for different values of the buoyancy parameter $C$: shear turbulence, laminarisation and convective turbulence. The latter two regimes correspond to the heat transfer deterioration regime and the heat transfer recovery regime, respectively (Jackson & Li 2002; Bae et al., Phys. Fluids, vol. 17, issue 10, 2005; Zhang et al., Appl. Energy, vol. 269, 2020, 114962). We confirm that convective turbulence is driven by a linear instability (Su & Chung, J. Fluid Mech., vol. 422, 2000, pp. 141–166) and that the deteriorated heat transfer within convective turbulence is related to a lack of rolls near the wall, which leads to weak mixing between the flow near the wall and the centre of the pipe. Having surveyed the fundamental properties of the system, we perform a nonlinear non-modal stability analysis, which seeks the minimal perturbation that triggers a transition from the laminar state. Given the differences between shear and convective turbulence, we aim to determine how the nonlinear optimal (NLOP) changes as the buoyancy parameter $C$ increases. We find that at first, the NLOP becomes thinner and closer to the wall. Most importantly, the critical initial energy $E_0$ required to trigger turbulence keeps increasing, implying that attempts to trigger it artificially may not be an efficient means to improve heat transfer at larger $C$. At $C=6$, a new type of NLOP is discovered, capable of triggering convective turbulence from lower energy, but over a longer time. It is active only in the centre of the pipe. We next compare the transition processes, from linear instability and by the nonlinear non-modal excitation. At $C=4$, linear instability leads to a state that approaches a travelling wave solution or periodic solutions, while the minimal seed triggers shear turbulence before decaying to convective turbulence. Deeper into the parameter space for convective turbulence, at $C=6$, the new nonlinear optimal triggers convective turbulence directly. Detailed analysis of the periodic solution at $C=4$ reveals three stages: growth of the unstable eigenfunction, the formation of streaks, and the decay of the streaks. The stages of the cycle correspond to changes in the linear instability of the turbulent mean velocity profile. Unlike the self-sustaining process for classical shear flows, where the streak is disrupted via instability, here, decay of the streak is more closely linked to suppression of the linear instability of the mean flow, and hence suppression of the rolls. Flow visualisations at $C$ up to $10$ also show similar processes, suggesting that the convective turbulence in the heat transfer recovery regime is sustained by these three typical processes.

Elastoinertial turbulence (EIT) is a chaotic flow resulting from the interplay between inertia and viscoelasticity in wall-bounded shear flows. Understanding EIT is important because it is thought to set a limit on the effectiveness of turbulent drag reduction in polymer solutions. Here, we analyse simulations of two-dimensional EIT in channel flow using spectral proper orthogonal decomposition (SPOD), discovering a family of travelling wave structures that capture the sheetlike stress fluctuations that characterise EIT. The frequency-dependence of the leading SPOD mode contains distinct peaks and the mode structures corresponding to these peaks exhibit well-defined travelling structures. The structure of the dominant travelling mode exhibits shift–reflect symmetry similar to the viscoelasticity-modified Tollmien–Schlichting (TS) wave, where the velocity fluctuation in the travelling mode is characterised by large-scale regular structures spanning the channel and the polymer stress field is characterised by thin, inclined sheets of high polymer stress localised at the critical layers near the channel walls. The travelling structures corresponding to the higher-frequency modes have a very similar structure, but are nested in a region roughly bounded by the critical layer positions of the next-lower-frequency mode. A simple theory based on the idea that the critical layers of mode $\kappa$ form the ‘walls’ for the structure of mode $\kappa +1$ yields quantitative agreement with the observed wave speeds and critical layer positions, indicating self-similarity between the structures. The physical idea behind this theory is that the sheetlike localised stress fluctuations in the critical layer prevent velocity fluctuations from penetrating them.

Recently, subcritical transition to turbulence in the quasi-two-dimensional (quasi-2-D) shear flow with strong linear friction (Camobreco et al., J. Fluid Mech., vol. 963, 2023, R2) has been demonstrated by the 2-D mechanism at $Re = 71\,211$, and the nonlinear Tollmien–Schlichting (TS) waves related to the edge state were approached independently of initial optimal disturbances. For 2-D plane Poiseuille flow, transition to the fully developed turbulence requires that the Reynolds number is several times larger than the critical Reynolds number $Re_c$ (Markeviciute & Kerswell, J. Fluid Mech., vol. 917, 2021, A57). In this paper, we observed the subcritical transitional flow in 2-D plane Poiseuille flow driven by the nonlinear TS waves by both linear and nonlinear optimal disturbances ($Re < Re_c$) with different quantitative edge states. The nonlinear optimal disturbances could trigger the sustained subcritical transitional flow for $Re \geqslant 2400$. The initial energy for nonlinear optimal disturbance is more efficient than the linear optimal disturbance in reaching the subcritical transitional flow for $2400 \leqslant Re \leqslant 5000$. Moreover, the initial energy of linear optimal disturbance is larger than the energy of its edge state. The nonlinear TS waves along the edge state are formed by the nonlinear optimal disturbances to trigger transitional flow, which agrees well with the main conclusions of Camobreco et al. (J. Fluid Mech., vol. 963, 2023, R2), while the required $Re$ of 2-D plane Poiseuille flow is much smaller.

The stimulation of instability and transport in the bottom boundary layer by internal solitary waves has been documented for over twenty years. However, the challenge of shallow slopes and a disparity of scales between the large-scale wave and the small-scale boundary layer has proven challenging for simulations. We present laboratory scale simulations that resolve the three-dimensionalisation in the boundary layer during the entire shoaling process. We find that the late stage, in which the incoming wave fissions into boluses, provides the most consistent source of three-dimensionalisation. In the early stage of shoaling, three-dimensionalisation occurs not so much due to separation bubble instability, but due to the interaction of vortices shed from the separation bubble with the overlying pycnocline. This interaction overturns the pycnocline, and creates bursts in kinetic energy and viscous dissipation, suggesting that the shed vortices induce turbulent motion and sediment resuspension in the water column above and behind the separation bubble.

Two-dimensional free-surface flow over localised topography is examined, with the emphasis on the stability of hydraulic-fall solutions. A Gaussian topography profile is assumed with a positive or negative amplitude modelling a bump or a dip, respectively. Steady hydraulic-fall solutions to the full incompressible, irrotational Euler equations are computed, and their linear and nonlinear stability is analysed by computing eigenspectra of the pertinent linearised operator and by solving an initial value problem. The computations are carried out numerically using a specially developed computational framework based on the finite-element method. The Hamiltonian structure of the problem is demonstrated, and stability is determined by computing eigenspectra of the pertinent linearised operator. It is found that a hydraulic-fall flow over a bump is spectrally stable. The corresponding flow over a dip is found to be linearly unstable. In the latter case, time-dependent simulations show that ultimately, the flow settles into a time-periodic motion that corresponds to an invariant solution in an appropriately defined phase space. Physically, the solution consists of a localised large-amplitude wave that pulsates above the dip while simultaneously emitting nonlinear cnoidal waves in the upstream direction and multi-harmonic linear waves in the downstream direction.

We present evidence revealing that an object with specific properties can exhibit multiple stable falling postures at low Reynolds numbers. By scrutinizing the force equilibrium relationship of a fixed object at various attack angles and Reynolds numbers, we introduce a methodology that can obtain the stable falling postures of the object. This method saves computational resources and more intuitively presents the results in the full parameter domain. Our findings are substantiated by free-fall tests conducted through both physical experiments and numerical simulations, which validate the existence of multiple stable solutions in accordance with the interpolation results obtained with fixed objects. Additionally, we quantify the abundance and distribution patterns of stable falling postures for a diverse range of representative shapes. This discovery highlights the existence of multiple stable solutions that are universally present across objects of different shapes. The implications of this research extend to the design, stability control and trajectory prediction of all free and controlled flights in both air and water.

Numerous studies have indicated that turbulence typically initiates along the boundary layer of the stationary disk within a rotor–stator cavity. To describe the transition process to turbulence on the stationary side of a closed rotor–stator cavity, a comprehensive approach combining global linear stability analysis with direct numerical simulation was adopted in the present study. The proposed model aligns with that of Yim et al. (J. Fluid Mech., vol. 848, 2018, pp. 631–647), who investigated the stability characteristics of the rotating-disk boundary layer in a rotor–stator cavity. In order to achieve a stable inflow for the stationary-disk boundary layer, we rotate the shroud together with the rotating disk. Through careful global stability analysis, the predominant spiral mode exhibiting the highest instability in the boundary layer of the stationary disk was discerned, corroborating observations from simulations. Initially, the spiral mode undergoes linear amplification, reaches a state of linear saturation and enters the nonlinear regime. Following nonlinear saturation in the flow field, a circular wave mode arises due to the influence of mean flow distortion. As the Reynolds number attained a sufficiently high level, the interplay between the downstream-propagating circular mode and spiral mode amplified disturbances in the boundary layer of the stationary disk, ultimately leading to the development of localised turbulence at the mid-radius of the rotor–stator cavity. Notably, the present study is the first to elucidate the coexistence of laminar–transitional–turbulent flow states in the stationary-disk boundary layer through direct numerical simulations.

Snap-through is a buckling instability that allows slender objects, including those in plant and biological systems, to generate rapid motion that would be impossible if they were to use their internal forces exclusively. In microfluidic devices, such as micromechanical switches and pumps, this phenomenon has practical applications for manipulating fluids at small scales. The onset of this elastic instability often drives the surrounding fluid into motion – a process known as snap-induced flow. To analyse the complex dynamics resulting from the interaction between a sheet and a fluid, we develop a prototypical model of a thin sheet that is compressed between the two sides of a closed channel filled with an inviscid fluid. At first, the sheet bends towards the upstream direction and the system is at rest. However, once the pressure difference in the channel exceeds a critical value, the sheet snaps to the opposite side and drives the fluid dynamics. We formulate an analytical model that combines the elasticity of thin sheets with the hydrodynamics of inviscid fluids to explore how external pressure differences, material properties and geometric factors influence the system's behaviour. To analyse the early stages of the evolution, we perform a linear stability analysis and obtain the growth rate and the critical pressure difference for the onset of the instability. A weakly nonlinear analysis suggests that the system can exhibit a pressure spike in the vicinity of the inverted configuration.

It is well known that inertia-free shearing flows of a viscoelastic fluid with curved streamlines, such as the torsional flow between a rotating cone and plate or the flow in a Taylor–Couette geometry, can become unstable to a three-dimensional time-dependent instability at conditions exceeding a critical Weissenberg ($Wi$) number. However, the combined effects of fluid elasticity, shear thinning and finite inertia (as quantified by the Reynolds number $Re$) on the onset of elasto-inertial instabilities are not fully understood. Using a set of cone–plate geometries, we experimentally explore the entire $Wi$–$Re$ phase space for a series of nonlinear viscoelastic fluids (with the dependence on shear rate $\dot{\gamma}$ quantified using a shear-thinning parameter $\beta _P(\dot {\gamma })$). We tune $\beta _P(\dot {\gamma })$ by varying the dissolved polymer concentration in solution. This progressively reduces shear thinning but leads to finite inertial effects before the onset of elastic instability, and thus naturally results in elasto-inertial coupling. Time-resolved rheometric measurements and flow visualization experiments allow us to investigate the effects of flow geometry, and document the combined effects of varying $Wi, Re$ and $\beta _P(\dot {\gamma })$ on the emergence of secondary motions at the onset of instability. The resulting critical state diagram quantitatively depicts the competition between the stabilizing effects of shear thinning and the destabilizing effects of inertia. We extend the curved streamline instability criterion of Pakdel & McKinley (Phys. Rev. Lett., vol. 77, no. 12, 1996, p. 2459) for the onset of purely elastic instability in curvilinear geometries by using scaling arguments to incorporate shear thinning and finite inertial effects. The augmented condition facilitates predictions of the onset of instability over a broader range of flow conditions, and thus bridges the gap between purely elastic and elasto-inertial curved streamline instabilities.

In the present paper, the sloshing flow in a cuboid tank forced to oscillate horizontally is investigated with both experimental and numerical approaches. The filling depth chosen is $h/L=0.35$ (with h the water depth and L the tank height), which is close to the critical depth. According to Tadjbakhsh & Keller (J. Fluid Mech., vol. 8, issue 3, 1960, pp. 442–451), as the depth passes through this critical value the response of the resonant sloshing dynamics changes from ‘hard spring’ to ‘soft spring’. The experimental tank has a thickness of $0.1L$, reducing three-dimensional effects. High-resolution digital camera and capacitance wave probes are used for time recording of the surface elevation. By varying the oscillation period and the amplitude of the motion imposed on the tank, different scenarios are identified in terms of free-surface evolution. Periodic and quasi-periodic regimes are found in most of the frequencies analysed but, among these, sub-harmonic regimes are also identified. Chaotic energetic regimes are found with motions of greater amplitude. Typical tools of dynamical systems, such as Fourier spectra and phase maps, are used for the regime identification, while the Hilbert–Huang transform is used for further insight into doubling-frequency and tripling-period bifurcations. For the numerical investigation, an advanced and well-established smoothed particle hydrodynamics method is used to aid the understanding of the physical phenomena involved and to extend the range of frequencies investigated experimentally.

The cubic interactions in a discrete system of four weakly nonlinear waves propagating in a conservative dispersive medium are studied. By reducing the problem to a single ordinary differential equation governing the motion of a classical particle in a quartic potential, the complete explicit branches of solutions are presented, either steady, periodic, breather or pump, thus recovering or generalizing some already published results in hydrodynamics, nonlinear optics and plasma physics, and presenting some new ones. Various stability criteria are also formulated for steady equilibria. Theory is applied to deep-water gravity waves for which models of isolated quartets are described, including bidirectional standing waves and quadri-directional travelling waves, steady or not, resonant or not.

This study extends our previous work (McCloughan & Suslov, J. Fluid Mech., vol. 887, 2020, A23), where the existence of a saddle-node bifurcation of steady axisymmetric electrolyte flows driven by the Lorentz force in a shallow annular domain was first reported. Here we perform further weakly nonlinear analysis over a wider range of the governing parameters to demonstrate that the previously reported saddle-node bifurcation is a local feature of a global fold catastrophe, which, in turn, is a part of cusp catastrophe occurring as the thickness of the fluid layer increases. The amplitude equation characterising multiple flow solutions in the finite vicinity of catastrophe points is derived. The sensitivity of its coefficients and solutions to the distance from the catastrophe points is assessed demonstrating the robustness of the used analytical procedure. The asymptotic flow solution past the catastrophe point is subsequently obtained and its topology is explored confirming the existence of the secondary circulation in the bulk of flow (two-tori background flow structure). The latter is argued to lead to the appearance of experimentally observable vortices on the fluid surface. The rigorous justification of this conjecture is to be given in Part 2 of the study.

An extended turbulent state can coexist with the stable laminar state in pipe flows. We focus here on short pipes with additional discrete symmetries imposed. In this case, the boundary between the coexisting basins of attraction, often called the edge of chaos, is the stable manifold of an edge state, which is a lower-branch travelling wave solution. We show that a low-dimensional submanifold of the edge of chaos can be constructed from velocity data using the recently developed theory of spectral submanifolds (SSMs). These manifolds are the unique smoothest nonlinear continuations of non-resonant spectral subspaces of the linearized system at stationary states. Using very low-dimensional SSM-based reduced-order models, we predict transitions to turbulence or laminarization for velocity fields near the edge of chaos.