Non-Symmetrical Collapse of an Empty Cylindrical Cavity Studied with Smoothed Particle Hydrodynamics

: The non-symmetrical collapse of an empty cylindrical cavity is modeled using Smoothed Particle Hydrodynamics. The presence of a nearby surface produces an anisotropic pressure ﬁeld generating a high-velocity jet that hits the surface. The collapse follows a different dynamic based on the initial distance between the center of the cavity and the surface. When the distance is greater than the cavity radius (detached cavity) the surface is hit by traveling shock waves. When the distance is less than the cavity radius (attached cavity) the surface is directly hit by the jet and later by other shock waves generated in the last stages of the of the collapse. The results show that the surface is hit by a stronger shock when distance between the center of the cavity and the surface is zero while showing more complex double peaks behavior for other distances.


Introduction
Cavitation is a phenomenon occurring in a liquid that undergoes rapidly changes in pressure.At first a bubble, or cavity, nucleates and growth over a nucleation site.After growing up to maximum value, the cavity collapses generating shock waves [1].
When the pressure field in the liquid is isotropic, the cavity preserves its symmetry and the collapse has spherical symmetry [2][3][4].However, when an anisotropic pressure field drives the collapse, the cavity does not preserve its symmetry and folds in a specific direction generating a high-velocity flow known as jet [1,5].The anisotropic pressure field is generated by the so-called anisotropic driver, which also defines the jet folding direction [6].The most common anisotropic drivers are rigid or free surface, gravity, presence of neighbor bubbles or a combination of the above.Another mechanism that makes a cavity folds during the collapse is the interaction with a shock waves that folds the cavity in the shock direction.This situation is known as shock-induced collapse, while the pressure driven collapse is called Rayleigh collapse [7,8].
From an engineering point of view, the high-speed jet generated in the Rayleigh collapse by the presence of a nearby solid surface is the scenario most addressed in the literature [8,9].In fact, the high-speed flow can hit the surface causing erosion and eventually loss of material.This phenomenon, known as cavitation erosion, is undesirable in industrial, military and power station equipment such as pump impellers, high-speed propellers and turbine blades [10][11][12].
Thanks to its Lagrangian nature, Smoothed Particle Hydrodynamics (SPH) can handle problems with large deformations better than mesh-based technique [13] and, for this reason, it is particularly suited for studying cavitation.Joshi et al. [14] developed a Smoothed Particle Hydrodynamics axisymmetric solver to simulate a shock-induced collapse of an empty cavity and the erosion process of the nearby surface.Pineda et al. [15] used a SPH-ALE method to study a gas filled cylindrical cavity Rayleigh collapse far and near a surface.Nair and Tomar [16] used SPH to simulate an oscillating gas filled cylindrical cavity under a variable isotropic pressure field.Albano and Alexiadis [17] developed a SPH model to study the Rayleigh collapse of a gas filled cylindrical cavity that takes in account the heat diffusion at the gas-liquid interface.Among these studies, only Joshi et al. and Pineda et al. take in account anisotropic collapses.Joshi et al.only for the shock-induced collapse case, while Pineda at al only for cavity detached from the surface.Despite its importance in practical applications, non-symmetrical Rayleigh collapse of surface attached cavities, to the best of our knowledge, has never been investigated with the SPH method.Moreover, cavities commonly nucleate on surfaces rather than away from it [18][19][20].This work, therefore, develops a SPH model for a non-symmetrical Rayleigh collapse of a cylindrical cavity which is considered the greatest causes of erosion and material loss [8,9].
Thanks to its particle nature SPH is part of the Discrete Multi-Physics framework where, coupled with other particle methods [38][39][40], is used to address multi-complex physics phenomena [41][42][43][44][45] to overcome the single weakness of each method.
The idea behind the SPH formulation lies in the integral representation of any continuum function f (r) depending on the three-dimensional position vector r where W is the smoothing function or kernel and h is the smoothing length.Dividing the domain in a finite number of computational particles with their own mass m = ρdr 3 , is possible approximate any continuum function f (r) as where m i , ρ i and r i are mass, density and position of the i-th particle.Equation ( 2) is known as particle approximation in SPH literature [13].We used the SPH particle approximation to discretize continuity, momentum and energy conservation equation: where v is the velocity vector with v ij = v i − v j , e internal energy and Π ij is the artificial viscosity introduced by Monaghan [24].The Monaghan artificial viscosity depends on a constant parameter called dimensionless dissipation factor, α, and from the speed of sound of particle i and j, c i and c j , with the following relationship. (4) To solve the set of Equations shown before, an Equation Of State (EOS) that links pressure P and density ρ is required.

