Evaluation of Cavitation Erosion Intensity in a Microscale Nozzle Using Eulerian–Lagrangian Bubble Dynamic Simulation

2019 ◽  
Vol 141 (6) ◽  
Author(s):  
Masoud Khojasteh-Manesh ◽  
Miralam Mahdi

In the present study, cavitation erosion is investigated by implementing an Eulerian–Lagrangian approach. Three-dimensional two-phase flow is simulated in a microscale nozzle using Reynolds-averaged Navier–Stokes (RANS) solver along with realizable k−ε turbulence model and Schnerr–Sauer cavitation model. The numerical results are in agreement with experimental observations. A modified form of Rayleigh–Plesset–Keller–Herring equation along with bubble motion equation is utilized to simulate bubble dynamics. Average values of mixture properties over bubble surface are used instead of bubble-center values in order to account for nonuniformities around the bubble. A one-way coupling method is used between Lagrangian analysis and RANS solution. The impact pressure resulted from bubble collapse is calculated for evaluation of erosion in diesel and soy methyl ester (SME) biodiesel in different situations. The results show that the initial size of the bubbles is an important factor for determining the intensity of erosion. So, the bubbles erosive power increases when their initial radius increases. It is also found that the intensity of erosion in diesel is much higher than that of biodiesel and this is because of the differences in fuels properties, especially in viscosity and vapor pressure. The effect of bubbles initial position on erosion intensity is also investigated in this study, and it is found that bubbles with the highest distance from sheet cavity termination have the highest contribution in erosion rate.

Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 208
Author(s):  
Peter Reinke ◽  
Jan Ahlrichs ◽  
Tom Beckmann ◽  
Marcus Schmidt

The volume-of-flow method combined with the Rayleigh–Plesset equation is well established for the computation of cavitation, i.e., the generation and transportation of vapor bubbles inside a liquid flow resulting in cloud, sheet or streamline cavitation. There are, however, limitations, if this method is applied to a restricted flow between two adjacent walls and the bubbles’ size is of the same magnitude as that of the clearance between the walls. This work presents experimental and numerical results of the bubble generation and its transportation in a Couette-type flow under the influence of shear and a strong pressure gradient which are typical for journal bearings or hydraulic seals. Under the impact of variations of the film thickness, the VoF method produces reliable results if bubble diameters are less than half the clearance between the walls. For larger bubbles, the wall contact becomes significant and the bubbles adopt an elliptical shape forced by the shear flow and under the influence of a strong pressure gradient. Moreover, transient changes in the pressure result in transient cavitation, which is captured by high-speed imaging providing material to evaluate transient, three-dimensional computations of a two-phase flow.


Lubricants ◽  
2021 ◽  
Vol 9 (4) ◽  
pp. 44
Author(s):  
Christian Ziese ◽  
Cornelius Irmscher ◽  
Steffen Nitzschke ◽  
Christian Daniel ◽  
Elmar Woschke

The vibration behaviour of turbocharger rotors is influenced by the acting loads as well as by the type and arrangement of the hydrodynamic bearings and their operating condition. Due to the highly non-linear bearing behaviour, lubricant film-induced excitations can occur, which lead to sub-synchronous rotor vibrations. A significant impact on the oscillation behaviour is attributed to the pressure distribution in the hydrodynamic bearings, which is influenced by the thermo-hydrodynamic conditions and the occurrence of outgassing processes. This contribution investigates the vibration behaviour of a floating ring supported turbocharger rotor. For detailed modelling of the bearings, the Reynolds equation with mass-conserving cavitation, the three-dimensional energy equation and the heat conduction equation are solved. To examine the impact of outgassing processes and thrust bearing on the occurrence of sub-synchronous rotor vibrations separately, a variation of the bearing model is made. This includes run-up simulations considering or neglecting thrust bearings and two-phase flow in the lubrication gap. It is shown that, for a reliable prediction of sub-synchronous vibrations, both the modelling of outgassing processes in hydrodynamic bearings and the consideration of thrust bearing are necessary.


2019 ◽  
Vol 30 (8) ◽  
pp. 3845-3877 ◽  
Author(s):  
Eric Goncalves Da Silva ◽  
Philippe Parnaudeau

