Internal combustion (IC) engines operating on fossil fuel oil provide about 25% of the world’s power (about 3000 out of 13,000 million tons oil equivalent per year—see Figure 1), and in doing so, they produce about 10% of the world’s greenhouse gas (GHG) emissions (Figure 2). Reducing fuel consumption and emissions has been the goal of engine researchers and manufacturers for years, as can be seen in the two decades of ground-breaking peer-reviewed articles published in this International Journal of Engine Research (IJER). Indeed, major advances have been made, making today’s IC engine a technological marvel. However, recently, the reputation of IC engines has been dealt a severe blow by emission scandals that threaten the ability of this technology to make significant and further contributions to the reduction of transportation sector emissions. In response, there have been proposals to replace vehicle IC engines with electric-drives with the intended goals of further reducing fuel consumption and emissions, and to decrease vehicle GHG emissions.

In this study,a novel numerical implementation for the adhesion of liquid droplets impacting normally on solid dry surfaces is presented. The advantage of this new approach, compared to the majority of existing models, is that the dynamic contact angle forming during the surface wetting process is not inserted as a boundary condition, but is derived implicitly by the induced fluid flow characteristics (interface shape) and the adhesion physics of the gas-liquid-surface interface (triple line), starting only from the advancing and receding equilibrium contact angles. These angles are required in order to define the wetting properties of liquid phases when interacting with a solid surface. The physical model is implemented as a source term in the momentum equation of a Navier-Stokes CFD flow solver as an "adhesion-like" force which acts at the triple-phase contact line as a result of capillary interactions between the liquid drop and the solid substrate. The numerical simulations capture the liquid-air interface movement by considering the volume of fluid (VOF) method and utilizing an automatic local grid refinement technique in order to increase the accuracy of the predictions at the area of interest, and simultaneously minimize numerical diffusion of the interface. The proposed model is validated against previously reported experimental data of normal impingement of water droplets on dry surfaces at room temperature. A wide range of impact velocities, i.e. Weber numbers from as low as 0.2 up to 117, both for hydrophilic (θadv=10°-70°) and hydrophobic (θadv=105°-120°) surfaces, has been examined. Predictions include in addition to droplet spreading dynamics, the estimation of the dynamic contact angle; the latter is found in reasonable agreement against available experimental measurements. It is thus concluded that theimplementation of this model is an effective approach for overcoming the need of a pre-defined dynamic contact angle law, frequently adopted as an approximate boundary condition for such simulations. Clearly, this model is mostly influential during the spreading phase for the cases of low We number impacts (We<˜80) since for high impact velocities, inertia dominates significantly over capillary forces in the initial phase of spreading.

A computational fluid dynamics cavitation model based on the Eulerian–Lagrangian approach and suitable for hole-type diesel injector nozzles is presented and discussed. The model accounts for a number of primary physical processes pertinent to cavitation bubbles, which are integrated into the stochastic framework of the model. Its predictive capability has been assessed through comparison of the calculated onset and development of cavitation inside diesel nozzle holes against experimental data obtained in real-size and enlarged models of single- and multi-hole nozzles. For the real-size nozzle geometry, high-speed cavitation images obtained under realistic injection pressures are compared against model predictions, whereas for the large-scale nozzle, validation data include images from a charge-coupled device (CCD) camera, computed tomography (CT) measurements of the liquid volume fraction and laser Doppler velocimetry (LDV) measurements of the liquid mean and root mean square (r.m.s.) velocities at different cavitation numbers (CN) and two needle lifts, corresponding to different cavitation regimes inside the injection hole. Overall, and on the basis of this validation exercise, it can be argued that cavitation modelling has reached a stage of maturity, where it can usefully identify many of the cavitation structures present in internal nozzle flows and their dependence on nozzle design and flow conditions.

The present paper focuses on the numerical simulation of the interaction of laser-generated bubbles with a free surface, including comparison of the results with instances from high-speed videos of the experiment. The Volume Of Fluid method was employed for tracking liquid and gas phases while compressibility effects were introduced with appropriate equations of state for each phase. Initial conditions of the bubble pressure were estimated through the traditional Rayleigh Plesset equation. The simulated bubble expands in a non-spherically symmetric way due to the interference of the free surface, obtaining an oval shape at the maximum size. During collapse, a jet with mushroom cap is formed at the axis of symmetry with the same direction as the gravity vector, which splits the initial bubble to an agglomeration of toroidal structures. Overall, the simulation results are in agreement with the experimental images, both quantitatively and qualitatively, while pressure waves are predicted both during the expansion and the collapse of the bubble. Minor discrepancies in the jet velocity and collapse rate are found and are attributed to the thermodynamic closure of the gas inside the bubble.

The collapse of a cavitation bubble near a rigid boundary induces a high-speed transient jet accelerating liquid onto the boundary. The shear flow produced by this event has many applications, examples of which are surface cleaning, cell membrane poration and enhanced cooling. Yet the magnitude and spatio-temporal distribution of the wall shear stress are not well understood, neither experimentally nor by simulations. Here we solve the flow in the boundary layer using an axisymmetric compressible volume-of-fluid solver from the OpenFOAM framework and discuss the resulting wall shear stress generated for a non-dimensional distance, $\unicode[STIX]{x1D6FE}=1.0$ ( $\unicode[STIX]{x1D6FE}=h/R_{max}$ , where $h$ is the distance of the initial bubble centre to the boundary, and $R_{max}$ is the maximum spherical equivalent radius of the bubble). The calculation of the wall shear stress is found to be reliable once the flow region with constant shear rate in the boundary layer is determined. Very high wall shear stresses of 100 kPa are found during the early spreading of the jet, followed by complex flows composed of annular stagnation rings and secondary vortices. Although the simulated bubble dynamics agrees very well with experiments, we obtain only qualitative agreement with experiments due to inherent experimental challenges.