Computational Set Up
In this work, we focus on the non-symmetrical collapse induced by an anisotropic pressure field.The anisotropy is generated by non-symmetric water domain shown in Figure 1.The domain is divided in three concentric regions, delimited by three different radii.In each region, different types of computational particles are used: The green region in Figure 1 does not refer to a different type of particle.We highlight it because, later, the pressure generated during the collapse will be monitored in the green region.As explained in next section, when a non-symmetric collapse is studied, it is necessary to quantify the anisotropy of the collapse.When the anisotropy is induced by the presence of a nearby solid surface, is common to use the stand-off, γ [6] defined as where d is the distance between the cavity center and the wall (see Figure 1) and R 0 its initial radius.The dynamic of a cylindrical cavity in a Rayleigh collapse, where the driving force is the pressure difference between the pressure in the liquid, P ∞ = P L , and the pressure in the cavity, p b , is described by a 2D Rayleigh-Plesset (2DRP) equation [46]: However, the 2DRP equation only describe the dynamic of a symmetrical collapse.Form our knowledge an equation to validate the dynamic of a non-symmetrical collapse has still to be developed.For this reason, as commonly done in the literature [15,17,46,47], the model presented as been validated for the case of symmetric collapse against Equation (6) [17].
In the simulations we use the Lucy Kernel [13] and the Tait EOS [13] with a smoothing length of h = 1.3 • dL where dL is the initial particle spacing.The dimensionless dissipation factor, time step and speed of sound were set as α = 1, t s = 1e − 10 [s] and c 0 = 1484 [m s −1 ].The sensitivity of the results to parameters such as the kernel function, EOS or h, was investigated in a previous publication for the case of symmetric collapse [17].In the next section, the particle resolution for the specific case of non-symmetric collapse is discussed.

Software for Simulation, Visualization and Post-Process
The simulations were run with the open-source code simulator LAMMPS [48,49].The visualization and data post-processing were generated with the open-source code OVITO [50].

Results
As mentioned in the introduction, collapsing cavities formed during cavitation generate high-pressure shock wave [1].When the collapse occurs near a solid surface, the symmetry of the cavity is not preserved and this produces a high-velocity flow known as re-entrant jet [5,7,8].The characteristics of the jet depend on γ (see Equation ( 5)).When γ > 1, the jet hits the opposite side of the cavity generating a high-pressure shock wave that travels in the direction of the solid surface.When γ ≤ 1, the jet hits directly the solid surface generating a water hammer pressure [51,52] that causes erosion.

Preliminary Results
Four preliminary simulations with different particle resolutions have been carried out with R 0 = 100 [µm], P ∞ = 5 MPa and γ = 1.2.These values are chosen to match those commonly used in computational studies [14,17,52].Figure 2 shows the total collapse time and the max jet speed (defined as V = |v|) for the different resolutions.Particle resolution is specified as the ratio between the initial particle spacing dL and the initial cavity radius R 0 .Both the collapsing time and max jet velocity converge for dL/R 0 ≥ 40.This resolution value was also independently found by Joshi et al. and Pineda et al.At this resolution, the model correctly reproduces jet formation [5,7]. Figure 3 shows that as the collapse proceeds, wall proximity produces a pressure difference between the top and the bottom of the cavity.This pushes down the upper-side of the cavity inducing the formation of a high-velocity/low-pressure jet.
In Section 3.3 we show simulations at higher resolution.In fact, at higher resolutions, we will be able to uncover more details of the velocity pattern occurring during the collapse.

