346 research outputs found
Cavitation Induction by Projectile Impacting on a Water Jet
The present paper focuses on the simulation of the high-velocity impact of a projectile impacting on a water-jet, causing the onset, development and collapse of cavitation. The simulation of the fluid motion is carried out using an explicit, compressible, density-based solver developed by the authors using the OpenFOAM library. It employs a barotropic two-phase flow model that simulates the phase-change due to cavitation and considers the co-existence of non-condensable and immiscible air. The projectile is considered to be rigid while its motion through the computational domain is modelled through a direct-forcing Immersed Boundary Method. Model validation is performed against the experiments of Field et al. [Field, J., Camus, J. J., Tinguely, M., Obreschkow, D., Farhat, M., 2012. Cavitation in impacted drops and jets and the effect on erosion damage thresholds. Wear 290–291, 154–160. doi:10.1016/j.wear.2012.03.006. URL http://www.sciencedirect.com/science/article/pii/S0043164812000968 ], who visualised cavity formation and shock propagation in liquid impacts at high velocities. Simulations unveil the shock structures and capture the high-speed jetting forming at the impact location, in addition to the subsequent cavitation induction and vapour formation due to refraction waves. Moreover, model predictions provide quantitative information and a better insight on the flow physics that has not been identified from the reported experimental data, such as shock-wave propagation, vapour formation quantity and induced pressures. Furthermore, evidence of the Richtmyer-Meshkov instability developing on the liquid-air interface are predicted when sufficient dense grid resolution is utilised
Recommended from our members
Numerical investigation of bubble dynamics using tabulated data
An explicit density-based solver of the compressible Euler equations suitable for cavitation simulations is presented, using the full Helmholtz energy Equation of State (EoS) for n-Dodecane. Tabulated data are derived from this EoS in order to calculate the thermodynamic properties of the liquid, vapour and mixture composition during cavitation. For determining thermodynamic properties from the conservative variable set, bilinear interpolation is employed; this results to significantly reduced computational cost despite the complex thermodynamics model incorporated. The latter is able to predict the temperature variation of both the liquid and the vapour phases. The methodology uses a Mach number consistent numerical flux, suitable for subsonic up to supersonic flow conditions. Finite volume discretization is employed in conjunction with a second order Runge-Kutta time integration scheme. The numerical method is validated against the Riemann problem, comparing it with the exact solution which has been derived in the present work for an arbitrary EoS. Further validation is performed against the well-known Rayleigh collapse of a pure vapour bubble. It is then used for the simulation of a 2-D axisymmetric n-Dodecane vapour bubble collapsing in the proximity of a flat wall placed at different locations from the centre of the bubble. The predictive capability of the incorporated Helmholtz EoS is assessed against the widely used barotropic EoS and the non-isothermal Homogeneous Equilibrium Mixture (HEM)
Recommended from our members
Smoothed particle hydrodynamics simulation of a laser pulse impact onto a liquid metal droplet
The impact of a laser pulse onto a liquid metal droplet is numerically investigated by utilising a weakly compressible single phase model; the thermodynamic closure is achieved by the Tait equation of state (EoS) for the liquid metal. The smoothed particle hydrodynamics (SPH) method, which has been employed in the arbitrary Lagrangian Eulerian (ALE) framework, offers numerical efficiency, compared to grid related discretization methods. The latter would require modelling not only of the liquid metal phase, but also of the vacuum, which would necessitate special numerical schemes, suitable for high density ratios. In addition, SPH-ALE allows for the easy deformation handling of the droplet, compared to interface tracking methods where strong mesh deformation and most likely degenerate cells occur. Then, the laser-induced deformation of the droplet is simulated and cavitation formation is predicted. The ablation pattern due to the emitted shock wave and the two low pressure lobes created in the middle of the droplet because of the rarefaction waves are demonstrated. The liquid metal droplet is subject to material rupture, when the shock wave, the rarefaction wave and the free surface interact. Similar patterns regarding the wave dynamics and the hollow structure have been also noticed in prior experimental studies
Recommended from our members
Heating Effects During Bubble Collapse Using Tabulated Data
An explicit density-based solver for the compressible Navier-Stokes equations able to simulate cavitating flows has been developed and utilised for the simulation of collapsing vapour bubbles. Phase-change is considered by employing the homogeneous equilibrium model (HEM). The wide variation of Mach numbers between the liquid, vapour and mixture regimes is tackled by a Mach consistent numerical flux, suitable for subsonic up to supersonic flow conditions. Time discretisation is performed using a second order low storage Runge-Kutta scheme. Thermodynamic closure is achieved by utilising the Helmholtz energy equation of state (EoS), making feasible simulation of conditions at subcritical and supercritical regions considering the variations of liquid and vapour temperatures during bubble collapse. In order to reduce the computational cost associated with the solution of the Helmholtz EoS at each time step, a tabulated data technique has been followed. The unstructured thermodynamic table, containing the thermodynamic properties derived from the Helmholtz EoS, has been constructed for n-dodecane, which has been the considered as the working fluid. The efficiency of the method is enhanced by a static linked-list algorithm for searching among the elements of the table. In addition, a finite element bilinear interpolation is used for approximating the unknown thermodynamic properties. After validating the numerical method, parametric studies considering 2-D axisymmetric vaporous bubble collapse in the proximity of a wall have been performed at conditions realised in micro-orifice flow passages. The temperature and pressure changes on the wall are estimated as function of the surrounding liquid pressure, the initial bubble radius and the location of the wall from the center of the initial bubble, revealing the expected range of variation as function on the set parameters
Recommended from our members
A tabulated data technique for cryogenic two-phase flows
Flashing flows of liquid oxygen (LOX) are prevalent in space applications, where LOX can be used as rocket engine propellant [1]–[3]. Towards this direction, the cryogenic flow in a converging-diverging nozzle has been investigated in the present study by utilising real fluid thermodynamics under the homogeneous equilibrium mixture (HEM) assumption for the LOX. A tabulated data method for the Helmholtz energy equation of state (EoS) has been developed in OpenFOAM (OF) [4] and has been incorporated into an explicit density based solver. Due to the wide variation of the speed of sound and consequently of the Mach number noticed in the liquid, vapour and mixture phases, a Mach consistent numerical flux has been employed suitable for subsonic up to supersonic flow conditions [5]. Since the Helmholtz EoS is computationally inefficient compared to simplified EoS, an ad-hoc thermodynamic table containing all the thermodynamic properties for the LOX has been created and stored prior entering the time loop [6], accompanied by a static linked-list algorithm for reducing the search time. Once the thermodynamic element of the table which satisfies the values of the density and internal energy as predicted from the numerical solution of the Navier-Stokes equations is identified, the unknown thermodynamic properties are approximated by a finite element interpolation [7]. The numerical method has been firstly validated against the Riemann problem at similar cryogenic flow conditions. Then, 2-D axisymmetric simulations of the phase-change process in a converging-diverging nozzle are performed and compared with prediction from other numerical tools as well as experimental data. It is concluded that the results are satisfactory while the applicability of the Helmholtz EoS to LOX simulations is demonstrated. This suggests that the proposed methodology can be utilized for the simulation of flashing flows
- …