The present work is devoted to the analysis of drop impact on a deep liquid pool, focusing on the high-energy splashing regimes caused by large raindrops at high velocities. Such cases are characterized by short time scales and complex mechanisms, thus they have received very little attention until now. The BASILISK open-source solver is used to perform three-dimensional direct numerical simulations. The capabilities of octree adaptive mesh refinement techniques enable capturing of the small-scale features of the flow, while the volume of fluid approach combined with a balanced-force surface-tension calculation is applied to advect the volume fraction of the liquids and reconstruct the interfaces. The numerical results compare well with experimental visualizations: both the evolution of crown and cavity, the emanation of ligaments, the formation of bubble canopy and the growth of a downward-moving spiral jet that pierces through the cavity bottom, are correctly reproduced. Reliable quantitative agreements are also obtained regarding the time evolution of rim positions, cavity dimensions and droplet distributions through an observation window. Furthermore, simulation gives access to various aspects of the internal flows, which allows us to better explain the observed physical phenomena. Details of the early-time dynamics of bubble ring entrapment and splashing performance, the formation/collapse of bubble canopy and the spreading of drop liquid are discussed. The statistics of droplet size show the bimodal distribution in time, corroborating distinct primary mechanisms of droplet production at different stages.

The present study focuses on two-dimensional direct numerical simulations of shallow-water breaking waves, specifically those generated by a wave plate at constant water depths. The primary objective is to quantitatively analyse the dynamics, kinematics and energy dissipation associated with wave breaking. The numerical results exhibit good agreement with experimental data in terms of free-surface profiles during wave breaking. A parametric study was conducted to examine the influence of various wave properties and initial conditions on breaking characteristics. According to research on the Bond number ( $Bo$ , the ratio of gravitational to surface tension forces), an increased surface tension leads to the formation of more prominent parasitic capillaries at the forwards face of the wave profile and a thicker plunging jet, which causes a delayed breaking time and is tightly correlated with the main cavity size. A close relationship between wave statistics and the initial conditions of the wave plate is discovered, allowing for the classification of breaker types based on the ratio of wave height to water depth, $H/d$ . Moreover, an analysis based on inertial scaling arguments reveals that the energy dissipation rate due to breaking can be linked to the local geometry of the breaking crest $H_b/d$ , and exhibits a threshold behaviour, where the energy dissipation approaches zero at a critical value of $H_b/d$ . An empirical scaling of the breaking parameter is proposed as $b = a(H_b/d - \chi _0)^n$ , where $\chi _0 = 0.65$ represents the breaking threshold and $n = 1.5$ is a power law determined through the best fit to the numerical results.

The effect of nonsinusoidal trajectory on the propulsive performances and the vortex shedding process behind a flapping airfoil is investigated in this study. A movement of a rigid NACA0012 airfoil undergoing a combined heaving and pitching motions at low Reynolds number (Re = 11,000) is considered. An elliptic function with an adjustable parameter S (flattening parameter) is used to realize various nonsinusoidal trajectories of both motions. The two-dimensional (2D) unsteady and incompressible Navier–Stokes equation governing the flow over the flapping airfoil are resolved using the commercial software starccm+. It is shown that the nonsinusoidal flapping motion has a major effect on the propulsive performances of the flapping airfoil. Although the maximum propulsive efficiency is always achievable with sinusoidal trajectories, nonsinusoidal trajectories are found to considerably improve performance: a 110% increase of the thrust force was obtained in the best studied case. This improvement is mainly related to the modification of the heaving motion, more specifically the increase of the heaving speed at maximum pitching angle of the foil. The analysis of the flow vorticity and wake structure also enables to explain the drop of the propulsive efficiency for nonsinusoidal trajectories.

Measurements are processed on a centrifugal pump model, which works with air and performs with the vane-island type diffuser of a real hydraulic pump, under five flow rates to investigate the internal flow characteristics and their influence on overall pump performance. The mean flow characteristics inside the diffuser are determined by using a miniature three-hole probe connected to an online data acquisition system. The flow structure at the inlet section of the diffuser is analyzed in detail, with a focus on the local pressure loss inside the vaneless gap and incidence angle distributions along the hub-to-shroud direction of the diffuser. Some existing calculations, including leakage effects, are used to evaluate the pressure recovery downstream of the impeller. Furthermore, particle image velocimetry (PIV) measurement results are obtained to help analyze the flow characteristics inside the vane-island diffuser. Each PIV measuring plane is related to one particular diffuser blade-to-blade channel and is analyzed by using the time-averaged method according to seven different relative positions of the impeller. Measurement results show that main loss is produced inside the vaneless part of the diffuser at low flow rates, which might have been caused by the strong rotor–stator interaction. When the impeller flow rate is greater than the diffuser design flow rate, a large fluctuating separated region occurs after the throat of the diffuser on the pressure side. Mean loss originates from the unsteady pressure downstream of the diffuser throat. For better characterization of the separations observed in previous experimental studies, complementary unsteady static pressure measurement campaigns have been conducted on the diffuser blade wall. The unsteadiness revealed by these measurements, as well as theirs effects on the diffuser performance, was then studied.

The present paper is devoted to the analysis of the various instabilities of cavitation attached to a two-dimensional (2D) profile. Time resolved stereo particle image velocimetry was conducted in a small-scale 2D Venturi type section, in different vertical planes in the streamwise direction, located at varying positions in the depth of the channel. These experiments enabled to obtain the time evolution of the three components of the velocity field in the cavitation area and to derive the time-averaged gradients in the spanwise direction. Test cases at various Reynolds numbers were conducted, maintaining either the pressure or the cavitation number constant, to discuss the impact of these parameters on the flow. Then, attention was focused on three distinct flow dynamics, namely, sheet cavitation, where no large-scale instability can be detected, single cloud cavitation, where a large cloud of vapor is shed periodically at the rear of the cavity, and multi-cloud cavitation, where the process is more complex, as more than one clouds are shed downstream. Data reveal that the structure and the structure of the re-entrant jet, which is one of the primary mechanisms of cloud cavitation, are more complex than reported in the previous studies. Although the jet can be detected as an intermittent low speed reverse flow in the streamwise direction, it is actually made of successive vortices about the channel depth, which are convected downstream while expanding in the vertical direction, causing the cavity lift and, thus, contributing to its final split and the cloud shedding.