Stand-Off and Pressure Analysis
To study the effect of γ on the collapse, eight different stand-offs, in the range [0, 1.4], are simulated with P ∞ = 5 and 50 MPa.The upper limit γ = 1.4 was chosen since this study focuses on strongly deformed collapses, which are known to produce cavitation erosion [1,8].In Figure 4, the dimensionless collapsing time and the maximal pressure at the wall are plotted for different γ.The dimensionless time is defined as the ratio between t γ , the collapsing time obtained in the simulation for a given γ, and t ∞ , the collapsing time for γ → ∞ [17].The dimensionless pressure is defined as the ratio between the maximal pressure at the wall and the characteristic pressure of a Rayleigh collapse [52].The pressure at the wall is calculated at the green circle shown in Figure 1 positioned below the center of the cavity with a diameter of 0.01 mm.The collapsing time shows a maximum for γ = 1 and a minimum for γ = 0.When γ = 1, the lower part of the cavity touches the wall, while the upper part is free to move (Figure 5c).During the collapse, the cavity loses its symmetry generating a re-entrant jet, which travels for the whole diameter of the cavity before reaching the surface.When γ = 0 the cavity reduces to a semi sphere.Therefore, the collapse is symmetric again (Figure 5h) and, in fact, t γ /t ∞ ≈ 1 (Figure 4a) as for γ → ∞   The pressure is shown in Figure 4b.The normalized maximum wall pressure is thought to scale with γ −1 [1,6].However, this was reported for values of γ > 1.4, which are outside the range investigated in this study.Other studies show results similar to ours in the range of investigation [9,53].

Pressure and Speed Developed in High-Resolution Collapse
For 3 different stand-offs (0.0, 0.6 and 1.0), pressure histories at the center of the surface (i.e., green circle in Figure 1) are plotted in Figure 6.For γ = 0, the pressure shows a single maximum of roughly 1200 MPa.This pressure is generated by a water hammer impact of the collapsing cavity on the surface.Since the symmetry is preserved, see Section 3.2, the collapse ends when the cavity impacts on the surface generating a high-pressure shock wave (Figure 7a).After the impact, the pressure at wall decreases as the shock moves into the liquid (Figure 7b,c).
When γ = 0.6 and 1 the pressure shows a double peak.As with γ = 0, the first peak is generated by the water hammer impact of the cavity with the surface (Figures 8a and 9a).After the impact, the jet splits into two high-speed/low-pressure lateral jets (Figure 8b) that will later impact with the sides of the cavity generating collapsing circles and a pair of "side" pressure waves (Figures 8c and 9b).
The collapse ends when the fluid fills the circles.When this happens a third pair of pressure waves is generated (Figures 8d and 9c).The waves will later merge at the center resulting in a second pressure peak (Figures 8e and 9d). Figure 10 shows a schematic representation of a wall-attached cavity collapse with the three pairs shock formation for γ = 0.6.It can be difficult to clearly identify these hydrodynamics patterns if the simulation is run at lower resolution (e.g., Joshi et al. 2019).< l a t e x i t s h a 1 _ b a s e 6 4 = " n f m B q t d e s J 5 I J y l s E V g f r g S 5 t V c = " Side pressure waves < l a t e x i t s h a 1 _ b a s e 6 4 = " q 5 c U     As shown in Figure 6, when γ = 1 the first peak is higher compared to the first peak of γ = 0.6.Before hitting the surface, the jet travels a longer distance for γ = 1 rather than γ = 0.6.Therefore, it is accelerated by the pressure gradient for longer resulting in a higher water hammer pressure at the wall [8,51,54].
However, for the second peak, the behavior is reversed, and the peak is higher for γ = 0.6.The water between the collapsing rings and the wall acts as a shield mitigating the shock wave [26,55].When γ = 1 the distance between the rings and the wall is higher than γ = 0.6 resulting in a greater mitigation effect and lower pressure peak.

