Numerical Solution on Spherical Vacuum Bubble Collapse Using MPS Method

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.

Author(s):  
Wenxi Tian ◽  
Suizheng Qiu ◽  
Guanghui Su ◽  
Yoshiaki Oka

In this study, the vacuum bubble collapse in liquid pool has been simulated using MPS code. The liquid is described using moving particles and the bubble-liquid interface was set to be vacuum pressure boundary without interfacial heat mass transfer. The motion and location of interfacial particles can be competent in configurating the topological shape of vacuum bubble. The time dependent bubble diameter, interfacial velocity and bubble collapse time were obtained under wide parametric range. The comparison with Rayleigh’s prediction showed a good agreement which validates the applicability and accuracy on MPS method in solving present momentum problems. The potential void-induced water hammer pressure pulse was also evaluated which is instructive for cavitaion erosion study. The bubble collapse with noncondensable gas has also been simulated and the rebound phenomenon was successfully captured which is similar with vapor-filled cavitation phenomenon. The present study exhibits some fundamental characteristics of vacuum bubble hydrodynamics and it is also expected to be instructive for further applications of MPS method to in complicated bubble dynamics.


2013 ◽  
Vol 135 (1) ◽  
Author(s):  
R. Fortes-Patella ◽  
G. Challier ◽  
J. L. Reboud ◽  
A. Archer

An original approach based on energy balance between vapor bubble collapse, emitted pressure wave, and neighboring solid wall response was proposed, developed, and tested to estimate the aggressiveness of cavitating flows. In the first part of the work, to improve a prediction method for cavitation erosion (Fortes-Patella and Reboud, 1998, “A New Approach to Evaluate the Cavitation Erosion Power,” ASME J. Fluids Eng., 120(2), pp. 335–344; Fortes-Patella and Reboud, 1998, “Energetical Approach and Impact Efficiency in Cavitation Erosion,” Proceedings of Third International Symposium on Cavitation, Grenoble, France), we were interested in studying the pressure waves emitted during bubble collapse. The radial dynamics of a spherical vapor/gas bubble in a compressible and viscous liquid was studied by means of Keller's and Fujikawa and Akamatsu's physical models (Prosperetti, 1994, “Bubbles Dynamics: Some Things we did not Know 10 Years Ago,” Bubble Dynamics and Interface Phenomena, Blake, Boulton-Stone, Thomas, eds., Kluwer Academic Publishers, Dordrecht, the Netherlands, pp. 3–15; Fujikawa and Akamatsu, 1980, “Effects of Non-Equilibrium Condensation of Vapor on the Pressure Wave Produced by Collapse of a Bubble in Liquid,” J. Fluid Mech., 97(3), pp. 481–512). The pressure amplitude, the profile, and the energy of the pressure waves emitted during cavity collapses were evaluated by numerical simulations. The model was validated by comparisons with experiments carried out at Laboratoire Laser, Plasma et Procédés Photoniques (LP3-IRPHE) (Marseille, France) with laser-induced bubble (Isselin et al., 1998, “Investigations of Material Damages Induced by an Isolated Vapor Bubble Created by Pulsed Laser,” Proceedings of Third International Symposium on Cavitation, Grenoble, France; Isselin et al., 1998, “On Laser Induced Single Bubble Near a Solid Boundary: Contribution to the Understanding of Erosion Phenomena,” J. Appl. Phys., 84(10), pp. 5766–5771). The efficiency of the first collapse ηwave/bubble (defined as the ratio between pressure wave energy and initial bubble potential energy) was evaluated for different bubble collapses. For the cases considered of collapse in a constant-pressure field, the study pointed out the strong influence of the air contents on the bubble dynamics, on the emitted pressure wave characteristics, and on the collapse efficiency. In the second part of the study, the dynamic response and the surface deformation (i.e., pit profile and pit volume) of various materials exposed to pressure wave impacts was simulated making use of a 2D axisymmetric numerical code simulating the interaction between pressure wave and an elastoplastic solid. Making use of numerical results, a new parameter β (defined as the ratio between the pressure wave energy and the generated pit volume) was introduced and evaluated for three materials (aluminum, copper, and stainless steel). By associating numerical simulations and experimental results concerning pitted samples exposed to cavitating flows (volume damage rate), the pressure wave power density and the flow aggressiveness potential power were introduced. These physical properties of the flow characterize the cavitation intensity and can be related to the flow hydrodynamic conditions. Associated to β and ηwave/bubble parameters, these power densities appeared to be useful tools to predict the cavitation erosion power.