Purpose The purpose of this paper is to quantify the relative importance of the multiphase model for the simulation of a gas bubble impacted by a normal shock wave in water. Both the free-field case and the collapse near a wall are investigated. Simulations are performed on both two- and three-dimensional configurations. The main phenomena involved in the bubble collapse are illustrated. A focus on the maximum pressure reached during the collapse is proposed. Design/methodology/approach Simulations are performed using an inviscid compressible homogeneous solver based on different systems of equations. It consists in solving different mixture or phasic conservation laws and a transport-equation for the gas volume fraction. Three-dimensional configurations are considered for which an efficient massively parallel strategy was developed. The code is based on a finite volume discretization for which numerical fluxes are computed with a Harten, Lax, Van Leer, Contact (HLLC) scheme. Findings The comparison of three multiphase models is proposed. It is shown that a simple four-equation model is well-suited to simulate such strong shock-bubble interaction. The three-dimensional collapse near a wall is investigated. It is shown that the intensity of pressure peaks on the wall is drastically increased (more than 200 per cent) in comparison with the cylindrical case. Research limitations/implications The study of bubble collapse is a key point to understand the physical mechanism involved in cavitation erosion. The bubble collapse close to the wall has been addressed as the fundamental mechanism producing damage. Its general behavior is characterized by the formation of a water jet that penetrates through the bubble and the generation of a blast wave during the induced collapse. Both the jet and the blast wave are possible damaging mechanisms. However, the high-speed dynamics, the small spatio-temporal scales and the complicated physics involved in these processes make any theoretical and experimental approach a challenge. Practical implications Cavitation erosion is a major problem for hydraulic and marine applications. It is a limiting point for the conception and design of such components. Originality/value Such a comparison of multiphase models in the case of a strong shock-induced bubble collapse is clearly original. Usually models are tested separately leading to a large dispersion of results. Moreover, simulations of a three-dimensional bubble collapse are scarce in the literature using such fine grids.


Author(s):  
Agus P. Sasmito ◽  
Tariq Shamim ◽  
Erik Birgersson ◽  
Arun S. Mujumdar

In open-cathode polymer electrolyte fuel cell (PEFC) stacks, a significant temperature rise can exist due to insufficient cooling, especially at higher current densities. To improve stack thermal management while reducing the cost of cooling, we propose a forced air-convection open-cathode fuel cell stack with edge cooling (fins). The impact of the edge cooling is studied via a mathematical model of the three-dimensional two-phase flow and the associated conservation equations of mass, momentum, species, energy, and charge. The model includes the stack, ambient, fan, and fins used for cooling. The model results predict better thermal management and stack performance for the proposed design as compared to the conventional open-cathode stack design, which shows potential for practical applications. Several key design parameters—fin material and fin geometry—are also investigated with regard to the stack performance and thermal management.


2012 ◽  
Vol 135 (1) ◽  
Author(s):  
Naoya Ochiai ◽  
Yuka Iga ◽  
Motohiko Nohmi ◽  
Toshiaki Ikohagi

Cavitation erosion is a material damage phenomenon caused by the repeated application of impulsive pressure on a material surface induced by bubble collapse, and the establishment of a method by which to numerically predict cavitation erosion is desired. In the present study, a numerical quantitative prediction method of cavitation erosion in a cavitating flow is proposed. In the present method, a one-way coupled analysis of a cavitating flow field based on a gas-liquid two-phase Navier–Stokes equation (Eulerian) and bubbles in the cavitating flow by bubble dynamics (Lagrangian) is used to treat temporally and spatially different scale phenomena, such as the macroscopic phenomenon of a cavitating flow and the microscopic phenomenon of bubble collapse. Impulsive pressures acting on a material surface are evaluated based on the bubble collapse position, time, and intensity, and the erosion rate is quantitatively predicted using an existing material-dependent relationship between the impulsive energy (square of the impulsive force) and the maximum erosion rate. The erosion rate on a NACA0015 hydrofoil surface in an unsteady transient cavitating flow is predicted by the proposed method. The distribution of the predicted erosion rate corresponds qualitatively to the distribution of an experimental surface roughness increment of the same hydrofoil. Furthermore, the predicted erosion rate considering the bubble nuclei distribution is found to be of the same order of magnitude as the actual erosion rate, which indicates that considering bubble nuclei distribution is important for the prediction of cavitation erosion and that the present prediction method is valid to some degree.


2014 ◽  
Vol 755 ◽  
pp. 142-175 ◽  
Author(s):  
Chao-Tsung Hsiao ◽  
A. Jayaprakash ◽  
A. Kapahi ◽  
J.-K. Choi ◽  
Georges L. Chahine

AbstractMaterial pitting from cavitation bubble collapse is investigated numerically including two-way fluid–structure interaction (FSI). A hybrid numerical approach which links an incompressible boundary element method (BEM) solver and a compressible finite difference flow solver is applied to capture non-spherical bubble dynamics efficiently and accurately. The flow codes solve the fluid dynamics while intimately coupling the solution with a finite element structure code to enable simulation of the full FSI. During bubble collapse high impulsive pressures result from the impact of the bubble re-entrant jet on the material surface and from the collapse of the remaining bubble ring. A pit forms on the material surface when the impulsive pressure is large enough to result in high equivalent stresses exceeding the material yield stress. The results depend on bubble dynamics parameters such as the size of the bubble at its maximum volume, the bubble standoff distance from the material wall, and the pressure driving the bubble collapse. The effects of these parameters on the re-entrant jet, the following bubble ring collapse pressure, and the generated material pit characteristics are investigated.


Author(s):  
Agus Pulung Sasmito ◽  
Tariq Shamim ◽  
Erik Birgersson ◽  
Arun Sadashiv Mujumdar