Conclusions
A SPH model is developed to study non-symmetrical Rayleigh collapse of an empty cylindrical cavity.When the cavity collapse near a solid surface the anisotropy of the pressure field induce the formation of the re-entrant jet.This anisotropy is quantified with the stand-off, γ.
Different collapses in the range 0 ≤ γ ≤ 1.4 are discussed showing that the model can correctly simulate the jet physics and the consequent pressure fields at the surface.
We study the collapse of attached cavities (γ ≤ 1), where the liquid exerts the maximum pressure over the surface.When γ = 0 the collapse is symmetric again and pressure shows a single max, which corresponds to the end cavity collapse.When 0 < γ ≤ 1 the

• 0 Figure 1 .Figure 1 .
Figure 1.Geometry of the simulation box Figure 1.Geometry of the simulation box.

Figure 2 .Figure 2 .
Figure 2. Resolution convergences for collapse time (a) and max jet speed (b) Figure 2. Resolution convergences for collapse time (a) and max jet speed (b). d

Figure 5 .
Figure 5. Cavity shapes for different stand off (dL/R 0 = 40 and P ∞ = 5 MPa).Each shape shows the cavity outline at a different time.

Figure 5 .
Figure 5. Cavity shapes for different stand-offs (dL/R 0 = 40 and P ∞ = 5 MPa).Each shape shows the cavity outline at a different time.
[mm] < l a t e x i t s h a 1 _ b a s e 6 4 = " p 1 / 7 b 6 Q I K b u c f R 5 I C c T I 7 I h / z o I = " > A A A B 9 n i c b V B N S 8 N A F H y p X 7 V + V T 1 6 W S y C p 5 I U Q Y 9 F L x 4 r W C u 0 s W y 2 L + 3 S z S b s b p Q S + j + 8 e F A Q r / 4 W b / 4 b N 2 0 O 2 j q w M M y 8 4 b 2 d I B F c G 9 f 9 d k o r q 2 v r G + X N y t b 2 z u 5 e d f / g T s e p Y t h m s Y j V f U A 1 C i 6 x b b g R e J 8 o p F E g s B O M r 3 K / 8 4 h K 8 1 j e m k m C f k S H k o e c U W O l h 5 Z C r V O F J A + N + 9 W a W 3 d n I M v E K 0 g N C r T 6 1 a / e I G Z p h N I w Q b X u e m 5 i / I w q w 4 h l P w 4 B y a c A 0 t a A M D B c / w C m / O k / P i v D s f 8 9 G S U 2 Q O 4 Q + c z x 8 4 P Z L J < / l a t e x i t > Pressure preak t e x i t s h a 1 _ b a s e 6 4 = " 2 P l z + 2 0 y