2011 ◽  
Vol 133 (4) ◽  
Author(s):  
Jinfeng Wu ◽  
Vijay K. Dhir

During phase change at the bubble-liquid interface, under subcooled boiling conditions, noncondensable gases dissolved in the liquid will be injected into the bubble along with vapor. Due to heat transfer into subcooled liquid, vapor will condense in the upper regions of the bubble while noncondensables will continue to accumulate. Subsequently, thermocapillary convection caused by nonuniform saturation temperature at the interface may occur. The aim of this work is to investigate the effects of noncondensables on heat transfer and bubble dynamics. The numerical results show that the effects of noncondensables on 5°C subcooled boiling of water are minor in terms of the equilibrium bubble diameter and overall Nusselt number. However, induced flow pattern around the bubble is altered, especially under reduced gravity conditions.


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.


2020 ◽  
Vol 2020 ◽  
pp. 1-14
Author(s):  
Chunhui Wang ◽  
Chunyu Guo ◽  
Fenglei Han

Modified 3D Moving Particle Semi-Implicit (MPS) method is used to complete the numerical simulation of the fluid sloshing in LNG tank under multidegree excitation motion, which is compared with the results of experiments and 2D calculations obtained by other scholars to verify the reliability. The cubic spline kernel functions used in Smoothed Particle Hydrodynamics (SPH) method are adopted to reduce the deviation caused by consecutive two times weighted average calculations; the boundary conditions and the determination of free surface particles are modified to improve the computational stability and accuracy of 3D calculation. The tank is under forced multidegree excitation motion to simulate the real conditions of LNG ships, the pressures and the free surfaces at different times are given to verify the accuracy of 3D simulation, and the free surface and the splashed particles can be simulated more exactly.


2015 ◽  
Vol 18 (4) ◽  
pp. 1363-1377 ◽  
Author(s):  
Davi Teodoro Fernandes ◽  
Liang-Yee Cheng ◽  
Eric Henrique Favero ◽  
Kazuo Nishimoto

Author(s):  
Kailun Guo ◽  
Ronghua Chen ◽  
Suizheng Qiu ◽  
Wenxi Tian ◽  
Guanghui Su ◽  
...  

Multiphase flow widely exists in the nature and engineering. The two-phase flow is the highlight of the studies about the flow in the vessel and steam explosion in nuclear severe accidents. The Moving Particle Semi-implicit (MPS) method is a fully-Lagrangian particle method without grid mesh which focuses on tracking the single particle and concerns with its movement. It has advantages in tracking complex multiphase flows compared with gird methods, and thus shows great potential in predicting multiphase flows. The objective of this thesis is to develop a general multiphase particle method based on the original MPS method and thus this work is of great significance for improving the numerical method for simulating the instability in reactor severe accident and two-phase flows in vessel. This research is intended to provide a study of the instability based on the MPS method. Latest achievements of mesh-free particle methods in instability are researched and a new multiphase MPS method, which is based on the original one, for simulating instability has been developed and validated. Based on referring to other researchers’ papers, the Pressure Poisson Equation (PPE), the viscosity term, the free surface particle determination part and the surface tension model are optimized or added. The numerical simulation on stratification behavior of two immiscible flows is carried out and results are analyzed after data processing. It is proved that the improved MPS method is more accurate than the original method in analysis of multiphase flows. In this paper, the main purposes are simulating and discussing Rayleigh-Taylor (R-T) instability and Kelvin-Helmholtz (K-H) instability. R-T and K-H instability play an important role in the mixing process of many layered flows. R-T instability occurs when a lower density fluid is supported by another density higher fluid or higher density fluid is accelerated by lower density fluid, and the resulting small perturbation increases and eventually forms turbulence. K-H instability is a small disturbance for two different densities, such as waves, at the interface of the two-phase fluid after giving a fixed acceleration in the fluid. Turbulence generated by R-T instability and K-H instability has an important effect in applications such as astrophysics, geophysics, and nuclear science.


