We use Bayesian inference, accelerated by adjoint methods, to construct a quantitatively accurate model of the thermoacoustic behaviour of a conical flame in a duct. We first perform a series of automated experiments on a ducted flame rig. Next, we propose several candidate models of the rig's components and assimilate data into each model to find the most probable parameters for that model. We rank the candidate models based on their marginal likelihood (evidence) and select the most likely model for each component. We begin this process by rigorously characterizing the acoustics of the cold rig. When the flame is introduced, we propose several candidate models for the fluctuating heat release rate. We find that the most likely flame model considers velocity perturbations in both the burner feed tube and the outer duct, even though studies in the literature typically neglect either one of these. Using the most likely model, we infer the flame transfer functions for 24 flames and quantify their uncertainties. We do this with the flames in situ, using only pressure measurements. We find that the inferred flame transfer functions render the model quantitatively accurate, and, where comparable, broadly consistent with direct measurements from several studies in the literature.

This work examines the influence of body porosity on the wake past nominally two-dimensional rectangular plates of fixed width $D$ in the moderate range of Reynolds numbers $Re = UD/\nu$ (with $U$ the incoming velocity and $\nu$ the kinematic viscosity) between 15 000 and 70 000. With porosity $\beta$ defined as the ratio of open to total area of the plate, it is well established that as porosity increases, the wake shifts from the periodic von Kármán shedding behaviour to a regime where this vortex shedding is absent. This change impacts the fluid forces acting on the plate, especially the drag, which is significantly lower for a wake without vortex shedding. We analyse experimentally the transition between these two regimes using hot-wire anemometry, particle-image velocimetry and force measurements. Coherence and phase measurements are used to determine the existence of regular, periodic vortex shedding based on the velocity fluctuations in the two main shear layers on either side of the wake. Results show that, independent of $Re$, the wake exhibits the classical Kármán vortex shedding pattern for $\beta <0.2$ but this is absent for $\beta >0.3$. In the intermediate range, $0.2<\beta <0.3$, there is a transitional regime that has not previously been identified; it is characterised by intermittent shedding. The flow alternates randomly between a vortex shedding and a non-shedding pattern and the total proportion of time during which vortex shedding is observed (the intermittency) decreases with increasing porosity.

The probability density function (PDF) for the free surface elevation in an irregular sea has an integral formulation when based on the cumulant generating function. To leading order, the result is Gaussian, whereas nonlinear extensions have long been limited to Gram–Charlier series approximations. As shown recently by Fuhrman et al. (J. Fluid Mech., vol. 970, 2023, A38), however, the second-order integral can be represented exactly in closed form. The present work extends this further, enabling determination of this PDF to even higher orders. Towards this end, a new ordinary differential equation (ODE) governing the PDF is first derived. Asymptotic solutions in the limit of large surface elevation are then found, utilizing the method of dominant balance. These provide new analytical forms for the positive tail of the PDF beyond second order. These likewise clarify how high-order cumulants (involving statistical moments such as the kurtosis) govern the tail, which is shown to get heavier with each successive order. The asymptotic solutions are finally utilized to generate boundary conditions, such that the governing ODE may be solved numerically, enabling novel determination of the PDF at third and higher order. Successful comparisons with challenging data sets confirm accuracy. The methodology thus enables the PDF of the surface elevation to be determined numerically, and the asymptotic tail analytically, to any desired order. Results are worked out explicitly up to fifth order. The theoretical probability of extreme surface elevations (typical of rogue waves) may thus be assessed quantitatively for highly nonlinear irregular seas, requiring only relevant statistical quantities as input.

We derive the spin Euler equation for ideal flows by applying the spherical Clebsch mapping. This equation is based on the spin vector, a unit vector field encoding vortex lines, instead of the velocity. The spin Euler equation enables a feasible Lagrangian study of fluid dynamics, as the isosurface of a spin-vector component is a vortex surface and material surface in ideal flows. We establish a non-blowup criterion for the spin Euler equation, suggesting that the Laplacian of the spin vector must diverge if the solution forms a singularity at some finite time. The direct numerical simulations (DNS) of three ideal flows – the vortex knot, the vortex link and the modified Taylor–Green flow – are conducted by solving the spin Euler equation. The evolution of the Lagrangian vortex surface illustrates that the regions with large vorticity are rapidly stretched into spiral sheets. The DNS result exhibits a pronounced double-exponential growth of the maximum norm of Laplacian of the spin vector, showing no evidence of the finite-time singularity formation if the double-exponential growth holds at later times. Moreover, the present criterion with Lagrangian nature appears to be more sensitive than the Beale–Kato–Majda criterion in detecting the flows that are incapable of producing finite-time singularities.