Figure 7 .Figure 7 .
Figure 7. Pressure and velocity field in the liquid and the spacial pressure trend at the wall (dL/R 0 =133, γ = 0 and P ∞ = 50 MPa)Figure 7. Pressure and velocity field in the liquid and the spatial pressure trend at the wall (dL/R 0 = 133, γ = 0 and P ∞ = 50 MPa).
r d I n Y a L s i w 2 Z q 7 8 T E x Z p P Y 4 C O x k x M 9 T L 3 k z 8 z + t k J r z 0 J y J O M 4 M x X y w K M 0 l Q m Y 9 L 5 Q y I 0 c W 8 K 4 E v Z W w o d M M W 5 s W y V b A l 3 + 8 i p p 1 a r 0 v E r p b a 1 S v 8 r 7 K M I x n M A Z U L i A O t x A A 5 r A 4 R G e 4 R X e n C f n x X l 3 P h a j B S f P H M E f O J 8 / i v S T / g = = < / l a t e x i t > 1st pressure peak < l a t e x i t s h a 1 _ b a s e 6 4 = " O n g+ z h c V N D 7 O Y d H 5 Y V H k E / v R v p 4 = " > A A A B 7 3 i c d V B N S 8 N A E J 3 4 W e t X 1 a O X x S J 4 C k m 1 t d 6 K X r w I F e w H t K F s t p t 2 7 W Y T d z d C C f 0 T X j w o 4 t W / 4 8 1 / 4 z a t o K I P B h 7 v z T A z z 4 8 5 U 9 p x P q y F x a X l l d X c W n 5 9 Y 3 N r u 7 C z 2 1 R R I g l t k I h H s u 1 j R T k T t K G Z 5 r Q d S 4 p D n 9 O W P 7 q Y + q 1 7 K h W L x I 0 e x 9 Q L 8 U C w g B G s j d S + w k y g W 6 p 7 h a J j n 1 W c k 3 I F Z a T i l G b k u F p G r u 1 k K M I c 9 V 7 h v d u P S B J S o Q n H S n V c J 9 Z e i q V m h N N J v p s o G m M y w g P a M V T g k C o v z e 6 d o E O j 9 F E Q S V N C o 0 z 9 P p H i U K l x 6 J v O E O u h + u 1 N x b + 8 T q K D q p c y E S e a C j J b F C Q c 6 Q h N n 0 d 9 J i n R f G w I J p K Z W x E Z Y o m J N h H l T Q h f n 6 L / S b N k u 2 X b u S 4 V a + f z O H K w D w d w B C 6 c Q g 0 u o Q 4 N I M D h A Z 7 g 2 b q z H q 0 X6 3 X W u m D N Z / b g B 6 y 3 T x 7 X k A g = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " d C m G H y 3 n y 8 O 1 u S E J o t P j T L L N Y U I = " > A A A B 9 H i c j V D L S g N B E O y N r x h f U Y 9 e B o P g K e w G Q Y 9 B L x 4 8 R D A m k C x h d t K b j J m d X W Z 6 h b D k N 7 x 4 U B C v f o w 3 / 8 b N 4 6 C i Y E F D U d V N F x U k S l p y 3 Q + n s L S 8 s r p W X C 9 t b G 5 t 7 5 R 3 9 2 5 t n B q B l 8 t O I u b f f g G 5 + 0 T e N O R x w = = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " n f m B q t d e s J 5I J y l s E V g f r g S 5 t V c = " > A A A B / X i c b V B N S 8 N A F N z U r 1 q / o j 1 6 W S y C p 5 I U Q Y 9 F L x 4 r W l t o Q 9 l s X t q l m 0 3 Y 3 V R C q H / F i w cF 8 e r / 8 O a / c d v m o K 0 D C 8 P M G 9 7 b 8 R P O l H a c b 6 u 0 t r 6 x u V X e r u z s 7 u 0 f 2 I d H D y / l a t e x i t > Side pressure waves < l a t e x i t s h a 1 _ b a s e 6 4 = " q 5 c U B p y m k f b B 6 J T b Q n Y Q T 6 p Y K E Y = " > A A A B + n i c d V D N S 8 M w H E 3 9 n P N r 0 6 O X 4 B A 8 l X S w u d 2 G u 3 i c 4 D 5 g K y P N 0 i 0 s T U u S K q C + P o p / J / 0 q r Z T s 9 F 1 t d K 6 z O M o g B N w C s 6 B A y 5 A C 1 y B D u g C A u 7 A A 3 g C z 9 a 9 9 W i 9 W K + r 1 j U r n z k G P 2 C 9 f Q L H a 5 R a < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " X U 6 h 6 f f X 4 E N 1 d q B F n L Y M 5 2 + a x X w = " > A A A B + X i c d V D L S s N A F L 2 p r 1 p f U Z d u B o v g K i Q V 0 y y L 3 b i s Y B / Q h j K Z T t q h k w c z k 0 I J / R M 3 L h R x 6 5 + 4 8 2 + c p h V U 9 M D A 4 Z z 7 m h O k n E l l 2 x 9 G a W N z a 3 u n v F v Z 2 z 8 4 P D K P T z o y y Q S h b Z L w R P Q C L C l n M W 0 r p j j t p Y L i K O C 0 G 0 y b S 7 8 7 o 0 K y J L 5 X 8 5 T 6 E R 7 H L G Q E K y 0 N T b P J B O F U I j 2 K 4 1 T S o V m 1 L a / u u W 4 d F c S 7 c l b E r b n I s e w C V V i j N T T f B 6 O E Z B G N F e F Y y r 5 j p 8 r P s V B M D 1 5 U B p m k K S Z T P K Z 9 T W M c U e n n x e U L d K G V E Q o T o V + s U K F + 7 8 h x J O U 8 C n R l h N V E / v a W 4 l 9 e P 1 O h 5 + c s T j N F Y 7 J a F G Y c q Q Q t Y 0 A j J i h R f K 4 J J o L p W x G Z Y I G J 0 m F V d A h f P 0 X / k 0 7 N c q 4 t + 6 5 W b d y s 4 y j D G Z z D J T h Q h w b c Q g v a Q G A G D / A E z 0 Z u P B o v x u u q t G S s e 0 7 h B 4 y 3 T x N P k / c = < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " V D + + + p S J J Y N V p Q 0 J j h C / c w K L + r g = " > A A A B / H i c b V B N S 8 N A F H z x s 9 a v W o 9 e F o v g q S R F 0 G P R i 8 c K 1 h b a U D a b 1 3 b p Z h N 2 N 2 I J / S t e P C i I V 3 + I N / + N m z Y H b R 1 Y G G b e 8 N 5 O k A i u j e t + O 2 v r G 5 t b 2 6 W d 8 u 7 e / s F h 5 a j 6 o O N U M W y z W M S q G 1 C N g k t s G 2 4 E d h O F N A o E d o L J T e 5 3 H l F p H s t 7 M 0 3