2015 ◽  
Vol 137 (5) ◽  
Author(s):  
Kyung Sung Kim ◽  
Moo-Hyun Kim ◽  
Jong-Chun Park

For oil/gas production/processing platforms, multiple liquid layers can exist and their respective sloshing motions can also affect operational effectiveness or platform performance. To numerically simulate those problems, a new multiliquid moving particle simulation (MPS) method is developed. In particular, to better simulate the relevant physics, robust self-buoyancy model, interface searching model, and surface-tension model are developed. The developed multiliquid MPS method is validated by comparisons against experiment in which three-liquid-sloshing experiment and the corresponding linear potential theory are given. The validated multiliquid MPS program is subsequently coupled with a vessel-motion program in time domain to investigate their dynamic-coupling effects. In case of multiple liquid layers, there exists a variety of sloshing natural frequencies for respective interfaces, so the relevant physics can be much more complicated compared with the single-liquid-tank case. The simulation program can also reproduce the detailed small-scale interface phenomenon called Kelvin–Helmholtz instability. The numerical simulations also show that properly designed liquid cargo tank can also function as a beneficial antirolling device.


2009 ◽  
Vol 629 ◽  
pp. 231-262 ◽  
Author(s):  
ERIC JOHNSEN ◽  
TIM COLONIUS

A high-order accurate shock- and interface-capturing scheme is used to simulate the collapse of a gas bubble in water. In order to better understand the damage caused by collapsing bubbles, the dynamics of the shock-induced and Rayleigh collapse of a bubble near a planar rigid surface and in a free field are analysed. Collapse times, bubble displacements, interfacial velocities and surface pressures are quantified as a function of the pressure ratio driving the collapse and of the initial bubble stand-off distance from the wall; these quantities are compared to the available theory and experiments and show good agreement with the data for both the bubble dynamics and the propagation of the shock emitted upon the collapse. Non-spherical collapse involves the formation of a re-entrant jet directed towards the wall or in the direction of propagation of the incoming shock. In shock-induced collapse, very high jet velocities can be achieved, and the finite time for shock propagation through the bubble may be non-negligible compared to the collapse time for the pressure ratios of interest. Several types of shock waves are generated during the collapse, including precursor and water-hammer shocks that arise from the re-entrant jet formation and its impact upon the distal side of the bubble, respectively. The water-hammer shock can generate very high pressures on the wall, far exceeding those from the incident shock. The potential damage to the neighbouring surface is quantified by measuring the wall pressure. The range of stand-off distances and the surface area for which amplification of the incident shock due to bubble collapse occurs is determined.


Author(s):  
Ronghua Chen ◽  
Lie Chen ◽  
Wenxi Tian ◽  
Guanghui Su ◽  
Suizheng Qiu

In the typical boiling water reactor (BWR), each control rod guide tube supports four fuel assemblies via an orificed fuel support piece in which a channel is designed to be a potential corium relocation path from the core region to the lower head under severe accident conditions. In this study, the improved Moving Particle Semi-implicit (MPS) method was adopted to analyze the melt flow and ablation behavior in this region during a severe accident of BWR. A three-dimensional particle configuration was constructed for analyzing the melt flow behavior within the fuel support piece. Considering the symmetry of the fuel support piece, only one fourth of the fuel support was simulated. The eutectic reaction between Zr (the material of the corium) and stainless steel (the material of the fuel support piece) was taken into consideration. The typical melt flow and freezing behaviors within the fuel support piece were successfully reproduced by MPS method. In all the simulation cases, the melt discharged from the hole of the fuel support piece instead of plugging the fuel support piece. The results indicate that MPS method has the capacity to analyze the melt flow and solidification behavior in the fuel support piece.


Sign in / Sign up

Export Citation Format

Share Document