In open-cathode polymer electrolyte fuel cell (PEFC) stacks, a significant temperature rise can exist due to insufficient cooling, especially at higher current densities. To improve stack thermal management whilst reducing the cost for cooling, we propose a forced air-convection open-cathode fuel cell stack with edge cooling (fins). The impact of the edge cooling is studied via mathematical model of the three-dimensional two-phase flow and associated conservation equations of mass, momentum, species, energy and charge. The model includes stack, ambient, fan and fins used for cooling. The model results predict better thermal management and stack performance for the proposed design as compared to the conventional open-cathode stack design, which shows potential for practical application. Several key design parameters — fin material and fin geometry — are also investigated with regards to the stack performance and thermal management.


1998 ◽  
Vol 120 (1) ◽  
pp. 179-185 ◽  
Author(s):  
Shuji Hattori ◽  
Hiroyuki Mori ◽  
Tsunenori Okada

In order to evaluate the quantitative cavitation-erosion resistance of materials, a pressure-detector-installed specimen was developed, which can measure both the impact load produced by cavitation bubble collapse and the volume loss simultaneously. Test specimens (pressure-detection rod) used were nine kinds of metals and were exposed to vibratory cavitation. A linear relation was obtained for all materials between the accumulated impact energy ∑Fi2 calculated from the distribution of impact loads and the volume loss, independent of test conditions. Impact energy accumulated during the incubation period and the energy for a unit material removal in steady-state period were obtained from the relation. These values are very Important concerning quantitative erosion resistance evaluation. That is, when the distribution of impact loads is acquired for different cavitation conditions, the volume loss can be estimated. This idea was applied to the venturi cavitation erosion. The experimental results for venturi test corresponds well with the prediction using these impact energy values. It was concluded that the quantitative impact energy values of materials can be determined independent of the apparatus and the test condition by using the newly developed pressure-detector-installed specimen.


2007 ◽  
Vol 135 (2) ◽  
pp. 549-566 ◽  
Author(s):  
Shu-Hua Chen

Abstract Three observational datasets of Hurricane Isidore (in 2002) were analyzed and compared: the Special Sensor Microwave Imager (SSM/I), the Quick Scatterometer (QuikSCAT) winds, and dropsonde winds. SSM/I and QuikSCAT winds were on average about 1.9 and 0.3 m s−1 stronger, respectively, than dropsonde winds. With more than 20 000 points of data, SSM/I wind speed was about 2.2 m s−1 stronger than QuikSCAT. Comparison of the wind direction observed by QuikSCAT with those from the dropsondes showed that the quality of QuikSCAT data is good. The effect of assimilating SSM/I wind speeds and/or QuikSCAT wind vectors for the analysis of Hurricane Isidore was assessed using the fifth-generation Pennsylvania State University–National Center for Atmospheric Research (PSU–NCAR) Mesoscale Model (MM5) and its three-dimensional variational data assimilation system. For the Hurricane Isidore case study, it was found that the assimilation of either satellite winds strengthened the cyclonic circulation in the analysis. However, the increment of the QuikSCAT wind analysis is more complicated than that from the SSM/I analysis due to the correction of the storm location, a positive result from the assimilation of wind vectors. The increase in low-level wind speeds enhanced the air–sea interaction processes and improved the simulated intensity for Isidore. In addition, the storm structure was better simulated. Assimilation of QuikSCAT wind vectors clearly improved simulation of the storm track, in particular during the later period of the simulation, but lack of information about the wind direction from SSM/I data prevented it from having much of an effect. Assessing the assimilation of QuikSCAT wind speed versus wind vector data confirmed this hypothesis. The track improvement partially resulted from the relocation of the storm’s initial position after assimilation of the wind vectors. For this case study, it was found that the assimilation of SSM/I or QuikSCAT data had the greatest impact on the Hurricane Isidore simulation during the first 2 days.


Author(s):  
Wen xi Tian ◽  
Sui-zheng Qiu ◽  
Guang-hui Su ◽  
Yuki Ishiwatari ◽  
Yoshiaki Oka

Single vacuum bubble collapse in subcooled water has been simulated using the moving particle semi-implicit (MPS) method in the present study. The liquid is described using moving particles, and the bubble-liquid interface was set to be the vacuum pressure boundary without interfacial heat mass transfer. The topological shape of the vacuum bubble is determined according to the location of interfacial particles. The time dependent bubble diameter, interfacial velocity, and bubble collapse time were obtained within a wide parametric range. Comparison with Rayleigh’s prediction indicates a good consistency, which validates the applicability and accuracy of the MPS method. The potential void-induced water hammer pressure pulse was also evaluated, which is instructive for the cavitation erosion study. The present paper discovers fundamental characteristics of vacuum bubble hydrodynamics, and it is also instructive for further applications of the MPS method to complicated bubble dynamics.


Sign in / Sign up

Export Citation Format

Share Document