Figure 8 .Figure 8 .
Figure 8. Pressure and velocity field in the liquid and the spacial pressure trend at the wall (dL/R 0 =133, γ = 0.6 and P ∞ = 50 MPa) Figure 8. Pressure and velocity field in the and the spatial pressure trend at wall (dL/R 0 = 133, γ = 0.6 and P = 50 MPa).
t e x i t s h a 1 _ b a s e 6 4 = " U l A E y U 2 s V C S S 1 S I 7 0 Q C T 9 m h p C N A = " > A A A B + 3 i c b V D L S g M x F L 1 T X 7 W + p r p 0 E y y C q z I p g i 6 L b l x W s L b Q D i W T 3 m l D M w + S j F J q P 8 W N C w V x 6 4 + 4 8 2 9 M 2 1 l o 6 4 H A 4 Z x 7 u D c n S K X Q x v O + n c L a + s b m V n G 7 t L O 7 t 3 / g l g / v d Z I p j k 2 e y E S 1 A 6 w o d M M W 5 s W y V b A l 3 + 8 i p p 1 a r 0 v E r p b a 1 S v 8 r 7 K M I x n M A Z U L i A O t x A A 5 r A 4 R G e 4 R X e n C f n x X l 3 P h a j B f P M E f O J 8 / i v S T / g = = < / l a t e x i t >

