In numerical linear stability investigations, the rates of change of the kinetic and thermal energy of the perturbation flow are often used to identify the dominant mechanisms by which kinetic or thermal energy is exchanged between the basic and the perturbation flow. Extending the conventional energy analysis for a single-phase Boussinesq fluid, the energy budgets of arbitrary infinitesimal perturbations to the basic two-phase liquid–gas flow are derived for an axisymmetric thermocapillary bridge when the material parameters in both phases depend on the temperature. This allows identifying individual transport terms and assessing their contributions to the instability if the basic flow and the critical mode are evaluated at criticality. The full closed-form energy budgets of linear modes have been derived for thermocapillary two-phase flow taking into account the temperature dependence of all thermophysical parameters. The influence of different approximations to the temperature dependence on the linear stability boundary of the axisymmetric flow in thermocapillary liquid bridges is tested regarding their accuracy. The general mechanism of symmetry breaking turns out to be very robust.

Gas turbine engines for fixed-wing or rotary-wing aircraft are operated in a variety of harsh weather environments ranging from arctic, volcanic zones, to desert conditions. Operation under these degraded conditions leads to the undesired entrainment of complex particulates resulting in drastic performance losses. Hence, there is a critical need to understand the governing mechanisms to inform the development of durable thermal and environmental barrier coatings. The objective of the current work is to present a novel multiscale physics-based approach to study two-phase flows that take into account the underpinning particle transport and deposition dynamics. Sessile droplet models are presented and used to compute the contact angle at high temperatures and compared with experiments. The study also investigates the sensitivity of deposition patterns to the Stokes number and the results identify local vulnerability regions. The analysis suggests that particle size distributions and the initial trajectories of the particles are critically important in predicting the final deposition pattern.

The dynamics of bubbles inertially collapsing in water near solid objects have been the subject of numerous studies in the context of cavitation erosion. While non-spherical bubble collapse, re-entrant jet dynamics and emitted shock waves have received significant interest, less is known about the temperatures thereby produced and their possible connection to damage. In this article, we use highly resolved numerical simulations of a single bubble inertially collapsing near a rigid surface to measure the temperatures produced in the fluid and estimate those in the solid, as well as to identify the responsible mechanisms. In particular, we find that elevated temperatures along the wall can be produced by one of two mechanisms, depending on the initial stand-off distance of the bubble from the wall and the driving pressure: for bubbles initially far from the wall, the shock generated by the bubble collapse is the source of the high temperature, while bubbles starting initially closer migrate towards the wall and eventually come into contact with it. A scaling is introduced to describe the maximum fluid temperature along the wall as a function of the initial stand-off distance and driving pressure. To predict the temperature of the solid, we develop a semianalytical heat transfer model, which supports recent experimental observations that elevated temperatures achieved during collapse could play a role in cavitation damage to soft heat-sensitive materials.

In two-phase flow simulations, a difficult issue is usually the treatment of surface tension effects. These cause a pressure jump that is proportional to the curvature of the interface separating the two fluids. Since the evaluation of the curvature incorporates second derivatives, it is prone to numerical instabilities. Within this work, the interface is described by a level-set method based on a discontinuous Galerkin discretization. In order to stabilize the evaluation of the curvature, a patch-recovery operation is employed. There are numerous ways in which this filtering operation can be applied in the whole process of curvature computation. Therefore, an extensive numerical study is performed to identify optimal settings for the patch-recovery operations with respect to computational cost and accuracy.

The effect of the Reynolds number on the behaviour of inertial particles in wall-bounded turbulent flows is investigated through large-scale direct numerical simulations (DNS) of particle-laden canonical channel flow spanning almost a decade in the friction Reynolds number, from $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\mathit{Re}_{\tau } = 150$ to $\mathit{Re}_{\tau } = 1000$. Lagrangian particle tracking is used to study the motion of six different particle sets, described by a Stokes number in the range $\mathit{St} = 1\text {--}1000$. At all Reynolds numbers a strong segregation in the near-wall region is observed for particles characterized by intermediate Stokes number, in the range $\mathit{St} =10\text {--}100$. The wall-normal concentration profiles of such particles collapse in inner scaling, thus suggesting the independence of the turbophoretic drift from the large-scale outer motions. This observation is also supported by the spatial organization of the suspended phase in the inner layer, which is found to be universal with the Reynolds number. The deposition rate coefficient increases with $\mathit{Re}_{\tau }$ for a given $\mathit{St}$. Suitable inner and outer scalings are proposed to collapse the deposition curves across the available ranges of Reynolds and Stokes numbers for the different deposition regimes.

The stability of a mixing layer made up of two miscible fluids, with a viscosity-stratified layer between them, is studied. The two fluids are of the same density. It is shown that unlike other viscosity-stratified shear flows, where species diffusivity is a dominant factor determining stability, species diffusivity variations over orders of magnitude do not change the answer to any noticeable degree in this case. Viscosity stratification, however, does matter, and can stabilize or destabilize the flow, depending on whether the layer of varying velocity is located within the less or more viscous fluid. By making an inviscid model flow with a slope change across the ‘viscosity’ interface, we show that viscous and inviscid results are in qualitative agreement. The absolute instability of the flow can also be significantly altered by viscosity stratification.