Air-blast atomizers are extensively used for a variety of purposes. Due to its complexity, the atomization mechanism has not been elucidated. In this study, a mechanistic model is proposed to predict the droplet diameter distribution based on the atomization process of a planar liquid film with co-current gas flows, and its validity is examined by comparing the estimated and measured droplet diameters using high-speed image analysis and laser measurement. As a result, using high-speed imaging, we clarified that the bag film rupture is caused not by the turbulence of the gas flow but by the impact of floating droplets on the liquid film of the expanding bag when the film is thin enough. The average thickness of the liquid film at the bag breakup is of the order of micrometres and varies greatly, resulting in a dispersed distribution of droplet diameters. After the film ruptures, the bag film shrinks towards its transversal and vertical rims due to surface tension, forming large-diameter ligaments. During the contraction process of the bag film, tiny droplets of the order of micrometers are formed at the edge of the perforation. Finally, the remaining ligaments with large diameters fragment into large droplets with submillimetre diameters. The good agreement between the measured and predicted droplet diameter distributions validated the mechanistic model.

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.

First-order sensitivities and adjoint analysis are used widely to control the linear stability of unstable flows. Second-order sensitivities have recently helped to increase accuracy. In this paper, a method is presented to calculate arbitrary high-order sensitivities based on Taylor expansions of the incompressible base flow and its eigenproblem around a scalar parameter. For the incompressible Navier–Stokes equations, general expressions for the sensitivities are derived, into which parameter-specific information can be inserted. The computational costs are low since, for all orders, a linear equation system has to be solved, of which the left-hand-side matrix stays constant and thus its preconditioning can be exploited. Two flow scenarios are examined. First, the cylinder flow equations are expanded around the inverse of the Reynolds number, enabling the prediction of the two-dimensional cylinder base flow and its leading eigenvalue as a function of the Reynolds number. This approach computes accurately the base flow and eigenvalue even in the unstable regime, providing, when executed subsequently, a mean to calculate unstable base flows. This case gives a clear introduction into the method and allows us to discuss its constraints regarding convergence behaviour. Second, a small control cylinder is introduced into the domain of the cylinder flow for stabilization. Higher-order sensitivity maps are calculated by modelling the small cylinder with a steady forcing. These maps help to identify stabilizing areas of the flow field for Reynolds numbers within the laminar vortex shedding regime, with the required number of orders increasing as the Reynolds number rises. The results obtained through the proposed method align well with numerically calculated eigenvalues that incorporate the cylinder directly into the grid.

This study experimentally investigates the wake structure of a porous square cylinder in terms of permeability over two decades of $Da$ (i.e. $2.4 \times 10^{-5} < Da < 2.9 \times 10^{-3}$). The porous cylinder, featuring a simple cubic lattice structure, was fabricated using an additive manufacturing technique. This unique method, combined with a periodic and scalable lattice structure, effectively isolates permeability from porosity, making it suitable for an in-depth parametric study. The key parameter, permeability, was directly estimated by measuring the pressure drop and superficial velocity for each porous case in an open-loop pipe flow system. The downstream flow fields were obtained using standard planar particle image velocimetry measurements in an open-loop wind tunnel. Based on the experimental data, structural modifications in the near wake were examined in relation to permeability, leading to the identification of four distinct flow regimes depending on $Da$. Additionally, the downstream flow adjustment length ($L_i$) was assessed by introducing a permeability-based source term into the momentum equation, facilitating the development of an analytical model for $L_i$. The present experimental data support this analytical model, and our results further confirmed that $L_i$ plays a crucial role as a characteristic length scale in the near wake.

Combustion instabilities in annular systems raise fundamental issues that are also of practical importance to aircraft engines and ground-based gas turbine combustors. Recent studies indicate that the injector plays a significant role in the stability of combustors by defining the flame dynamical response and setting the inlet impedance of the system. The present investigation examines the effects of combinations of injectors of two different types ($U$ and $S$) on thermoacoustic instabilities in a laboratory-scale annular combustor and compares different circumferential staging strategies. The combustor operates in a stable fashion when all injection units belong to the $S$-family, but exhibits large amplitude pressure oscillations when all these units are of the $U$-type. When the system comprises a mix of $U$- and $S$-injectors, it is possible to determine the number of $S$-injectors leading to stable operation. For a fixed proportion of $U$- and $S$-injectors, some arrangements give rise to stable operation while others do not. Results also show that introducing symmetry-breaking elements affects the system's modal dynamics. These experimental observations are interpreted in an acoustic energy balance framework used to derive an expression for the growth rate as a function of the describing functions of the flames formed by the different injectors and their respective azimuthal locations. Growth rates are determined for the different configurations and used to explain the various observations, estimate the system damping rate and predict the location of the nodal line when the standing mode prevails.