Figure 9 .Figure 9 .
Figure 9. Pressure and velocity field in the liquid and the spacial pressure trend at the wall (dL/R 0 =133, γ = 1 and P ∞ = 50 MPa)Figure 9. Pressure and velocity field in the liquid and the spatial pressure trend at the wall (dL/R 0 = 133, γ = 1 and P ∞ = 50 MPa).
t e x i t s h a 1 _ b a s e 6 4 = " 1 f 6u 7 z x v H C J 2 X r G 9 f A w 9 G S R H r M k = " > A A A B 9 X i c d V D L S s N A F J 3 U V 6 2 v q k s 3 g 0 V w F Z I Y 2 r o r u n F Z 0 d p C G 8 p k c t M O n T y c m R R K 6 X e 4 c a E g b v 0 X d / 6 N k 7 a C i h 4 Y O J x z D / f O 8 V P O p L K s D 6 O w s r q 2 v l H c L G 1 t 7 + z u l f c P 7 m S S C Q o t m v B Ed H w i g b M Y W o o p D p 1 U A I l 8 D m 1 / d J n 7 7 T E I y Z L 4 V k 1 S 8 C I y i F n I K F F a 8 m 5 Y A F g H p M w E 9 M s V y z y r 3 V N I l f P 0 U / 0 / a j m m 7 p m 1 f O 5 X G x b K P I j p C x + g U 2 a i G G u g K N V E L U X S P H t A T e j b G x q P x Y r w u R g v G M n O I f s B 4 + w T F I Z K T < / l a t e x i t > Side pressure < l a t e x i t s h a 1 _ b a s e 6 4 = " d Q J 3 y u P V d 0 Z A y S o I r o O d H a A I x w c = " > A A A B 7 X i c d V B N S 8 N A E N 3 4 W e t X 1 a O X x S J 4 C t k Y 2 n o r e v F Y w d h C G 8 p m u 2 m X b j Z h d 1 M p o b / B i w c F 8 e r / 8 e a / c d N W U N E H A 4 / 3 Z p i Z F 6 a c K e 0 4 H 9 b K 6 t r 6 x m Z p q 7 y 9 s 7 u 3 X z k 4 v F N J J g n 1 S c 2 o p w J 6 m u m O e 2 k k u I 4 5 LQ d j q 8 K v z 2 h U r F E 3 O p p S o M Y D w W L G M H a S P 4 9 n l D V r 1 Q d + 6 J R c 7 0 a d G z H q S M X F c S t e + c e R E Y p U A V L t P q V 9 9 4 g I V l M h S Y c K 9 V F T q q D H E v N C K e z c i 9 T N M V k j I e 0 a 6 j A M V V B P j 9 2 B k + N M o B R I k 0 J D e f q 9 4 k c x 0 p N 4 9 B 0 x l i P 1 G + v E P / y u p m O G k H O R J p p K s h i U Z R x q B N Y f A 4 H T F K i + d Q Q T C Q z t 0 I y w h I T b f I p m x C + P o X / k 7 Z r I 8 9 G 6 M a t N i + X e Z T A M T g B Z w C B O m i C a 9 A C P i C A g Q f w B J 4 t Y T1 a L 9 b r o n X F W s 4 c g R + w 3 j 4 B A Z y P S w = = < / l a t e x i t > waves < l a t e x i t s h a 1 _ b a s e 6 4 = " k g 7 e o N d N 2 H E s q E l s g H e D + E V i m 2 w = " > A A A B + 3 i c d V D L S g M x F M 3 U V 6 2 v q S 7 d B I v g q k z G o a 2 7 o h u X F a w t t E P J p L d t a O Z B k l F K 7 a e 4 c a E g b v 0 R d / 6 N m b a C i h 4 I H M 6 5 h 3 t z g k R w p R 3 n w 8 q t r K 6 t b +

Figure 10 .
Figure 10.Schematic representation of a wall-attached cavity collapse (γ = 0.6): blue surface represents the cavity, black arrows represent the liquid flow, red circles represent the shock waves propagation at different collapse stages.

Figure 10 .
Figure 10.Schematic representation of a wall-attached cavity collapse with the three pairs shock formation for γ = 0.6.The blue lines represents the cavity, black arrows represent the liquid flow, red circles represent the shock waves propagation at different collapse stages.