The onset of cavitating conditions inside the nozzle of liquid injectors is known to play a major role on spray characteristics, especially on jet penetration and break-up. In this work, we present a Direct Numerical Simulation (DNS) based on the Lattice Boltzmann Method (LBM) to study the fluid dynamic field inside the nozzle of a cavitating injector. The formation of the cavitating region is determined via a multi-phase approach based on the Shan-Chen equation of state. The results obtained by the LBM simulation show satisfactory agreement with both numerical and experimental data. In addition, numerical evidence of bubble break-up, following upon flow-induced cavitation, is also reported.

Air-induced drag reduction was investigated on a 12.9 m long flat plate test model at a free stream speed of $6. 3~\mathrm{m} ~{\mathrm{s} }^{- 1} $. Measurements of the local skin friction, phase velocity profiles (liquid and gas) and void fraction profiles were acquired at downstream distances to 11.5 m, which yielded downstream-distance-based Reynolds numbers above 80 million. Air was injected within the boundary layer behind a 13 mm backward facing step (BFS) while the incoming boundary layer was perturbed with vortex generators in various configurations immediately upstream of the BFS. Measurements confirmed that air layer drag reduction (ALDR) is sensitive to upstream disturbances, but a clean boundary layer separation line (i.e. the BFS) reduces such sensitivity. Empirical scaling of the experimental data was investigated for: (a) the critical air flux required to establish ALDR; (b) void fraction profiles; and (c) the interfacial velocity profiles. A scaling of the critical air flux for ALDR was developed from balancing shear-induced lift forces and buoyancy forces on a single bubble within a shear flow. The resulting scaling successfully collapses ALDR results from the current and past studies over a range of flow conditions and test model configurations. The interfacial velocity and void fraction profiles were acquired and scaled within the bubble drag reduction (BDR), ALDR and transitional ALDR regimes. The BDR interfacial velocity profile revealed that there was slip between phases. The ALDR results showed that the air layer thickness was nominally three-quarters of the total volumetric flux (per unit span) of air injected divided by the free stream speed. Furthermore, the air layer had an average void fraction of 0.75 and a velocity of approximately 0.2 times the free stream speed. Beyond the air layer was a bubbly mixture that scaled in a similar fashion to the BDR results. Transitional ALDR results indicate that this regime was comprised of intermittent generation and subsequent fragmentation of an air layer, with the resulting drag reduction determined by the fraction of time that an air layer was present.

The dynamics of isolated air bubbles crossing the horizontal interface separating two Newtonian immiscible liquids initially at rest are studied both experimentally and computationally. High-speed video imaging is used to obtain a detailed evolution of the various interfaces involved in the system. The size of the bubbles and the viscosity contrast between the two liquids are varied by more than one and four orders of magnitude, respectively, making it possible to obtain bubble shapes ranging from spherical to toroidal. A variety of flow regimes is observed, including that of small bubbles remaining trapped at the fluid–fluid interface in a film-drainage configuration. In most cases, the bubble succeeds in crossing the interface without being stopped near its undisturbed position and, during a certain period of time, tows a significant column of lower fluid which sometimes exhibits a complex dynamics as it lengthens in the upper fluid. Direct numerical simulations of several selected experimental situations are performed with a code employing a volume-of-fluid type formulation of the incompressible Navier–Stokes equations. Comparisons between experimental and numerical results confirm the reliability of the computational approach in most situations but also points out the need for improvements to capture some subtle but important physical processes, most notably those related to film drainage. Influence of the physical parameters highlighted by experiments and computations, especially that of the density and viscosity contrasts between the two fluids and of the various interfacial tensions, is discussed and analysed in the light of simple models and available theories.

Homogeneous condensation of vapours mixed with a carrier gas in the stagnation point boundary layer flow near a cold wall is considered. There is a condensation region near the wall with supersaturated vapour. Assuming that the surface tension times the molecular area is much larger than the thermal energy far from the wall, droplets are nucleated exclusively in a narrow nucleation layer where the Zeldovich flux of clusters surpassing the critical nucleus size is at a maximum. The vapour condenses in the free molecular regime on the droplets, which are thermophoretically attracted to the wall. Unlike the narrow condensation region for heterogeneous condensation on solid particles, in the case of homogeneous condensation the condensation region is wide even when the rate of vapour scavenging by droplets is large. A singular perturbation theory of homogeneous vapour condensation in boundary layer flow approximates very well the vapour and droplet density profiles, the nucleation layer and the deposition rates at the wall for wide ranges of the wall temperature and the scavenging parameter . A key point in the theory is to select a trial vapour number density profile among a one parameter family of profiles between an upper and a lower bound. The maximum of the Zeldovich flux for supercritical nuclei provides the approximate location of the nucleation layer and an approximate droplet density profile. Then the condensate number of molecules and the vapour density profile are calculated by matched asymptotic expansions that also yield the deposition rates. For sufficiently large wall temperatures, a more precise corrected asymptotic theory is given.

The simulation of multiphase flows is an outstanding challenge, due to the inherent complexity of the underlying physical phenomena and to the fact that multiphase flows are very diverse in nature, and so are the laws governing their dynamics. In the last two decades, a new class of mesoscopic methods, based on minimal lattice formulation of Boltzmann kinetic equation, has gained significant interest as an efficient alternative to continuum methods based on the discretisation of the NS equations for non ideal fluids. In this paper, three different multiphase models based on the lattice Boltzmann method (LBM) are discussed, in order to assess the capability of the method to deal with multiphase flows on a wide spectrum of operating conditions and multiphase phenomena. In particular, the range of application of each method is highlighted and its effectiveness is qualitatively assessed through comparison with numerical and experimental literature data.